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L’indium utilisé dans la cathode creuse était de 
l’indium naturel contenant 96% de 115]n et 4 % 
de 3]n ; on n’a pas trouvé d’indication de compo- 
santes résolues dues à l’isotope le plus rare. 


La source lumineuse. — La figure 1 montre le 


dessin du tube à décharge à cathode creuse. La 
partie extérieure du tube est la cathode et est 


Fic: 1. — Cathode creuse refroidie à l’eau. 


constituée par un tube d’acier de 8 em de long et 
3,5 em de diamètre. L’extrémité inférieure du tube 
est fermée par une pièce de base plate au centre 
de laquelle se trouve un petit cylindre de 15 mm 
de long et 5 mm de diamètre qui constitue la 
cavité cathodique. Un tube de verre (ou silice 
. fondue) de 10 cm de long s’adapte dans la cathode 
en acier, reposant sur la base ; l’extrémité supé- 
rieure du tube est fermée par une fenêtre plane 
scellée. À cette extrémité supérieure du tube de 
silice sont fixés deux tubes latéraux ; l’un est 


connecté à un tube contenant du charbon activé, 
l’autre sert au passage d’un fil fin pour la con- 
nexion de l’anode. Cette anode est constituée par 
un cylindre creux d’aluminium s’adaptant exac- 
tement dans le tube de silice et se terminant à 1 em 
au-dessus de l’extrémité inférieure de la cathode. 
L'ensemble est rendu étanche au moyen d’un joint 
de picéine entre l'extrémité supérieure du tube 
d’acier et le tube de verre ; le passage du fil d’anode 
est également scellé à la picéine. 

La lampe.à cathode creuse est remplie d’un gaz 
porteur (néon ou hélium) par le tube latéral, qui est 
connecté au charbon activé. La circulation du gaz 
porteur n’est pas nécessaire ; le gaz est purifié et 
maintenu dans l’état de pureté nécessaire par le 
charbon refroidi avec de l’azote liquide. Après 
réglage de la pression du gaz porteur à sa valeur 
optimum, on peut fermer le robinet et la lampe 
peut être séparée du système de pompage et de 
remplissage. Avec un bon robinet et des joints de 
picéine corrects, la lampe peut fonctionner des 
dizaines d'heures sur une période de plusieurs mois 
sans nouveau remplissage: 

Environ 100 mg d’indium sont placés dans la 
cavité cathodique ; avec du néon comme gaz 
porteur, les raies de résonance du spectre d’arc de 
l’indium sont émises avec une bonne intensité pour 
des courants de quelques milliampères, tandis que 
le multiplet du premier spectre d’étincelle apparaît 
fortement au-dessus de 50 mA. L’intensité de ce 
multiplet dépend de la pression du néon, mais la 
pression optimum n’est pas critique; avec un 
courant de 100 mA la pression la plus favorable est 
d'environ 5 mm de mercure. Le tube est refroidi au 
moyen d’eau glacée, dans un grand récipient, la 
cathode plongeant de 5 cm dans l’eau. Avec des 
courants inférieurs à 50 mA les raies de réso- 
nance 4 101 et 4 511 À sont exemptes de renver- 
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sement et d’élargissement par auto-absorption ; les 
composantes hyperfines de 4511 À, distantes 
de 37 mK, sont bien résolues, ce qui implique une 
largeur de raie n’excédant pas 25 à 30 mK ; la 
largeur Doppler de ces raies, à 300 °K est approxi- 
mativement 25 mK. Aux courants plus élevés, 
nécessaires pour l’émission des raies d’étincelle, ces 
raies montrent, comme il se doit, un élargissement 
par auto-absorption. La largeur des raies du mul- 
tiplet 4 638-4 685 À est nettement plus grande que 
celle des raies de résonance ; les raies les plus 
proches qui peuvent être clairement séparées, même 
avec un courant limité à 50 mK sont distantes 
d'environ 60 mK, ce qui indique une largeur de 40 
à 50 mK, soit à peu près deux fois la largeur 
Doppler. Cet élargissement peut être dû à l'effet 
Stark produit par le champ électrique des ions 
d’indium voisins ; mais alors il faut admettre que 
les raies d’étincelle sont émises par une région de la 
cavité cathodique différente de celle qui émet les 
raies d’arc, puisque ces dernières sont exemptes de 
cet élargissement. Le courant normalement utilisé 
est 100 mA et les raies d’étincelle ne sont pas nota- 
blement plus larges qu'aux courants plus faibles. 
Pour l'observation de la raie faible 4 652 À le 
courant est poussé jusqu’à 140 mA. Des courants 
plus élevés sont à éviter pour les poses longues à 


cause de la pulvérisation de l’indium qui forme une. 


couche conductrice de métal sur la paroi interne du 
tube de verre, produisant un court-circuit entre la 
cathode et l’anode. 


L’interféromètre et le spectrographe. — C’est la 
méthode photographique qui a été utilisée. Le sys- 
tème optique était constitué par un spectrographe 
produisant la première séparation des raies, soit un 
instrument à transmission à deux prismes de courte 
distance focale et de grande ouverture relative, soit 
un instrument Littrow de grande distance focale 
et de haute dispersion, puis un étalon Fabry-Perot 
à montage externe pour donner la haute résclution. 
Le spectregraphe à transmission a l’avantage, grâce 
à sa grande ouverture relative, d'utiliser toute la 
lumière transmise par l’étalon, même quand celui-ci 
est utilisé à pleine ouverture, ce qui est essentiel 
dans le cas, par exemple, de la raie faible 4 652 À 
qui exige un temps de pose d’une heure et demie 
avec l’étalon travaillant à f/15. Cependant, en 
raison de sa dispersion limitée (25 À/mm), il ne 
pouvait être utilisé pour observer la raie 4 657 À 
de faible intensité relative par rapport à la raie 
voisine 4 656 À, puisque même avec une fente assez 
fine pour séparer les deux raies, les franges de la 
raie la plus faible étaient noyées dans le halo dû à la 
raie intense. Dans ce cas le spectrographe Littrow, 
ayant une dispersion cinq fois plus grande était 
utilisé ; l'ouverture relative 1/35 était suffisante 
pour permettre de photographier toutes les raies 
du multiplet 5 D — 4 3F à l'exception de la raie 
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interdite 4 652 À ; elle ne permettait pas non plus 
de photographier la raie d’intercombinaison 
4 674 À ;: 55D,—41F,. 

L’étalon était du type classique, avec des lames 
et des cales en silice fondue. Il était éclairé à 
travers un collimateur, sur la fente duquel était 
projeté l’image de la partie brillante de la cavité 
cathodique ; un objectif achromatique de 50 em de 
focale était utilisé pour former l’image du système 
de franges sur la fente du spectrographe. Les lames 
de l’étalon étaient recouvertes de 7 couches diélec- 
triques alternées de cryolite et de sulfure de zine ; 
à pleine ouverture la limite de résolution était 
approximativement un trentième d'ordre, et 
pouvait être poussée jusqu'à un cinquantième 
d'ordre en diaphragmant. Le contraste était suf- 
fisant pour observer sans difficulté une composante 
d'intensité 500 fois plus faible que l'intensité totale 
dans la raie 4 681 À. Un certain nombre de dis- 
tances des lames furent utilisées, comprises entre 
0,812 mm et 4,122 mm ; la distance des lames était 
déterminée approximativement en mesurant l’épais- 
seur des cales avec un palmer puis à partir de cette 
valeur, correcte à 3 ou 4 u près, une valeur précise 
était obtenue interférométriquement par mesure 
des diamètres des anneaux donnés par les compo- 
santes principales des raies vertes, bleu-vertes et 
bleues du cadmium. L’inclinaison de l’étalon était 
réglée de façon à obtenir dans le spectre de franges 
une dispersion suflisante pour ne pas introduire de 
perte de résolution, puisque dans l'observation 
photographique, si la dispersion n’est pas adéquate, 
la limite de résclution est déterminée non par le 
pouvoir de résolution mais plutôt par la dispersion. 
Une dispersion de 4 mm par ordre a été trouvée 
convenable, correspondant à environ 0,1 mm de 
séparation entre deux composantes à la limite de 
résolution ; si ce critère était appliqué à un réseau 
concave de 15 cm utilisé en montage Rowland 
normal, cela nécessiterait un rayon de courbure 
de 30 mètres, ou le double dans un montage 
Wadsworth. La supériorité de résclution de l’étalon 
par rapport à celle du réseau utilisé par Paschen et ! 
Campbell était due surtout à la plus grande dis- 
persion, soit environ 5 mm par À avec les plus 
petites cales et plus du double avec les plus grandes. 

Pour toutes les raies, au moins deux valeurs de 
l’épaisseur de l’étalon ont été utilisées, afin de 
s'assurer que les composantes observées étaient 
attribuées à l’ordre d’interférence correct. Dans le 
cas des raies 4 685 et 4 657 À, la largeur occupée 
par la structure hyperfine était, sauf pour les plus 
courtes épaisseur de l’étalon, supérieure à l’inter- 
valle spectral de l’étalon, et il était nécessaire de 
travailler avec la fente du spectrographe Littrow 
suffisamment fine pour que la dispersion sépare les 
composantes se superposant dans les ordres d’inter- 
férence adjacents. Avec une épaisseur de 1,66 mm, 
la fente devait être réduite à environ 0,6 mm; 
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mais, bien que l’aspect du système de franges soit, 
dans ces conditions, très déplaisant, il n’y avait pas 
de difficulté spéciale à mesurer les positions des 
composantes, et il n’y avait pas d'augmentation 
significative des erreurs statistiques. La déviation 
moyenne des mesures simples des composantes 


… fortes, bien résolues de leurs voisines était d’envi- 


ron 5 mK avec une épaisseur d’étalon de 2,5 mm, 
et 10 mK avec la petite épaisseur, 1,1 mm ; et 
puisque au moins 9 mesures différentes étaient 
faites, l'erreur statistique était réduite à 2 mK. 
Pour les composantes faibles, voisines des compo- 
santes fortes, la déviation moyenne était à peu près 
le double, et pouvait être associée avec des erreurs 
systématiques allant jusqu’à 10 mK avec la plus 
grande épaisseur d’étalon ou 20 mK avec la plus 
petite. La comparaison des valeurs trouvées avec 
les différentes épaisseurs fournit une bonne esti- 
mation de la précision; comme on peut le voir 
d’après les tableaux ci-dessous, dans les cas favo- 


rables les valeurs sont cohérentes à 2 ou 3 mK près, 


mais dans les cas défavorables — composantes 
faibles voisines de composantes intenses, ou 
mélanges non résolus — il y a des différences de 
l’ordre de 10 mK. 

Les nombres décrivant les intensités relatives 
des composantes sont des estimations visuelles, avec 
une échelle linéaire; Dans un intervalle de 10: 1 elles 
sont obtenues par comparaisons avec des rapports 
étalons donnés par un échelon absorbant. Dans 
l’intervalle de 10 à 5 elles sont probablement 
valables à un cinquantième de l’intensité donnée 
près ; en dessous de cela, la précision relative 
décroît et tombe à environ 50 % pour l'intensité 
la plus faible 1. Pour estimer l'intensité des compo- 
santes d'intensité plus faible que 1 : 10, des poses 
successives ont été faites pour lesquelle l'ouverture 
de l’étalon était réduite dans le rapport de 10à 1 et, 
pour la raie 4 681 À, aussi de 100 à 1. De la sorte, 
on peut avoir une estimation de l’ordre de grandeur 
de l’intensité relative des composantes très faibles. 


Tableau des mesures des structures hyperfines. 
— Le Tableau I montre les raies observées et 


… donne les indications relatives à leurs intensités 


relatives (totales), à la largeur totale de la structure 
hyperfine, au nombre de composantes observé 
réellement et calculé, dans le présent travail, et 
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dans celui de Paschen et Campbell. Dans tous les 
cas, le nombre de composantes observées est en 
dessous du nombre calculé ; cela peut provenir de 
composantes trop faibles pour être observées, ou du 
fait que deux ou plusieurs composantes appa- 
raissent comme un mélange sielles sont trop voisines 
pour être résolues. Seules les composantes qui ne 
sont pas de tels mélanges sont convenables pour en 
déduire la structure des termes : les nombres entre 
parenthèses donnés à droite des nombres de compo- 
santes observées montrent le nombre de compo- 
santes qui ne présentent pas ces mélanges. À cet 
égard le progrès réalisé par les nouvelles mesures 
est notable, le nombre total des composantes non 
mélangées étant passé de 36 à 62. 

Les tableaux suivants (tableaux IT)contiennent 
les résultats des mesures de structure hyperfine des 
raies précédentes. Les valeurs moyennes des 
mesures faites avec différentes épaisseurs d’étalons 
sont données séparément ainsi qu'une moyenne 
générale, obtenue d’après les valeurs trouvées avec 
les différentes épaisseurs, affectées de poids propor- 
tionnels au nombre de mesures faites et inver- 
sement proportionnels à la déviation moyenne. Ces 
valeurs sont données dans l’avant-dernière et la 
dernière horizontale des tableaux, la déviation 
moyenne étant la moyenne des différentes valeurs 
pour toutes les composantes mesurées. 

Les différentes composantes sont identifiées par 
des lettres dans la première colonne verticale de 
chaque tableau, pour faciliter la comparaison des 
nombres donnés dans chaque tableau avec les 
reproductions des courbes microphotométriques. 
On peut voir sur les interférogrammes que les 
composantes marquées de la lettre æ dans les 
courbes relatives aux raies 4 656 et 4 652 À sont 
déplacées par rapport au reste des composantes 
dans la direction de la dispersion du spectrographe 
et ne sont donc pas des composantes des raies en 
question mais sont dues à d’autresraies très voisines 
(la plus forte composante de la raie 4 657 À dans 
le premier cas, et une raié du néon dans le dernier 
cas). De même les deux composantes a et o de ia 
raie 4 685 qui se superposent sur les courbes micro- 
photométriques sont séparées par la dispersion du 
spectrographe et peuvent être mesurées séparément 
au comparateur. 

Les valeurs trouvées par Paschen et Campbell 


TABLEAU I 
RaïE INTENSITÉ LARGEUR NOMBRE DE COMPOSANTES 
DE LA S. H. F. CaLc. Ogs. (J.) Oss. (P. et C.) 
4.685 À 55D,—45%F, 20 4 313 mK 19 45 (14) 45 (14) 
4681 À 55D;—45F, 4 000 1 667 mK 21 13e (0) & (0) 
MI657 A5 SD, — 4 5F, 50 3 574 mK 13 42 (12) 49 + (12) 
& 656 À 535D,—45F, 300 2 090 mK 15 44 (9) 6) 
& 652 À 53D,— 45F, 1 /2 1415 mK 45 LOMME (7)) 20010) 
4 638 À 5D,—45F, 300 4 667 mK 9 get L(7) Ja) 
& 674 À 55D;—41F, 5 2 247 mK 19 9 (4) 5 (0) 
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TABLEAU II 
x 4 685 À 55D,—45%F, 


COMPOSANTE 081241851662 Moyenne * INTENSITÉ P.et"G 
a 0 0 0 0 10 0 10 
b 310 303 287 295 2 DO 
C 870 863 840 854 k 833 4 
d TAS0 41406489 1:450 3 4474 4 
e 1380 4403 1 379 10357 3 1 467 2 
Î 4 870 1899 1 889 1 893 EE 1 909 6 
£ se 2412 5 02H95 24925 LE 2 065 d 
h 2270002 -306-22809 2 296 4 2315 2 
] 260027064025 0 2 756 5 2 74h 5 
k 3020 3049 3 038 3 042 k 301 
l 3420 3468 3431 3 447 5 2120710000 
m 3,9001031628 281599 3 609 3 1002 
3 950 3 950 3 
n 3 9700231993 : 3 967 6 
4 040 4 040 2 
0 4290 4297 4290 4 293 3 4 296 2 
k B42 À 
Nombre k 40 12 
D. M. 20 40 10 
* Seulement 1,135 et 1,662. 
À 4 681 À 55D,—435F, 
COMPOSANTE 1,662 2,461 RC) DA MOYENNE INTENSITÉ P:retics 
a — 620 — 624 as — 624 1/8 
b NAN sa 2520 1/4 
C — 82 — 81 — 81 Fe — 81 4/2 LARENEE 
d 0 0 0 0 0 1 1/2 30 2 
Qu ra ERA AE AU 421 149 1420 120 4/2 Me : 
Î 225 226 223 225 3 
£ 9311 288 302 300 d Sr d 
h AAA 405 410 409 3 SCO 
ïl 538 541 037 539 2) RE: 
k 672 676 674 674 20 670 20 
l 849 846 850 848 5 AE 5e 
TS MATE E TRS 968 975 969 971 6 
99517 
n RASE 4 041 1048 1 040 1 043 7 desire 
Nombre 3 10 6 8 
D. M. 20 5 3 3 


sont ainsi données à titre de comparaison ; leurs 
intensités ont été normalisées à notre échelle déci- 
male. . 

Il est évident, en comparant les résultats obtenus 
avec différentes épaisseurs d’étalon que les mesures 
faites sur les raies 4 681, 4 656 et 4 638 À sont 
nettement plus exactes que les autres. Ces raies 
sont plus facilement mesurables car le domaine 
spectral couvert est plus petit et un étalon plus 
épais peut être utilisé ; de plus ce sont des raies 
fortes et on peut élargir la fente du spectrographe 
sans risque d’être gêr.é par d’autres raies. On peut 
se fier aux valeurs des moyennes à 3 mK près. Ces 
mesures de Paschen et Campbell ont été données 
à 1 mK près : mais le dernier chiffre n’a proba- 
blement aucun sens (les valeurs sont obtenues à 
partir des mesures exprimées à { mÀ près, ce qui 
est équivalent à 5 mK). Leurs composantes séparées 
de moins que 150 mK environ n’apparaissent pas 


résolues, et leurs résultats pour les raies 4 681 et 
4 652 ont peu de valeur. Pour les raies plus lar- 
gement espacées, où.les plus petites séparations 
sont de 200 à 300 mK, les mesures de Paschen et 
Campbell sont beaucoup plus satisfaisantes, et sont 
probablement bonnes à + 10 ou + 20 mK près. 
Pour ces raies, les nouvelles mesures, en raison des 
plus faibles épaisseurs d’étalon et des plus faibles 
largeurs de fente de spectrographe utilisées sont 
nettement moins exactes que pour les raies 4 681, 
4 652 et 4 638 À ; elles sont probablement bonnes 
à + 10 mK près. 

Si on tient compte du fait que leur pouvoir 
résolvant était plus faible, les mesures de Paschen 
et Campbell sont en bon accord avec les nouvelles 
mesures ; mais il y a des différences importantes 
sur deux points. Pour la raie 4 685 À, la compo- 
sante à 4412 mK observée par Paschen et 
Campbell est certainement absente dans les spectres 
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COMPOSANTE 1,135 1,662 Moyenne 

a 0 0 0 8 0 10 

b 377 379 378 8 362 9 

C 682 679 681 A 687 2 

d 1048 1071 4 059 2 4 045 & 

e 1290203 1 396 10 1380 10 

Î 1 690 4703 1 696 k 1769940413 

£ 2304102:523 2 314 9 21292 8 

k 2709202602 21097 3 2 607 A 

] 2 869 2 879 2 874 44/2 NN here 

k 3 094 3091 3 092 5 3 078 6 

l 3349 3 356 3 352 2 313072 

m 3 07Æ0 3 574 3 574 5 SD0 ND 

Nombre 42 42 
D. M. 10 10 
MIE GOOLA SD, 54 SF, 
COMPOSANTE 4,662 2,462 2,606 MOYENNE INTENSITÉ Pret 
a 0 0e 0 10 ŒuA0 
b 809 302 304 305 3 219802 
ce 690 686 685 688 7 GHITRNAT 
d 922 920 923 922 5 941 2 
e 1 080 1089 1 087 1 085 4 ser lee 
f 1253401495 27246205 1 254 5 49251 6 
£ 1426 1420 1 422 4 422 9 à AS 0 Nm 
k 1069224522 400152 1523 A RUES 
] 1447412700 %4,709 4 710 3 HER 
k 1 828 1835 1 818 1 824 6 ALES 
l 2 090 2 099 2 096 3 2507/6105 
Nombre 10 5 16 
D. M. 6] & A 


obtenus à l’interféromètre ; et on voit difficilement 
comment elle aurait pu être séparée dans les 


Fic. 2. — à 4 685 À. Épaisseur de l’étalon: 1,135 mm. 


spectres donnés par le réseau, puisqu'elle est 
donnée comme étant séparée de la composante 
voisine de 116 mK, ce qui est trop petit pour la 


résolution du réseau, mais facile à résoudre sur 
l'interférogramme. L’autre différence importante 
concerne la raie 4 638 À. D’abord, les composantes 
bien résolues a et b ne sont pas résolues dans le 
spectre de réseau, ensuite la composante d, qui est 
séparée de sa voisine de 200 mK environ manque 
dans le spectre de réseau et enfin une composante 
supplémentaire est donnée à 946 mK, c’est-à-dire 
seulement à 81 mK de sa voisine 1 027 mK, et 
ainsi trop proche pour être résolue dans le spectre 
de réseau. 

La reproduction des courbes microphoto- 
métriques des structures hyperfines est intéres- 
sante, car elle montre le caractère satisfaisant de la 
résolution des composantes. Dans le cas de la raie 
4 681 À, deux courbes ont été reproduites. La 
première est un interférogramme fait avec une 
épaisseur de 4,122 mm, séparation bien appropriée 
à la mesure des composantes fortes, serrées. La 
seconde est faite avec une très longue pose et une 
épaisseur de 2,461 mm. Les composantes faibles a 
et b sont très clairement visibles, mais les compo- 
santes fortes, de 7 à L sont très fortement sur- 
exposées ; cependant d et f peuvent être mesurées 
de façon satisfaisante, ce qui rend possible de relier 


450 


les nombres d’ondes de a et b à ceux des compo- 
santes fortes de g à l, puisque d et f peuvent aussi 
être bien mesurées sur des poses faites avec l’épais- 
seur de 4,122 mm. Il faut noter que le zéro dans 


m np d e F qe 1 k t mn  d 
A f\ 
FA 
| fut. 4e 8 
1 BR 
Vf A: ur 
\ N 11: 
À HAE UV 
\ J°} Sn, / ÿ | 
i 154 / : 
| LEP FA : | 
| 7 A1 M) À 
4 # V W i A 
‘4 7 | 
|  - 
ü 200 400 600 600 000 1200 
; mk 
m7 (0 = + 624 mK) é 
É ie ROUE £ x 
Fi. 3. — x 4 681 À. Épaisseur de l’étalon : 4,122 mm. 


les courbes montrant les composantes a et. b se 
rapporte à la composante 4, tandis que dans les 
autres courbes le zéro se rapporte à la compo- 
sante d et deviendrait 624 mK si le zéro se rap- 
portait à a. La composante c n’est pas bien résolue 
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dans la courbe reproduite et se présente seulement 
comme un épaulement à la gauche de d ; elle est 
mieux résolue sur les clichés moins surexposés et 
aussi sur ceux qui sont faits avec une épaisseur 
de 3,595 mm. ; 


Structure hyperfine des termes. — Puisque tous 
les termes du multiplet 4 %F,,,,,et deux de 5 5D,,,, 
interviennent dans deux raies, les transitions corres- 
pondant aux composantes hyperfines observées 
peuvent être déterminées sans ambiguïté. Dans le 
tableau IV sont données les transitions possibles 
pour toutes les raies observées et les attributions 
des composantes hyperfines. 

A partir de cette analyse on peut calculer les 
intervalles des termes hyperfins pourvu que chacun 
des termes hyperfins intervienne dans au moins 
ure composante observée sans mélange ; et cela 
est le cas pour 29 des 30 intervalles hyperfins, la 
seule exception étant le niveau Æ# — 1/2 du 
terme 4 ?F,, de sorte que le seul intervalle non 
déterminé est celui qui sépare les niveaux F# = 3/2 
et F — 1/2 de ce terme. De plus, sur les 29 inter- 
valles déterminés, 18 sont déterminés à partir da 
deux raies différentes, et dans bien des cas, à partir 
de deux ou plusieurs composantes hyperfines dans 
une ou deux raies. La comparaison des valeurs 
trouvées à partir des différentes raies et compo- 
santes donne ure estimation de la précision ; la 


m 
cd  f 9 KL n ae badon 


FiG. 4. — x 4 681 A. Épaisseur de l’étalon : 2,461 mm 
(temps de pose prolongé pour photographier les composantes faibles). 


déviation moyenne est moindre que + 3 mK, et la 
plus grande déviation est -- 7 mK. On peut donc 
se fier aux valeurs données pour les intervalles des 
termes à mieux que 10 mK. 

Les intervalles du terme singulet 4 1F, ont été 
trouvés à partir de la structure de la raie d’inter- 
combinaison 4 674 À. les intervalles du terme 5 D, 
étant connus à partir des structures des raies de 
triplet. Les niveaux F — 13/2 et F — 11/2 sont 


observés seulement dans des composantes qui se 
présentent comme des mélanges avec d’autres com- 
posantes et, par conséquent, seuls les trois inter- 
valles les plus bas sont déduits de composantes non 
mélangées. Ces derniers sont probablement exacts à 
environ 10 mK près, alors que l'exactitude des trois 
termes supérieurs est nettement moindre. Cepen- 
dant les composantes mélangées à partir desquels 
ils ont été déduits sont fines, leur largeur apparente 
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TABLEAU III 


F 5 SD, 5 D, 58D, 
15/2 Vi ee 
; DEA (A RSR AREA NL Een n en TRE ETS 
1572 
747. (740) SD UDC) RARE enr ere 
11/2 
632 (626) 337. (335) — 94 (917) 
9/2 
5170 1512) 285 (274) — 954 (— 751) 
7/2 
390 (399) DD AND 1S) MR ne TE 
5/2 
BRL AO ENS ES AS NT SERRE 
3/2 
3 416 19299 — 1 668 
SP. etc 3 536 1 191 — 1 768 
F LS, LS LSF, &LF, 
17/2 
Be (SD. 2 ER ENS ARE EE 
15/2 
624 (644) = 300 EL Ge) ENCRES 420 (432) 
13/2 
552 (558) 035 1107) — 692 (— 768) 380 (385) 
11/2 
500 (472) 169 (—164) 681 (649) 330 (317) 
9/2 
430 (386) 10004196) 568 (531) 250 (259) 
7/2 
339 (300) 160 (104) #70) (#13) 210 (202) 
5/2 
190 (215) es: BEC MÉTIER TRE 140 (146) 
3/2 
CS DR de ion ET Penn io 
1/2 
> 3.348 + æ — 894 21361 1 730 
& 3 430 
ppt C: 3 536 —_ 884 2007 1 768 


n'étant pas plus grande que celle des composantes 
non mélangées ; les valeurs données pour ces 
intervalles sont probablement exactes à un tiers 
près de la largeur des mélanges, soit environ 20 mK. 

Dans le tableau IIT sont donnés, en mK, les 
intervalles hyperfins pour tous les termes qui inter- 
viennent dans les raies dont la structure a été 
mesurée. Les nombres donnés entre parenthèses 
sont ceux que l’on devrait trouver si la règle des 
intervalles était exactement suivie, calculés d’après 
la largeur totale mesurée pour la structure du 
terme. À l’exception de l'intervalle 3/2 — 17/2 des 
termes 4 %F,, le tableau est complet. Cet inter- 
valle est désigné par x ; la valeur de x est proba- 
blement comprise entre 100 et 150 mK. La seule 
transition dans laquelle le niveau intervient est 
dans la raie 4 681 À, et une valeur de x de 110 mK 
ferait coïncider exactement la composante corres- 
pondante avec la composante très forte corres- 
pondant à la transition 17/2 — 15/2, tandis que 
toute valeur comprise entre les limites précédentes 


rendrait cette composante suffisamment proche de 
la composante forte pour être masquée par elle. 

Les nombres donnés en bas de chaque tableau 
sont les écarts totaux des niveaux, c’est-à-dire les 
sommes de tous les intervalles. Les valeurs de 
Paschen et Campbell sont aussi données à titre de 
comparaison. Elles diffèrent nettement dans le cas 
des termes 5 5D,, 5 3D, et 4%F,, mais pour les 
quatre autres termes l’accord est satisfaisant, les 
valeurs données par Paschen et Campbell étant 
probablement exactes à 30 mK près. La valeur 
élevée obtenue par eux pour le terme 5 #D, est une 
conséquence directe du fait qu’ils n’ont pas observé 
la première composante, a, de la raie 4 368 À. La 
valeur élevée, qu’ils donnent pour le terme 5 DS 
est liée à la présence de la composante à 4 412 mK 
dans la raie 4 685 À ; celle-ci n’apparaît pas dans 
nos interférogrammes ; de plus l’accord excellent 
entre les positions observées et calculées des compo- 
santes mélangées dans la raie 4 681 À ne serait pas 
conservé si le niveau inférieur de ce terme était 
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déplacé de 100 mK vers le bas, comme cela serait 
exigé par la valeur plus élevée de l'intervalle total ; 
la raie 4 681 À n’était pas résolue de façon satis- 
faisante par Paschen et Campbell. La valeur donnée 
par eux pour le terme 4 °F, est simplement celle 
qu’ils donnent pour 5 #D, ; les niveaux de ce terme 
ne peuvent être déterminés qu’à partir d’une réso- 
lution satisfaisante de la raie 4 681 À ou des raies 
infra-rouges 9 217 ou 9 213 À dont aucune n’a été 
résolue convenablement par Paschen et Campbell. 


Comparaison des positions et des intensités 
mesurées et calculées des composantes hyperfines. 
— Dans le tableau IV sont donnés les nombres 
d’ondes correspondant aux transitions possibles, 
calculés d’après les intervalles des termes donnés Fi. 


TABLEAU IV (1) 


À 4 685 À 
TRANSITION Cazc. Mes. 
HRSREVOSISD)S (o] Int Ce} Int 
15/2 45/2 0 19,4 0 17 
13/2 15/2 300 3,4 295 31/2 
15/2 13/2 850 3,3 851 71/2 
13/2 13/2 4150 41,1 1450 6 
112840) 15589 01005 119796 
13/2 11/2 1897 5,5 1893 9 
11/2 11/2 2132 5,3 2425 11/2 
9/2 ,41/2 2 301 6,4 2 296 8 
11/2 9/2 2764 6,92 2756 9 
9/2 9/2 2933 4,7 ES D PEN. 
7/2 Nos 3038 6,2 3042 61/2 
9/2: 7/2 3450 6,2 3447 8 
7/20 7/9 3555 0,1 CRE es 
5/2 7/2 3160005; 30095 
7/2 5/2 31040, 39950 4 
5/2 5/2 3995 0,3 ART me 
3/2 5/2 £030 3,1 4040 4 
5/2. 3/2 2975 3,1 
3/2 3/2 4310 26 A 
À 4 681 
TRANSITION GAL C. Cazc. (Bl.) Mes. 
LS 25 D: () Int. 0] Int. Comp. Cj Int 
19/2 004572 0 0,2 Résolue a (2002 
11/2 13/2 298 0,5 » b 297 0,5 
CDR TE) 545 0,9 » c 543 4 
45/2 15/2 624 9,7 » d 624 3 
2/2 9/2 747 1,3 » e 74h A 
13/2 13/2 850 P#A » f 849 51/2 
52: IS 925 4,7 » £ 924, 9 
41/2 11/2 41-045 5,2 » h ‘1033. 4/2 
3/2. 5/2 4 125 2,0 5 
9/2 9 UE SA 1162 7,4 j 1163 4 
7/2 7/2 À 264 4,9 
17/2 -15/2 1 302 20,0 
5/2 5/2 1 315 3.9 1298 28,8 k 1298 38 
4/2 3/2 1405—7x 9,2 
3/2  - 3/2 À 405 2,5 EVA RES A 
19/2 04972 1 474 15,1 Résolue l 1747920412 
2 3 
RE ue . ds l 1597 12,0 m AA 
7/2 5/2 4 65% 9,7 
41/2 9/2 4 677 7 | 1107720415,1 n 1667 13 
9/2% 7/2 1 694 4,8 | 
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TABLEAU IV (2) 


À 4 657 À 
TRANSITION CALC. Mes. 
LOF CS ED; Co] Int. CJ Int 
43/2 13/2 0 21,5 DATE 
19/204741/2 380 6,5 378 13 
4152 13/2 692 6,5 681 7 
41/2 11/2 4072 8,1 1059 3 
11/2 9/2 1409 9,4 1396 16 
9/2 11/2 1703 9,4 1696 7 
DIDeE 11979 2040 1,2 PR LE 
9/2 7/2 25329 93 29140 1% 
pale) 9/2 2608 9,3 2597 5 
a ee 99 2893 0,4 2 ÉTRN 2 
ma )2 3113 6,3 3092 7 
5/2 7/2 38631 6,3 3352 4 
5/2 5/2 3 083 5,7 3 974 8 
À 4 656 À 
TRANSITION CALC: MEs 
LISE AS ÉDE Co} Int. C) Int 
45/2 43/2. 0 22,8 0 21 
13/2 13/2 300 4,6 305 5 
11/2 13/2 535 0,5 Li Lee 
13/2 11/2 680 15,4 683 14 
41/2 11/2 95. 7,2 922 10 
9/2 41/2 108% 1,5 1085222 
11/2 9/2 1252 9,5 1954 9 
9/2 9/2 1421 7,9 1422 9 
7/2... 9/2 1526 2,7 1523 &4 
. Le 1 706 4,9 4 710 5 
122 5/2 1 811 7,0 
5/2 7/2 1861 41 DORE 
2 2e 2 031 1,8 RC ASE : 
5/2 5/2 2081 4,5 
3/2 5/2 2116 5,7 L LEE 


- dans le tableau III. Les intensités calculées sont 


24 #0 


aussi données. A titre de comparaison les nombres 
d’ondeset les intensités mesurées (réduites en pour- 
centage comme pour les intensités calculées) sont 
aussi donnés. Dans le cas des raies 4 681, 4 652 et 
4 674 À, dans lesquelles un grand nombre de tran- 
sitions donne des composantes si serrées qu’elles ne 
sont pas résolues sur les interférogrammes, on 
donne les positions des centres de gravité des 
mélanges qui en résultent, afin de permettre la 
comparaison avec les structures mesurées ; on a 
supposé que les raies séparées de moins de 50 mK 
sont mélangées. Les positions des centres de gravité 
sont calculées d’après les intensités calculées ; cette 
façon de faire est correcte pour la raie 4 681 À, 
pour laquelle l'accord entre intensités calculées et 
mesurées est bon ; mais pour les raies 4 652 et 
4674 À, pour lesquelles il y a des différences 
notables entre les intensités calculées et mesurées, 
il reste une possibilité d’erreurs de l’ordre de 10 
à 20 mK dans le cas des mélanges les plus larges. 

L'accord entre les positions calculées et mesurées 
est excellent pour les composantes résolues et pour 
les mélanges. Dans la plupart: des cas les diffé- 


rences sont moindres que 10 mK, et dans le cas des 
raies 4 681, 4 656 et 4 638 À, qui peuvent être 
mesurées plus exactement en raison de la plus faible 
largeur totale de la structure qui permet d'utiliser 
un étalon plus épais, la différence moyenne est 
moindre que 3 mK. 


a De d e F g °h x a bcde fig 0h >x 
AA 
MA : 
À 0 À LA 
} ‘ f # À la ‘ À V / 
PA HR. ; À (#2) [2 
Po k f V \ w Ÿ É À 
À : | L à f Lt Î #4 
y UN 
Ni 
| | | | | | AUS | | 
c) 400 600 1200 1600 0 400 800 200 1600 | 
2022 


ARE] 
| 


Fic. 6. — à 4 652 À. Épaisseur de l’étalon : 2,460 mm. 
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TABLEAU IV (3) 


JOURNAL DE PHYSIQUE 


À 4 652 À 
TRANSITION CALc. Cazc. (Bl.) Mes 
ANR SIDE (} Int. (y Int. Comp. c Int. 
5/2 7/2 DTA 
3/2 5/2 304857 DHEA a DOTE 
To Ole EURE 24) 
9/2004472 ARTMIASS 
5/9: 45/2 220 4,5 220 1065 ù 
11/2 13/2 267 0,5 : 
72 7/2 339:-27,0 Résolue b 330 : 12 
9/2 9/2 48h 7,9 » 490 12 
7/2 5/2 559. 4,8 » ANR 
14/2072 647 FI » d 649 12 
9/2 7/2 769 4,9 
HE en 793 95 e 795 121 
41/2 9/2 984 9,5 Résolue Î 079 
JON AMOOMETOSE » £ 1201720 
45/2048/2 1443 22,8 » k 4445 20 
À 4 638 À 
TRANSITION CALC. MES 
LACS RAIN D Co] Int. Co) Int 
92 2e 0 4,9 0 6 
11/2 9/2 123 14,2 119 144 
43/2 11/2 345 28,8 339 28 
CAEN AE 568 9,8 562 10 
9/2 9/2 754 12,9 750 12 
11/2, 41/2 1037 9,8 
5/2 9/2 1038 12,0 1085 20 
719: 00972 1322 : 6,2 1320 8 
9/2 41/2 11608 229 1667 2 
À 4 674) À 
TRANSITION CALC. Cazc. (Bl.) Mes. 
LRO EDS Cy Int o Int. Comp. Cy Int 
13/22245)/2 Omar Résolue a Pt 2 
15/2 45/2 420 19,4 
192 139 70 58 431 24,9 b 497 +022 
13/2 13/2 850 11,1 
9/2 11/2 887 64 864 17,5 € 8614 13 
11/2 11/2 1217 5,3 
NET 6. 1269 6,4 4251 15,1 412149530045 
HN Ph EU 17000 SE / 
9/2 9/2 2 MS Ce EE 
512 20 978 ‘4576 5,1 | 1584 11,8 e 1576 9 
13/2 11/2 1597 5,5 
3/2 5/2 41 826%%"3;1 Résolue Î 1830-0138 
5/2 5/2 1966 0,3 » ER 
9/2 7/2 2036 6,2} ñ 2 
3/2 3/2 2106 26 2 057 8,8 £ 2045 43 
7/2 5/2 D TE | Résolue h 2 185 9 
5ÿ2 319 29246 3,1 » j D'OATUASE 


Dans le cas des raies 4 681, 4 652 et 4 674 À, 
une proportion relativement grande des compo- 
santes sont groupées en mélanges assez serrés ; cela 
donne une explication satisfaisante de la diver- 
gence entre le nombre de composantes calculées et 
le nombre plus petit de composantes observées, de 
largeur apparente non supérieure à la largeur de 
composantes simples ; et les positions calculées des 
centres de gravité de ces mélanges sont en excellent 


accord avec les mesures, surtout si on considère la 
possibilité d'erreurs dues à la différence entre les 
intensités calculées et mesurées. 

La comparaison entre les intensités calculées et 
mesurées n’est pas sans intérêt. Pour la raie 4 638 À 
lPaccord est excellent, bien au delà de la précision 
expérimentale ; pour la raie 4 681 À aussi, l’accord 
est dans les limites d'erreurs, sauf pour la compo- 
sante la plus forte et sa plus proche voisine ; mais 


N°3 


cette composante est tellement plus intense que les 
autres que l’estimation de son intensité est plus 
difficile, ainsi que celle de sa voisine. Il apparaît 
donc que pour cette raie aussi il y a un bon accord 


/ N4 À } a 
Y MW | 


Fic. 7.— 214 638 À. Épaisseur de l’étalon’: 2,460 mm. 


entre les intensités observées et calculées. Pour la 
raie 4 656 À, l'accord entre les intensités calculées 
et observées est bon aussi, la plus grande diver- 


_ gence n’excédant pas 30°%, sauf pour le mélange 


à 1 824 mK, qu’on pourrait s'attendre à voir faible, 


parce que la composante plus faible était presque 


résolue. 

Mais pour les autres raies il y a des différences 
qui sont bien au delà des limites d'erreurs expéri- 
mentales. Ainsi pour la raie 4 685 À, qui est 
presque exempte de mélanges, l’intensité calculée 
de la composante à 850 mK est moindre que la 
moitié de la valeur observée, alors que celle de sa 
voisine, à 4 150 mK est presque deux fois plus 
grande ; ces intensités observées sont en accord 
avec celles de Paschen et Campbell. De même, 
pour la raie 4 657 À, la différence entre les inten- 
sités calculées et observées des composantes à 0 et 
380 mK est, de loin, plus grande que l'erreur expéri- 
mentale ; et ici aussi l'accord est excellent avec les 
mesures de Paschen et Campbell. Pour la raie 
interdite 4 652 À, l'intensité des six premières 
composantes est beaucoup plus faible que les 
valeurs calculées ; les trois premières forment un 
mélange ayant environ un tiers de l'intensité 
voulue, alors que le second groupe de trois, qui 
devrait aussi former un mélange, n'apparaît pas 
du tout ; ceci est dû, en partie, à la grande largeur 
sur laquelle ce mélange est étalé, mais si les inten- 
sités étaient conformes au calcul, la composante la 
plus forte devrait être à elle seule, environ deux 
fois moins forte que les composantes à 339, 484 et 
649 mK. La raie d’intercombinaison, 4 674 À, 
montre aussi des différences notables entre inten- 
sités calculées et mesurées, la plus importante étant 
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pour la composante à 1 836 mK dont l'intensité 
calculée est 3,1, contre 13 mesurée (la valeur de 
Paschen et Campbell est probablement plus faible 
puisqu'ils n’ont pas observé cette composante, mais 
leur résolution de cette structure est très incom- 
plète). 


1600 2000 mK À 


Fic. 8. — à 4 674 À. Épaisseur de l’étalon : 2,460 mm. 


Il faut remarquer que les valeurs expérimentales 
des intensités sont les moyennes prises sur toutes 
les plaques mesurées. Il y a des différences notables 
même entre les différentes plaques faites avec la 
même épaisseur d’étalon (auquel cas il n’y a pas 
d'influence de la différence de dispersion). Ces dif- 
férences ne se présentent que dans les raies pour 
lesquelles il y a désaccord entre les intensités 
calculées et mesurées:elles sont de l’ordre de --25 % 
et peuvent présenter des inversions entre intensités 
de composantes voisines ; par exemple sur une 
pose deux composantes peuvent avoir les inten- 
sités 5 et 6, alors qu’elles sont 6 et 5 sur une autre 
pose. Mais ces variations dans les intensités sont 
beaucoup plus petites que la plus grandes des dif- 
férences entre intensités calculées et mesurées. 


Diseussion relative aux largeurs totales et aux 
rapports d’intervalles dans la structure hyperfine 
des termes. — Les configurations électroniques des 
termes 5 D et 4 °F sont respectivement 5s5d $D et 
5s4f 3F ; en conséquence la structure hyperfine de 
ces termes est presque entièrement due à l’élec- 


tron 5s pour lequel la valeur moyenne de (1/rà) 
est relativement grande alors que la contribution 
des électrons 5d ou 4f est beaucoup plus faible. 
L'ordre de grandeur de l’écart total dû à des 
derniers est donné par 2M.u,/u,.2/Z, où M est 
l’écart du multiplet, u, et u, les moments magné- 
tiques respectifs du noyau et de l’électron et Z le 
numéro atomique ; les valeurs de SM sont respec- 
tivement de 220 et 22 mK pour les termes 5 #D et 
43F, ce qui donne des écarts hyperfins totaux de 
l’ordre de 0,02 mK et 0,002 mK. Ces quantités sont 
négligeables et si l’on considère seulement l’inter. 
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action du moment magnétique nucléaire avec l’élec- 
tron 5s, les écarts totaux des termes devraient être 
(en admettant la valeur 3418 mK pour le 
terme 5 #)),) de — 1 708 mK pour le terme 5 5D,, 
3 418 mK pour 4%F, et 2275 mK pour 4 °F, 
Dans le cas des termes 5 51), et 4 3D,, les valeurs 
mesurées diffèrent de ces valeurs de 40 mK et 
85 mK. Ces divergences ne sont pas grandes, mais 
elles sont bien supérieures à l’erreur expérimentale 
probable ; elles sont dues en partie au fait que la 
variation d'énergie de l’électron 5s en fonction de 
l’orientation nucléaire affecte aussi l’énergie de 
l’électron 5d et, à un moindre degré, celle de l’élec- 
tron 4f. 

Les intervalles, dans les termes 5 3D,, ,,, sont, 
ainsi que le montre le tableau ITT, en bon accord 
avec ceux qu’exige la règle des intervalles (pour les 
valeurs observés des écarts totaux). Pour les 
termes 5 #D,, et 5 5D)., les différences entre les 
valeurs calculées et mesurées sont bien inférieures 
à l'incertitude expérimentale, la plus grande dif- 
férence étant 9 mK et la différence moyenne infé- 
rieure à 6 mK. L’accord n’est pas aussi bon pour le 
terme 5 #D,, où la plus grande différence vaut 
16 mK et la moyenne 9 mK ; mais ces différences 
ne sont pas suffisamment supérieures à l'incertitude 
expérimentale pour qu'on puisse les considérer 
comme significatives. On doit s’attendre, naturel- 
lement, à ce que ces écarts soient conformes à la 
règle des intervalles malgré le moment quadru- 
polaire exceptionnellement grand de In puisque 
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la structure hyperfine, étant due presque entiè- $ 


rement à l’électron 5s, ne doit pas être influencée 
par ce moment quadrupolaire. 

Pour les termes 4 3F,,,, les différences entre les 
intervalles observés et ceux conformes à la règle des 
intervalles sont largement supérieures à l'erreur 
expérimentale probable; et, de plus, il y a une 
régularité dans ces différences. 

Ces divergences peuvent être attribuées à une 


perturbation mutuelle entre les termes, les écarts, 


du multiplet (45F, >=— 43F, — 10/7 K4et 


ASF, — 43F, — 5,2 K) n'étant pas beaucoup” 


plus grands que les écarts de structure hyperfine. 
C’est ainsi que, pour les grandes valeurs de F les 
intervalles sont trop grands dans 4 °F, et trop 
petits dans 4 %F,, alors que pour les petites valeurs 
de F c’est l’inverse qui a lieu. 
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Pour le terme 4 1F, les intervalles suivent la M 


règle des intervalles dans les limites de l’exactitude 
expérimentale ; on doit en conclure qu’il n’est pas 
perturbé par les termes de triplet, bien que le plus 
proche d’entre eux ne soit distant que de 34,3 K. 


L'auteur est heureux de profiter de cette publi- 
cation pour remercier le Pr Jacquinot d’avoir mis 
à sa disposition les moyens du laboratoire et de 


l’intérêt qu'il a constamment porté à ces recherches, 


ainsi que M. Blaise des discussions profitables qu'il 
a eues avec lui. 


Manuscrit reçu le 30 novembre 1956. 
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RÉSOLUTION DE L'ÉQUATION DE SCHRODINGER DES ATOMES A DEUX ÉLECTRONS 
II. MÉTHODE RIGOUREUSE. ÉTATS S SYMÉTRIQUES. 


Par G. MUNSCHY et P. PLUVINAGE, 


Institut de Physique de la Faculté des Sciences de Strasbourg. 


Le présent travail fait suite à un article paru 
il y a environ deux ans [1], et qui était intitulé : 
« Approximations systématiques dans la résolution 
de l'équation de Schrüdinger des atomes à deux 


électrons. — I. Principe de la méthode. États S 


symétriques. » La référence" correspondante sera 
désignée par PI. Le long délai qui s'est écoulé 
depuis cette première publication a été consacré à 
la recherche d’une justification mathématique du 


. procédé d’approximation utilisé. Il est apparu que, 


moyennant quelques précautions dans la mise en 
œuvre de l’idée directrice, on pouvait résoudre 
rigoureusement le problème posé. Le principe de la 
solution a déjà été indiqué brièvement dans une 
note aux Comptes Rendus [2] mais, depuis, un 


. nouveau progrès, important par ses conséquences 


pratiques, a été réalisé et nous donnons ici une 
mise au point définitive. 


1. Transformations de l’équation de Schro- 
dinger. — Les notations sont, en général, les mêmes 
qu'en PI. Nous rappelons les plus importantes. 
L’atome se compose d'un noyau de charge Ze et de 
deux électrons. En fonction des variables classiques 
de Hylleraas r;, ro, ri, OÙ S$, t, re les variables 
que nous introduisons sont 


S = Fr +E ro D At nt (1) 
S S 
Le domaine de b et de + est un triangle & 
DDR —p<T< 0. (2) 


Nous nous bornons encore à l'étude des états S. 
L'élément de volume de l’espace de configuration, 
réduit au domaine des variables s, o et +, est 


du = s5 pfi — +?) ds de dr. (3) 


Un facteur 8r?, sans importance, est omis 
dans dv. 

On peut donner du problème à résoudre, soit 
l'énoncé intégral 


sf(ve + Fy— Eÿ) do = 0 


soit l'énoncé différentiel de Schrüdinger en fonction 
de l'énergie cinétique T'et de l'énergie potentielle F 


Hÿ= (T + V) ÿ = Ey. (5) 


(4) 


En PI sont exposés les avantages du changement 
de variable et d'inconnue 


DES ei —— FH = 0) + >0 (6) 
® — et (x, p, 7). (7) 

L'équation (5) devient 
HQ = EQ + &ee-ai [at p(1 — )] Co. (8) 


© est un opérateur dans lequel nous distinguons 
deux termes 


O — es — GB, (9) 
( a? f 2 T” ra ° 12 / 
À = — pli — r?) + — ex pf1 — +?) Tee? (10) 
& & = 
OR: fer 
GG —=={(kZo—1 + +) (11) 


Il reste à exprimer l'opérateur énergie ciné- 
tique 7° en fonction des variables +, 9 et +, et à 
effectuer le produit T'exp (— x/2) Le résultat peut 
se mettre sous forme condensée en posant, comme 
en PI, 


L = (1 + +?) D[o(t — p°) D, + {1 — p°) d:] 
+ Det + r) (1 — 6°) D, + of — +1) d.] (12) 
M = T1 — p°) d, + pTÛL — 7°) de (13) 
L'opérateur & s'écrit alors 
A = — x? pf1 — +?) (Oise — de) 
—$ + AIN — xf4 — 2 D,) IN — Aro de + 2x0 (14) 


Il ne dépend plus de & ; l'opérateur (9) est donc 
fonction linéaire de &. 

La fonction (7) est solution de l’équation de 
Schrüdinger (5) si la fonction 9 vérifie 


(15) 


Op = (< — GB) o = 0. 


Pour que la solution soit physiquement accep- 
table, il faut que © soit partout régulière sauf 
pour æ — ©, où, cependant, elle ne doit pas aug- 
menter plus vite que exp (x/2). 


2. Choix des fonctions de base. — Les singula- 
rités importantes à distance finie del équation (15), 
obtenues en annulant les coefficients des dérivées 
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secondes, sont toutes situées sur la frontière du 
domaine de définition de x, p et +. Ce sont 


æ— 0, pet + quelconque, (16) 
p — 0; T0, + quelconque, (17) 
pu, TL, iriquelconque: (18) 


Nous allons définir un système de fonctions, 
partout régulières sauf pour x — co, qui satisfont 
à des équations différentielles ou aux dérivées 
partielles présentant ces mêmes singularités, et qui, 
par suite, permettront un développement relati- 
vement simple de +. Ce sont des produits de fonc- 
tions de + par des fonctions F de b et +. Les calculs 
effectués en PI ont montré que les premières sont, 
pour les états symétriques, soit des polynômes de 
Laguerre Li(x) soit des fonctions æLi(x). Mais la 
détermination des F n’avait pas été effectuée d’une 
manière satisfaisante et c’est à ce défaut que nous 
allons remédier(?) . 

Le point de départ est une équation aux valeurs 
propres étudiée par P.-Appel et J. Kampé de 
Ferrié [3] qui s’écrit originalement 


dE[E(L — Ë) de F — End, F] 
+ Dali — 1) 0nF— End, F]=RKF 


EE 0 An OMC ENREAE 
Le changement de variables 

so Ni 

et Chat SU TIRE MT. 


permet de se ramener au domaine de définition 
de pet et l’équation devient 


EF = D,[efi — 0) d, F + 71 — bp) d. F] 
moelle) DOME (p — +?) d- F) HCHs (19) 


Ses singularités coïncident avec (17) et (18). 
Appel et Kampé de Ferrié ont démontré que le 
spectre de X est donné par 


k = — mm + 2) IN OMINNASENSE (20) 


Les fonctions propres sont des polynômes en p 
et +, qu'on détermine par une seconde condition. 
Au lieu de s’en tenir à la condition posée par Appel 


et Kampé de Ferrié, il est plus avantageux ici 


d'imposer aux polynômes d’être fonctions propres : 


d’un second opérateur G commutant avec Æ, et de 
choisir pour G l’un des opérateurs auto-adjoints 
envisagés par G. Munschy [4]. C’est l'opérateur 
désigné primitivement par &, et qui s'écrit, avec 
les variables p et 7, 


G = d:(p? — 7°) d+. (21) 

Alors les fonctions F seront définies par cn 
KF =— mim +2)F, (22) 
GE = gr. (23) 


() La notation F(9, +) remplace la notation Y(e,+) de PI. 
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La seconde équation se ramène à celle des poly- 


nômes de Legendre P, quand on prend +/9 comme 


variable. Posant alors 


F=pnP, () J(e) 


et substituant dans (22), on trouve pour J 


p(1 — 6).dés J + [2n + 2 — (2n +3) pd, J 
+ (m— 2) (m En 2) J = 0,N(28) 


C’est l'équation d’une famille de polynômes de 


Jacobi qui, avec les notations de Morse et 
Feshbach [5], s’écrivent 


Jm-nl2n + 2, 2n + 2]e) 


2 1)! din 
— ® _— ; OT 2H dons [omtm+i(t — p}"”-"] 
(25) 
n < m. IS (26) 


Aïnsi les fonctions F sont définies par 


T 


FANS Dh Pi () Jn-nl2% + 2,.2n 20 CAR 


Kg = mimi (28) 


GEmn = — nn + 1) Fyn- (29) 


Elles ont en 7 la parité de n. Les fonctions paires 
seront notées F#%,. Elles permettent de former les 
fonctions d’onde des états symétriques. Les fonc- 
tions impaires, qui seront désignées par F%, servent 
pour les états antisymétriques. Les fonctions Æ 
forment un système orthogonal avec la densité 1. 


1 +e 
L de / Fu ER dr 
0 —? 


\ ï 


LR ({m— n) ! (2n +1 
(m + 1) (2r + il 


(RENTE CE 


] * Sim Sim (80) 


Les premières polynômes ont les expressions 


suivantes 
F9 = 1, Pas 

F9 = 15e, FR = rs 

Fo = —Se, FR = RSS 
phni pt ste EE 


Les fonctions de base utilisées pour les états S. 


symétriques seront 
(Noop) 2 Lil) Foo 


(Nnnp) */? æLi(x) Fbnle, +), 
î RAS 


pi="0,40228 
Umnp(T-P; 7) = 


; n pair < m. 
(31) 
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* Pour faciliter les calculs qui vont suivre, les 


facteurs de normation sont définis par 


PET gt e-2(13) de [7 & [ * Fe F5) dr 
ue WP: api 
PE (2) 
“A = fi x e#(xLi) j° de fa ef se, Nadar 
(p + 5)! [(m—n) !(2n +1)1/? . (33) 
pl{m +1) (2n +1) [(m+n + 1) 1? 


.! 3. Résolution de l’équation de Schrôdinger. — 


les équations (9) à 


Cherchons s’il existe pour la fonction inconnue o, 
un développement de la forme 
m_ co 


QUE 5 Z ZE œmnp Umnp(Z; Ps T) 


(34) 
m=0n=0 2=0 


(n pair ; »m’ plus grand entier pair < m). 


Les coefficients «,,, doivent être déterminés de 
façon que l’action sur + de l’opérateur © défini par 
(14), soit nulle. Si les a, étaient 
toutes orthogonales avec la même fonction de poids, 
il suflirait de calculer les éléments de matrice de © 
et de développer © suivant les v,,,. Ici une petite 


complication surgit parce que la fonction de poids 
est 2er pour. les uw, et æte-* pour les u,,, 
quand » n’est pas nul. Nous allons montrer que 


cette difficulté peut être aisément surmontée. 


Les propriétés des #,,, entraînent des simpli- 
fications notables dans le résultat de leur trans- 
formation par ©. Avec des notations claires par 
elles-mêmes, 


OLD oo — Corel Fo + DES Ofup &LS Far, (35) 
Ï aTs 
PR Or 2 Le Fo + DECO, xls. FO. (36) 
8 Grs 


Ces équations résultent de la façon dont les u,» 
dépendent de x et sont déjà implicitement con- 
tenues dans PI. Ce qui est nouveau ici, c’est d’abord 
l’orthogonalité, avec la fonction de poids a* e7, 
des fonctions xL£ F0, xL5 FS aux fonctions Lé Fi 


.2L, F5. En adoptant la notation 


(37) 


co 1 +p 
ly, Ou) = feat de [def you ds, 
0 ai 


et en utilisant les normations (32) et (33), nous 


obtenons 
Ofon = (Noop) (Ls Fo: OL; F0); 
Oop — (Nars)—* (xLe Far» OLS Fo), (38) 
One — (None) —* (LS Foo OZL$ Far), 
OU — (Nour)? (2LS FE), OxLs Fh). 
Le progrès essentiel sur le travail exposé en PI 


est que la plupart des intégrales sont nulles à cause 


des relations de récurrence, entre les F,,. Ces rela- 


' 


sn “HR 


tions seront établies dans la thèse de doctorat de 
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G. Munschy. On y démontrera que le nombre des 
intégrales non nulles pour une valeur fixée, soit des 
indice inférieurs, soit des indices supérieurs est 
soixante-trois au plus. 

Si nous opérons sur les fonctions normées 
nous obtenons 


oi's Nynnp FT 9 
mn? = niet (Ugrss Ounnp)s 
grs 


MAD 


Les nouveaux coefficients 


j np = (uors, Oumnp) (39) 


sont symétriques. Pour établir cette propriété, il 
faut démontrer que © est auto-adjoint avec la 
fonction de poids 4e. Or, d’après (3) et (6), 
l’élément de volume de l’espace de configuration 
est proportionnel à x° p(1 — +?). 

Pour deux fonctions 4 et « telles que 


d => Ce ®; 
il résulte de (8) et de l’hermiticité de Æ, que 


oO —= e—*/2 X> 


fu ERA 00 


fe x? (xOp — PO x) dx de dr = 0, 


22e) 


ce qui prouve la propriété énoncée. 

Substituons maintenant le développement (34) 
dans l’équätion (15), et annulons les coefficients 
de Züoop €t de Umnp Nous obtenons le système 
linéaire 

Fo = Oour Goo + ZE Onap Xqrs = 0, 
Fe (40) 
Fmnp = >? One %oos + 222 Omnp dns = 0. 
s ars 


Le nombre des inconnues par équation est 
soixante-trois au plus. D’après (9), chaque coef- 
ficient est fonction linéaire de e et l’on peut poser 

Omno = A Anny CHERE Brines 

Les matrices À et B associées aux opérateurs & 
et G sont symétriques. La forme (10) montre que 
la matrice À est celle d’une forme quadratique 
définie positive mais on ne peut rien dire de B. 

Voici les premières équations (40) 


184 il 8V30 8V30 
+ (Re?) caen (OR e+ ee 


5 5 105 105 
199V6 263V6 
ps ; Le, = 0 
( 105.77 405 ) aan ie 


460 
D UNS 
T'920 —= 9 Xooo 
8V30 8V30 1 328 75 
105 © Ÿ qus ) 200 F 63 ue 7220 
8vs 3v5 
24) 200 += 0. 


L’équation en #,,59 comprend douze termes, celle 
en F»p, vingt-quatre. Nous pouvons nous donner 
arbitrairement un certain nombre d’inconnues 49 
et calculer toutes les autres de proche en proche. Par 
exemple, dans la première équation, trois inconnues 
sont arbitraires et la quatrième s'exprime en fonc- 
tion des premières. Nous déterminons bien, ainsi, 
pour une valeur quelconque de e, une solution de 
l’équation de Schrôdinger mais cette solution pré- 
sentera en général pour æ — co une singularité qui 
empêchera la normation. Il y a exception si le 


déterminant est nul 
D = |Ae— B| = 0: (41) 


Toute racine positive e; de cette équation corres- 
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pond à une solution physiquement acceptable. 
Elle détermine une fonction +, une fonction 
d’onde {; et un niveau d'énergie £, 


Tr Tite ll 
: 1274 2 
Vi —= TES p; |s, —;: 

g UE S 


JE = — à. 


Dans la pratique, on opèrera avec un nombre 
fini N de termes dans le développement (34) et on 
annulera les N premières formes (40). Les racines e; 
du déterminant ®,, du système seront des valeurs 


approchées des s;. Les fonctions propres approchées 


correspondantes vérifieront, d’après leur cons- 
truction même, 


(Ÿs, Hÿ'i) = —ef(p, d'). 


La valeur moyenne de l’énergie pour la fonc- 


(42) 


_tion ; est donc — £;? et, comme dans la méthode 


de variation, — e; est une approximation par excès 
de la valeur asymptotique — &;. 

Les questions d’approximation numérique seront 
traitées dans un prochain article. 


Manuscrit reçu le 12 janvier 1957. 
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NOUVELLES MESURES DES RAPPORTS GYROMAGNÉTIQUES DU FER ET DU NICKEL 


Par ANDRÉ J. P. MEYER et Sxezron BROWN (1), 
Laboratoire Pierre-Weiss, Institut de Physique, Strasbourg. 


I. Introduction. — Le rapport gyromagnétique 
d’un corps magnétique est défini par le rapport de 
son moment cinétique total à son moment magné- 
tique total. Dans une expérience gyromagnétique 
on observe la variation de l’une de ces grandeurs 
entraînées par la variation de l’autre : o = AJ/AM. 
En ramenant le phénomène à l’échelle atomique : 


e — A(Japin de Jorbite) 
A(Mspin cie Morbite) 


et en introduisant 


Jon Le Jan ne 
Mpin e Morb. e- 
“onobtient 
55 AMorb 
—| — 1 —— 
2 5) ; gi A(Mapin + Morb.) 


Suivant la contribution relative du moment orbi- 
tal au moment total le rapport gyromagnétique 


que bo NC : 
exprimé en unités (F sera donc un nombre égal 
ou supérieur à 1 si Moy et Mo Sont de même sens. 
En 1953, les résultats expérimentaux connus 
pour les ferromagnétiques étaient peu concordants : 
suivant les auteurs, les rapports gyromagnétiques 


se groupaient près de p _ = 1,00 (Chattock et 


Bates [1], Sucksmith et Bates [2], Coeterier [3], 


Meyer [4]) soit en valeurs légèrement supérieures 
à 1,00 (Barnett [5], Ray Chaudhuri [6], Kikoin et 
Gubar [7], Barnett et Kenny [8], Scott [91, [10]) 


. sans concordance marquée entre elles. 


La précision des expériences récentes permettait 
de supposer qu’une erreur systématique devait 
affecter l’un de ces groupes de résultats. C’est afin 
de mettre en évidence cette erreur que fut décidée 
la construction d’un nouvel appareil de mesure de 
l’effet Einstein-de Haas bénéficiant de l’expérience 
du laboratoire de S. J. Barnett, et de celui de 
Strasbourg. C’est au cours de la réalisation de ce 
nouvel appareil de mesure que put être mise en 
évidence l'erreur systématique qui affectait le 


groupe des résultats à p (£) & 1,00f10]. Il s’agis- 


sait’ d’une erreur commise dans la détermination 


(2) A présent à Fresno State College, Californie. 


de la valeur angulaire de l’amplitude de résonance 
dans l’effet Einstein-de Haas. Cette amplitude est 
généralement déterminée par lecture sur une échelle 
de la déviation d’un faisceau lumineux réfléchi par 
un miroir porté par le système oscillant, l’échelle se 
trouvant à une distance connue du miroir. Cette 
détermination est exacte à condition qu'il n’y ait 
pas de lentille entre le miroir et l’échelle. Cependant 
dans la plupart des dispositifs expérimentaux on 
utilise un miroir plan devant lequel on dispose une 
lentille qui est traversée deux fois par le faisceau 
lumineux. Dans ce cas la distance optique à utiliser 
dans la détermination de la déviation angulaire est 
célle de la lentille à l’échelle. Il semble que toutes 
les mesures de l’effet Einstein-de Haas à l'exception 
de celles utilisant la méthode de zéro, soient 
affectées de cette erreur qui se traduit par un 
abaissement de quelques pour cent des rapports 
gyromagnétiques. 

Il fut heureusement possible de corriger les prin- 
cipales mesures antérieures affectées de cette erreur 
(Meyer, Scott) et de constater enfin pour plusieurs 
substances, dont le fer, un bon accord entre les 
différents résultats. Les désaccords subsistants 
pouvaient être attribués à des erreurs non systé- 
matiques d’une expérience très délicate ou à des 
défauts des substances comme dans le cas du nickel 
par exemple (voir chapitre V). 

Après la découverte de cette erreur systématique 
la mise au point de l’appareil fut achevée, en vue 
de procéder à de nouvelles mesures sur les prin- 
cipaux ferromagnétiques pour lesquels la concor- 
dance des résultats antérieurs laissait encore à 
désirer. 


II. L’appareil de mesure. — Le nouvel appareil 
de mesure de l’effet KEinstein-de Haas utilise comme 
l’ancien la méthode de résonance forcée dont nous 
rappelons brièvement le principe. L’aimantation de 
la substance est inversée synchroniquement, lors de 
chaque passage à la position d'équilibre, par un 
dispositif électronique. Un spot lumineux étroit est 
réfléchi par le miroir porté par le système oscillant 
et passe périodiquement sur une fente disposée 
devant une cellule photoélectrique ; son empla- 
cement coïncide avec la position d’équilibre du 
spot. La cellule photoélectrique actionne l’inver- 
seur du courant de la bobine d’aimantation. Celle-ci 
a été choisie mobile afin d'éviter les interactions 
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indésirables entre les composantes horizontales du 
moment de la substance et du champ de la bobine. 

Dans ce cas, le rapport gyromagnétique se déduit 
des données expérimentales par la formule suivante 


tenant compte de l'effet gyromagnétique propre. 
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des électrons de conduction du courant dans la 
bobine d’aimantation : 


Fic. 2. — Vue d'ensemble du nouvel appareil de mesure de l’cffet Einstein-de Haas. 


TO 
M : moment magnétique de la substance étudiée ; 


M' : moment magnétique de la bobine ; 
2mc 


: rapport gyromagnétique des électrons de 


conduction = 1,1370.10-7 ; 
2A : double amplitude de résonance ; 
l : longueur du rayon optique ; 
: coefficient de frottement du système 
HCLSIT : 


avec : 


I : moment d'inertie ; 
à : décrément logarithmique. 
T : période du système pendant la mesure de à. 


A l'exception de la bobine d’aimantation mobile, 
le principe de l’appareil est identique à celui déjà 


utilisé au laboratoire P. Weiss et à la description 
détaillée duquel on pourra se référer [4] ; dans la 
suite nous limiterons la description du nouvel 
appareil à ses points essentiels. On trouvera dans la 
littérature de nombreuses discussions de la méthode 
de mesure et des multiples erreurs expérimentales 
à éliminer (Coeterier [3], Meyer [4], Barnett et 
Kenny[8]). 


19 LE SYSTÈME OSCILLANT. — Celui-ci est repré- 
senté en coupe sur la figure 1. La substance (a) à 
étudier se présente sous forme de tige de 153 mm 
de long et de 4,7 mm de diamètre maximum, aux 
extrémités taillées en cônes de 600, dont le cof- 
ficient de champ démagnétisant est d'environ 0,04. 
Ces tiges sont introduites par le bas dans un tube 
de quartz (b) portant la bobine d’aimantation (c) : 
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99 spires/cm sur 14 cm de long d’un fil de cuivre 
émaullé de 0,19 mm de diamètre, en deux couches. 
La tige est bloquée entre deux minces disques élas- 
tiques de cuivre-béryllium (d) solidaires chacun 
d’un écrou. : 

Le haut du système oscillant est conçu de façon à 
permettre l’alignement de l’axe de la tige à étudier 
avec la verticale, c’est-à-dire réduire au minimum 
la composante horizontale du moment magné- 
tique ; à cet effet la tête de fixation (e) du fil de 
suspension peut être déplacée au moyen de trois vis 
à faible pas (f) dans un plan perpendiculaire à 
l’axe de rotation. Le fil de suspension (g) est un 
fil de tungstène de 0,02 mm de diamètre et de 
15 cm de long recuit, sous une tension équivalente 
à celle donnée par le poids du système, dans une 
atmosphère d'hydrogène près de 1 000 ©C : on évite 
ainsi les déplacements lents de la position d’équi- 
libre par suite de fluage du fil. Ce fil est fixé à la 
tête de torsion (h), qui peut tourner autour de son 
axe. Un tube de quartz (i), mince et transparent, 
porte deux miroirs plans (j) permettant les mesures 
dans deux azimuths différents de 1800. Deux fils 
torsadés faisant suite à un petit inverseur AB 
montésur la flasque supérieure delabobine amènent 
le courant à deux plaquettes (k). Deux spiraux 
d’or pur (l) convenablement formés et recuits 
assurent le passage du courant vers les bornes 
fixes (m) solidaires de la tête de suspension. 

Le bas (n) du système oscillant peut être pris 
dans un système de blocage amovible, visible sur la 
figure 2 ; il porte un système support (o) pour les 
disques servant à la mesure du moment d’inertie. 
L'ensemble est fixé à une pierre qui repose sur une 
voûte du bâtiment par l'intermédiaire de 4 blocs 
épais de caoutchouc. Ce dispositif assurait au sys- 
tème oscillant une stabilité suffisante : de nuit 
entre 24 et 5 heures les oscillations parasites du 
système au repos étaient d'environ 0,1 mm. 

Le courant d’aimantation est fourni par des 
accumulateurs et est stabilisé par un élément fer- 
hydrogène pour éviter toute variation accidentelle 
du courant pendant les expériences. L’inversion 
automatique du courant est assurée par deux relais 
polarisés Siemens, actionnés par un appareillage 
électronique identique à celui utilisé dans les expé- 
riences antérieures effectuées au laboratoire [4]. 
Un inverseur d’asymétrie I-Il est monté entre 
l'inverse: r automatique etla bobine d’aimantation. 
La position de la fente placée devant la cellule 
photoélectrique concorde avec la position d’équi- 
libre du système. Une vis, micrométrique, permet 
de la déplacer durant l’expérience en cas de dépla- 
cement du zéro. La largeur de son ouverture est de 
quelques dixièmes de millimètres. 


20 MESURES DES MOMENTS. — Celles-ci se font 
par mesure de la f. e. m. induite lors d’une inversion 
de moment dans une bobine p (fig. 1 et fig. 2) de 
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À 
champ connu, qui entoure le système oscillant. 


L’enroulement de cette bobine est porté par un 
tube de plexiglass fileté, sa longueur est de 35,3 cm, 
son diamètre de 2,7 em. Le nombre mesuré de 
spires/cm est de 33,348 + 0,006. La constante 
moyenne de la bobine sur la longueur de la tige est 
alors HJi — 417,63 + 0,07. Cette bobine est fixée 
à une pièce comportant trois fenêtres de passage (q). 


1 


des rayons lumineux. Celle-ci s’applique hermé-) 
tiquement contre lesupport de la tête de suspension 


et, pour permettre l’accès au système oscillant, est w 
amovible vers le bas, guidée dans son déplacement 
par un rail fraisé dans une planchette fixée à la 
pierre supportant toute l’installation (fig. 2). 

L’étalonnage de la mutuelle Hartmann et Braun 
de 0,01 Henry qui a servi d’étalon pour nos mesures 
a été confié au Laboratoire Central des Industries, 
Electriques et àla Physikalisch-Technische Bunde- 
sanstalt ; les deux étalonnages concordent à 0,04 % 
près donnant une valeur moyenne de 0,009958 
Henry. La précision des mesures du moment doit 
être de 0,1 %. 


39 COMPENSATION DES CHAMPS TERRESTRES. — 
Celle-ci est assurée pour la composante NS par un 


1 
{ 


jeu de bobines de Helmholtz de 126 em de dia-, 
mètre. Une bobine rectangulaire (100 X 63 em) 
assure la compensation de la composante EO rési-\ 
duelle. La composante verticale est compensée par, 


deux autres bobines de Helmholtz de 145 em de 
diamètre. Celles-ci peuvent chacune pivoter autour 
d’un axe fixe pour permettre l’accès à l’appareil. 
Enfin quatre petites bobines situées dans un plan 
perpendiculaire à l’axe de la substance à étudier 
assurent la compensation des gradients de 
champ [11]. La compensation des champs hori- 
zontaux est effectuée par l’intermédiaire des indi- 


cations d’un magnétomètre à cadre mobile [4], - 


celle de la composante verticale à l’aide d’une 
bobine tournante, celle des gradients à l’aide d’un 
magnétomètre dipolaire [11]. 

Deux magnétomètres disposés à deux mètres du 
système oscillant indiquent les varitions du-champ 
terrestre en intensité et en direction pendant l’expé- 
rience. Deux graphiques de correspondance entre 
les indications de ces deux magnétomètres et les 
valeurs des champs dans l’espace de mesure, préala- 
blement établis, permettent d'assurer une compen- 
sation continue pendant l'expérience. 


40 MESURE DU MOMENT D’INERTIE. — Elle est 
classique. Deux jeux, chacun de trois disques de 
moment d’inertie différents et de même poids per- 
mettent de déterminer le moment d'inertie à envi- 
ron 0,1 % près ; cette précision a été difficile à 
obtenir. 


99 MESURE DU DÉCRÉMENT LOGARITHMIQUE. — 
Elle a été faite par enregistrement photographique 


No:3 


- des oscillations libres du système, la substance 

- étant maintenue dans le champ de mesure. Une 
moyenne de 10lectures différentes permet d’obtenir 
une précision d'environ 0,1 à 0,2 #. 


60 MESURE DE LA DÉVIATION ANGULAIRE. — 
- L’amplitude de résonance du système est mesurée 
sur une échelle de lecture située à 4 mètres du sys- 
 tème oscillant. L’exiguité du local disponible a 
exigé l’utilisation d’un miroir fixe dont la planéité a 
| été vérifiée par interférométrie. Dans le montage le 
rayon lumineux passe deux fois par la lentille 
disposée devant le miroir tournant, le rayon optique 
à utiliser dans le calcul de la déviation angulaire 
. est donc dans ce cas celui qui va de l’échelle au 
centre optique de la lentille. Sa longueur a dû être 
mesurée de façon indirecte en utilisant les deux 
miroirs plans montés sur le système oscillant : on 
fait tourner celui-ci d’un angle qui est déterminé 
par la déviation observée par lunette sur une échelle 
“située à distance exactement connue de l’un des 
miroirs. Sur l’échelle de mesure on observe en 
même temps le déplacement d’un spot produit par 
la rotation dusecond miroir: la déviation angulaire 
du système étant connue on en déduit la longueur 
du rayon optique qui était de 390,8 cm + 0,1. 


III. Les substances étudiées. — 1° La tige de 
fer étudiée a été usinée dans un cylindre de fer 
Armco ; son poids est de 20,62 grammes. L'analyse 
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effectuée par les services de l’I. R.S. TI. D. indique 
la répartition suivante des impuretés : 


G 0,02 
Si 0,02 
Mn 0,042 © 
S 0,029 
P 0,007 


% 


L’aimantation spécifique à saturation de ce fer 
est en u. e. m./c. g.s. de cor — 222,4 ; à tempé- 
rature ambiante o$,8sx — 217,7 avec une dureté 
magnétique a = 8. 

La tige de nickel a été coulée [4] à partir d’un 
nickel ex carbonyle de l’I. G.-Farben. Son poids 
est de 16,16 grammes. L’analyse de l’I. R. S. I. D. 
a montré que ce nickel était pur à 99,9 % au moins. 
Son aimantation à saturation est de 58,89, à com- 
parer à celle de 58,90 mesurée sur un nickel spec- 
troscopiquement . pur de Johnson-Matthey 
à 99,99 % de pureté ; à température ambiante 
Coo,sssexk — 02,33 avec a — 12: Le point de Curie est 
de 360 °C. 


IV. Les mesures. — Les mesures ont été effec- 
tuées entre 24 heures et 5 heures lorsque les pertur- 
bations mécaniques et électromagnétiques sont les 
plus faibles. Avant chaque expérience il convient 
de réduire au minimum ou de vérifier la compo- 
sante horizontale du moment magnétique de la 
tige. La qualité du réglage obtenu est jugée d’après 


TABLEAU I 


DÉTAILS NUMÉRIQUES DE 2 EXPÉRIENCES GYROMAGNÉTIQUES COMPLÈTES. UNITÉS C. G. S. 


Muem 


: MESURES Caneo M'uem fl 
uem 
œ 1457,9 70,70 11,11 3,196 
1 À 
B 1457,9 70,70 » » 
Fer 
æœ  1460,6 70,83 » » 
22h 
B 1460,6 70,83 » » 
\ 
/ x 474,80 29,38 » 2,936 
1 À 
B 476,79 29,50 ; » 
Nickel 
œ 476,32 29,47 » » 
A 
B 476,32 29,47 » » 


e e 
9 / ve = 
à Fe L E p (2) g (2) 
0,01172 3,289 390,8 276% 1,039 
Il : 7,64 1,039 | 
0,011 72 3,289 » 12766 1,042 
II : 7,67 1,044 
\ 1,039 
0,01128 3,264 » 2789 1,040 
II : 7,88 1,039 
0,01128 3,264 » 15486: 11,038 
11: 286. | 1,039 
0,01149 3,148 » I : 2,90 1,110 
II : 2,945 1127 
0,01134 3,149 » I : 2,865 1,076 
1292289 1,085 
1.080 
0,01116 3,149 " I : 2,93 1,108 
1154 2:79 1,030 
0,01134 3,148 » 1:2,805 1,052 | 
I1:2,800 1,054 / 
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l'amplitude de résonance atteinte par le système 
après une cinquantaine d’inversions du courant de 
la bobine d’aimantation en avance de phase de 900 
sur l’oscillation dans un champ terrestre non com- 
pensé ; après réglages cette amplitude n’atteignait 
que quelques centimètres. Avant chaque mesure la 
substance étudiée a été soigneusement désaimantée 
par application d’un champ alternatif décroissant. 
L'expérience est faite par l'observation de Pampli- 
tude de résonance qui est précédée et suivie de 
mesures du moment magnétique et décrément loga- 
rithmique. Il s’est avéré indispensable de répéter 
l'expérience pour différentes combinaisons de posi- 
tion des inverseurs d’asymétrie et d’orientation 
d’azimuth du système oscillant. 

Si A et B désignent les deux positions de l’inver- 
seur de courant monté sur le système oscillant, 
I et IT celle de l’inverseur suivant l’inverseur auto- 
matique, & et (x + 1800) — 6 deux azimuths du 
système oscillant, on peut schématiser la séquence 
expérimentale d’une nuit par : 


Aol, Acll, A(B)II, A(B)I 


la moyenne de ces quatre expériences éliminant 
l'effet I-II et l'effet x-B. Afin d’éliminer l'effet 
A-B il faut procéder ensuite à une expérience : 
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Comme cette séquence est mesurée au cours 


d’une autre nuit, les conditions expérimentales ne … 


sont pas exactement identiques pour les expé- 
riences À et B : leur moyenne ne peut pré- 
tendre éliminer strictement l'effet ; à fin de plus de 
sûreté il est donc nécessaire de répéter les expé- 
riences. 

Au tableau I figurent tous les chiffres afférents à 


deux mesures complètes, l’une du fer, l’autre du , 
nickel. Toutes les grandeurs sont exprimées en | 


unités c. g. s. Au tableau, en plus des données 
nécessaires au calcul de 9 on a fait figurer l’aiman- 
tation spécifique atteinte par l'échantillon au cours 


des expériences. Dans les tableaux IT et III sont \ 


rassemblés les résultats de l’ensemble des expé- 
riencés. 


V. Discussion des résultats. — 19 FER. — La 
moyenne des quatre expériences complètes est : 


(a) = 1,035 + 0,002 (?). L'examen des résul-_ 


tats partiels figurant au tableau II montre que 


dans le cas du fer les effets d’azimuths et de pertur- . 


bations périodiques asymétriques sont nuls en 
moyenne et faibles pour les mesures individuelles, 
Ce résultat paraît donc très sûr. L'erreur absolue 


N°87 


CT NE DER 


1 


Bal, Boll, B(8)II, B(B)I. maxima provenant des mesures secondaires 
* TABLEAU I | 
RAPPORTS GYPROMAGNÉTIQUES MESURÉS DU FER EN UNITÉS () , 
Mesures Pois «I « II & BI BI 8 aB 
1 A 1 1,039 45039 1,039 4,042 1,044 1,043 1,041 
2 B 1 1,040 1,039 1,039 1,035 L03 081085 1,037 
3B 1 1,035 1,032 1,033 1,037 10872087 1,035 
k B dl 1,044 1,040 1,042 1,031 1,030 1,031 1,036 
5 B il 1,031 1,029 1,030 1,031 4,029 1,030 1,030: 
6 À 1 1,035 1,038 1,036 1,030 1,031 1,030 1,033 
7 A 1 1,033 1,034 1,033 1,033 1,037 1,035 1,034 
8 À 1 1,025 1,028 1,027 1,041 1,044 1,042 1,034 
1:4,085. © «1,035 À 4,036 
T:14,035  B 1,035 B 1,034 
Moyenne finale : 1,035 + 0,002 
(mesures de W, I, l) qui est la même pour toutesles L'examen des résultats partiels figurant au 


expériences et peut donc entâcher le résultat final 
est estimée à environ 0,2 % dans le cas le plus 
défavorable. L’erreur totale maxima serait donc 
de + 0,004. Ce résultat est en parfait accord avec 
ceux de Barnett [5] et de Barnett et Kenny [8], 
avec les mesures antérieures effectuées au labo- 
ratoire P. Weiss [4], [10], ainsi qu’avec les pre- 
mières mesures, de Scott [9], [10]. 


tableau IIT montre que l’effet I-IT est ici encore 
négligeable en moyenne bien que déjà important 
pour une mesure individuelle. Par contre, l’eftet 
d’azimuth moyen et l'effet A-B moyen atteignent 
2 %, du résultat final. II nous semble que ce résultat 
soit en rapport avec la dureté magnétique du nickel 
déjà notable par rapport à celle du fer : les champs 
de dispersion sont peut-être plus importants le long 
| de la tige occasionnant des effets de perturbation 
29 NicxeL. — La moyenne pondérée de 9 expé- (2) Une série de 4 expériences préliminaires avait donné 


riences complètes est : o (£) — 1,078 + 0,003. 1 — 1,037. 
mc, 
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À / TABLEAU III 


; m 
RAPPORTS GYROMAGNÉTIQUES MESURÉS DU NICKEL EN UNITÉS () 
(4 


MESURES Poips I œIl 
1 A 1 1,110 15427 
2 A 4} 1,107 1,124 
3 B 4 1,108 1,030 
kB 1 1,061 4,041 
5B 0,5 1,100 1,078 
6 B - 0,5 1,056 4,038 

7 A 0,5 1,087 AIDE 
8 A 0,5 1,084 1,086 
9 A 0,5 1,086 1,082 
10B 0,5 1,081 1,083 
11 B 0,5 1,088 1,092 
12 B 0,5 1,065 4,070 
13 À 0,5 1,094 1,086 
14 A 0,5 1,072 1,066 
15 A 1 1,113 1,119 
16 À 1 1,120 1,128 
17 B 4 4,104 1,097 
18B 4 41,082 1,080 

I: 1,079 œ 

F5 245077 8 


Moyenne pondérée finale 


lorsqu'ils entrent en interaction avec les compo- 
. santes résiduelles des champs terrestres (effet d’azi- 
muth) ou des circuits d'alimentation de la bobine 
(effets A-B de permutation des connections de la 
bobine). L’erreur absolue maxima due aux mesures 
- secondaires est la même que pour le fer : 0,2 %. 
L’erreur totale maxima possible serait donc d’envi- 
ron + 0,005. 

Des résultats antérieurs connus, seul celui des 
premières expériences corrigées de Scott [10] con- 
corde à peu près avec le nôtre. Sur des déchets de 
tournure du nickel étudiée par cet auteur nous 
avons mesuré une aimantation spécifique à satu- 
ration de 58,83 u. e. m./c. g. s. et un point de 
Curie ® = 3499 à comparer à o — 58,90 et © — 3600 
pour un nickel pur. 

Les mesures antérieures effectuées au laboratoire 
sur le nickel [4] ne sont pas utilisables par suite 

de la non élimination des importants effets d’asy- 
métrie «-f et A-B. 

M. Barnett avait bien voulu mettre à notre dis- 
position la tige de nickel qu’il avait utilisé dans ses 
propres mesures de l'effet Einstein-de Haas. Ce 
nickel possède une aimantation spécifique à satu- 

ration de 71.04 u. e. m./c. g. s. et un point de 
Curie à 285 0C. Le nickel III de Barnett était pur 
_à 98,6 %. Les propriétés des substances étudiées 
par M. Barnett s’écartant sensiblement de celles du 
nickel pur, il semble difficile de retenir sa valeur 
de p pour le nickel. 

Au tableau IV figure l’ensemble des résultats à 
retenir concernant les mesures des rapports gyro- 
magnétiques du fer et du nickel. 


C 


& BI BII B aË 
1,118. 1,076 1,085 1,080 1,099 

,11& 4,076 1,066 1,071 1,092 

,069 1,052 1,054 1,053 41,061 
1,051 41,061 1,053 1,057 1,054 
1,089 1,042 1,032 1,037 1,063 
1,047 1,089 1,064 1,077 1,062 
1,099 1,045 1103 1,074 1,086 
1,085 1,049 1,109 1,079 1,082 
1,08 1,045 1,074 1,060 1,072 
1,082 14,045. 1,033 41,039 1,060 
1,090 1,048 1,061 4,054 1,072 
1,068 : 1,052 1,068 1,060 1,064 
1,090 1,098 1,084 1,091 1,090 
1,069 1,113. 1,102 1,108 1,088 
1,416 1,059 1,071 1,065 1,090 
142% 1,062 1,057 1,060 1,092 
1,400 1,071 - 1,067 1,069 1,084 
1,084 1,069 1,059 1,064 1,073 
1,090 À 1,090 
1,066 PB 1,066 

: 1,078 + 0,003 


D’après les plus récentes mesures de Scott [12] 
les rapports gyromagnétiques du Fe, du Ni et des 
alliages Fe-Ni dépendraient de l’aimantation de ces 


TABLEAU IV 


ENSEMBLE DES RÉSULTATS VALABLES DES MESURES 
DES RAPPORTS GYROMAGNÉTIQUES DU FER ET DU NICKEL. 


me 
UNITÉS (=) : 
(3 


FER NICkREL 
PARNPET TORRES 1,035 
+,003 
SOON ALI OMR LE RER 1,038 1,089 
+,001 +,002 
ME VER AO AP EMRE Er 1,033 
+,004% 
Barnerr et Kenny, 1952 1,037 
+,003 
SCORE Mere ne 1,068 — 1,042 1,110 — 1,093 
+,003 +,001 
Meyer et Brown, 1956 . 1,035 1,078 
+,004 +,005 


substances. À aimantation croissante le rapport 
gyromagnétique du fer diminuerait à partir de 1,068 
pour tendre vers une valeur limite de 1,042 pour les 
aimantations supérieures à 80 u. e. m./cem° ; pour 


, e : à : 
le nickel p (EL passerait dans les mêmes condi- 
c 


tions de 1,110 à 1,093. Dans nos expériences les 
aimantations du fer et du nickel étaient respec- 
tivement d'environ 560 et 260 u. e. m./cm$ donc 
proches de la moitié de la saturation ; c’est-à-dire 
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qu'aux rapports gyromagnétiques que nous avons 
trouvés correspondraient les valeurs les plus basses 
de Scott qui pour le nickel sont supérieures aux 
nôtres d'environ 1 %. Dans les expériences de 
Barnett et Kenny les aimantations utilisées sont 
encore supérieures aux nôtres. Nous n’avons pu 
effectuer des mesures de £ en fonction de l’aiman- 
tation par suite du manque de sensibilité de notre 
appareil que nous avons réalisé en vue de l’obten- 
tion d’une précision élevée. A propos de ces résul- 
tats de Scott il faut encore remarquer que l'effet 
Barnett qui est une expérience apparemment réa- 
lisée dans des conditions correspondant à des 
-champs très faibles donnent des résultats similaires 
à ceux obtenus par mesure de l’effet Einstein-de 
Haas avec des aimantations bien supérieures à 
80 u. e. m./cm°. 
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VI. Conelusions. — De nouvelles mesures de 
rapport gyromagnétique donnent pour le fer 


e() — 1,035 + 0,004 et pour le nickel 
mc < 


€ 
© FE) 1078 + 


l'erreur systématique affectant certaines expé- 
riences antérieures (Meyer, Scott) l’accord est très » 
bon pour les résultats obtenus sur le fer par divers | 


0,005. Après correction de | 


observateurs. Le rapport gyromagnétique du nickel | 


semble être effectivement plus élevé qu’admis . 
jusqu’à présent, les résultats de certaines expé- 
riences antérieures (Barnett, Meyer) n'étant pas 
utilisables pour des raisons qu’on a données. On n’a … 
pu vérifier si, comme l’admet Scott, les rapports 
gyromagnétiques sont fonction de l’aimantation. 


Manuscrit reçu le 10 décembre 1956. 
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UNE FORME SPINORIELLE DES ÉQUATIONS DE L'ÉLECTROMAGNÉTISME 


Par G. BODIOU, 


Faculté des Sciences de Marseille. 


Le présent mémoire est l’application au cas 


général des principes que nous avons déjà énoncés 


lors d’un travail sur l’onde plane [6]. Nous pensons 
que la correspondance entre bivecteurs et spineurs 
gér éralisés du second rang qui conduit à remplacer 
les équations classiques de l’électromagnétisme par 
certames de leurs combinaisons linéaires ne pré- 
sente qu’un intérêt formel. Notre description, en 
termes de spineurs du premier rang, entre plutôt 
dans le cadre général de la méthode de fusion 
de M. L. de Broglie, mais en cherchant à diminuer 
le plus possible le nombre de degrés de liberté de 
l'ensemble des spineurs décrivant l’onde électro- 
magnétique, de telle sorte que la correspondance 
soit très stricte, réciproque si possible, de façon à 
éviter l'apparition de grandeurs et d’équations 
parasites, c’est-à-dire non-électromagnétiques. 
Notations : Un tenseur, X, étant considéré, la 
matrice du 4° ordre qui lui correspond bi-univo- 
quement est notée de même : X ; mais un produit 
de telles matrices est distingué par le signe X 
(multiplié par): X x Y. Le dual, ou supplé- 
mentaire, du tenseur X, est noté D X ; son conjugué 
dans l’espace-temps est noté X’ ; si X est extrait, 
par une loi déterminée, du produit tensoriel des 
spineurs Ÿ, et Ÿ,, nous le noterons X(Y,, F,); 
on sait qu'un tel produit se décompose en un 
ensemble de tenseurs, savoir: un scalaire : 
S(Y;,, F2) ; un vecteur : V(F,, F 2e un bivecteur : 
B(F;, #2) ; un trivecteur : T(F,, YF) ; un quatre- 
vecteur, ou pseudo- -scalaire : Q(F,, F,). Nous 
notons vE le spineur ayant pour matrice :w X #. 


La matrice transposée de X et X. 

Un spineur Ÿ étant donné, nous appelons 
« lacune > associée à ce spineur, et notons F—, le 
spineur symétrique du conjugué, T’, de , par 
rapport à l’origine. F est la somme de deux semi- 
spineurs d'espèces différentes : get © : F = 9 + ©. 
Pour que Ÿ soit confondu avec sa lacune, F — T-, 
il faut et suffit que : © — 9, c’est-à-dire que : 
WW = o + 9. La condition de confusion d’un 
spineur avec sa lacune peut être, sans modifier les 
résultats ci-dessous, définie par : F — + V— (équi- 
valence forte des spineurs [6]), au lieu de l'être par 
l'identité : # —Y— ; sa condition nécessaire et 
suflisante est alors que la « densité de spin », 
T(F, Ÿ”), soit nulle [6]. Il est important de remar- 


_quer que cette notion de « spineur confondu avec 


sa lacune » ne présente pas l’invariance de jauge de 


pr emière espèce : 


RE = 1 (HP) 


F— + Ÿ—, n'implique pas : 
, pOur : LE 1: 


I. Un bivecteur réel est le dual de la somme des 

« densités de moment électromagnétique » de deux 

spineurs confondus avee leurs lacunes. — Démon- 

trons qu'il existe une infinité, à un paramètre, de 

décompositions d’un bivecteur B, donné, réel, sous 
la forme : 

B = B(V,,F,) + B(P,, F) 


avec : 
ir = Gt, Va — "T3 = gs + a» () 

(1) équivaut à : 
B = DB(F,, #5) + DB(F,, Fi), (1 


qui fait apparaître B comme le dual de la somme 
des « densités de moment électromagnétique » 
attachées aux spineurs Ÿ, et F.. 


E, 1. — Pour le prouver, nous utiliserons, d’abord 
un raisonnement géométrique qui évite le choix 
d'une représentation particulière des spineurs et 
souligne le caractère intrinsèque des problèmes 
étudiés. Ce raisonnement utilise une image à trois 
dimensions que nous avons déjà définie [6] : elle 
consiste à faire correspondre au cône isotrope de 
sommet 0 dans l’espace à tee dimensions (M), 
cône dont l’équation est : gf X* XF = (, la qua- 
drique (A) qui admet cette équation pour équation 
homogène dans un espace auxiliaire (m) rapporté 
à des axes orthonormaux. Le carré tensoriel d’un 
semi-spineur est un bivecteur totalement isotrope, 
et tous deux paramètrent un biplan totalement 
isotrope, ou son image dans (m) qui est une géné- 
ratrice de (A), de l’un ou de l’autre système. 

Nous ferons correspondre, à un bivecteur. quel- 
conque de (M), le système de vecteurs glissants, 
défini par ses éléments de réduction en 0, dans (m) : 

A=BS, B—B%, C—B#, 
L=—B"“, MR NEURE 


Pour qu'un bivecteur soit simple il faut et il 
suffit que son système image soit équivalent à un 
vecteur unique. 

Si, de plus, ce bivecteur simple est isotrope, le 
vecteur unique image est tangent à (A). 

La décomposition d’un tel bivecteur simple et 
isotrope en semi-bivecteurs, a pour image la décom- 
position de son vecteur image en deux vecteurs 
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glissant respectivement sur les génératrices de (A) 
qui passent par le contact. Les bivecteurs tota- 
lement isotropes, dont ces deux vecteurs sont les 
images, déterminent, au signe près, les semi- 
spineurs dont ils sont les carrés tensoriels. Ce raison- 
nement fonde la correspondance entre une onde 
électromagnétique plane et un spineur, utilisée 
en [6]. ; 

Mais le bivecteur pour lequel il faut démontrer 
la décomposition (1), avec les conditions (2), est un 
bivecteur réel quelconque, ni simple, ni isotrope. 

Transformons (1). Cette relation s’écrit : 


B = Blo:, ps) + Bo ge) + Blois pi) + Bo2; Pa ). 


Or, By, +), où + est un semi-spineur de 
17e espèce, est un bivecteur anti-dual, et B(9—, ®=), 
où © est. un semi-spineur de 2€ espèce, est auto- 
dual ; d’où : 

DB = — (B(o1, @1) + Box 62)) 
+ (Blps pi) + Boss Pa Je 


La relation (1) équivaut donc aux deux condi- 
tions : 


B— DB — 2(B(o:, ®1) 2e B(? Pe)); 
et : B + DB = 2(B(o;, o1) + Blps pe )), 


qui, étant conjuguées l’une de l’autre, se réduisent 
à une seule, par exemple la première. 

Sous cette nouvelle forme, le problème équivaut 
” à la décomposition du système de vecteurs glissants 
image de (B — DB), en deux vecteurs glissants sur 
des génératrices de même système de la qua- 
drique (A). Ces génératrices seront les images des 
biplans totalement isotropes des semi-spineurs w;, 
et ®2. Or on sait que la décomposition d’un système 
de vecteurs glissants en deux vecteurs est déter- 
minée par le support de l’un d’eux, ici par la 
donnée d’une génératrice. Or (B — DB) est anti- 
dual, le premier bivecteur de décomposition, corres- 
pondant à un biplan dont l’image est une géné- 
ratrice, est aussi anti-dual, donc le second bivec- 
teur est anti-dual, donc son biplan a bien pour 
image une génératrice de même système que la 
première. 

Il y a donc une infinité à un paramètre de telles 
décompositions, ce paramètre étant celui de la 
génératrice arbitraire. 

Pour achever de déterminer la décomposition, 
il suffirait d'imposer une seconde relation invo- 
lutive entre ces deux génératrices ; comme elles 
appartiennent à un même système, l’image d’une 
correspondance par retournement d’espace-temps 
entre les biplans de +, et w, serait une telle relation. 


I, 2. — Pour donner une preuve analytique, 
nous choisirons une représentation élémentaire des 
spineurs qui mette commodément en évidence les 
semi-spineurs des deux espèces : par exemple celle 
de E. Cartan [1]. 
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Posons : 
®1 = [lai il, 
avec b, = u,a,, ou u, est le paramètre du biplan 


totalement isotrope de +, ou de sa génératrice 4 


image ; de même: 6, = [4 b]|, avec: bs = Us Qo: 
La décomposition de (B— DB) écrite plus haut 
équivaut, [1, $ 154], au système : 


CE CRT 
ai ui, +-di us — P 
2 DES 
di + a = Q 
2 2 7, 
Gi Uy + Go Us — 


où P, Q, R, sont déterminés par B. 

Si l’on exclut le cas où (4, — u,) = 0, quicorres- 
pond à : @: = kw, les deux dernières équations 
forment un système de Cramer aux inconnues d 
et a? ; la condition de leur compatibilité avec la 
première s’écrit : 

P— Ru + u) + Quiuw = 0, 


qui exprime la correspondance involutive entre w; 
et w, prévue géométriquement. Ayant choisi %, et 
u, liés par cette condition, le système de Cramer 
déterminera a? et a2, c’est-à-dire o, et o,, au signe 
près ;1l y a donc bien une infinité, à un paramètre, 
u, par exemple, de solutions pour (1). 


I, 3. — Remarques : Dans notre mémoire [6] 
nous avons utilisé la décomposition suivante de 
l’onde plane «complexe » : 

2Ti 
B(M) — B, exp Ée OP.OM) 
2 B(Yoexp + OP.OM, Y exp + OP.OM) ; 


où B, est un bivecteur simple, isotrope, réel, donné, 
LOU PIE ENS. 


On en déduit, pour l’onde plane réelle, 
B; cos 27 OP.OM, la décomposition : 


B; cos E OP.LM 


a B(Y cos? OP.OM, Ÿ cos à OP.OM) 


+ B (ri sin + OP.OM, Yisin 7 OP. OM) 


hk 
= BF, Pit BF, 
avec : 


LR cos 5 OP.OM ER Th. | 
Y, — oi sin + OP.OM = +%, 


qui est du type (1), F, et Ÿ, étant « fortement 
équivalents » à Let Z. 


II. Forme spinorielle des équations de Maxwell- 
Lorentz. — II, 1. — Écrivons les équations de 


 No3 


Maxwell-Lorentz, dans le cas 6ù le courant est J : 
di Bi = Ji 
dt DB = 0. 


Elles équivalent à l’équation : 
di(B — DB)5 — Ji, 
jointe à l'équation conjuguée. C'est-à-dire, en uti- 
lisant la décomposition de (B — DB) : 
dE Bulons o1) + D Bulo» Pa) — Jul2. 
Mais :: 
_Bule, = X OX NX vi X @ 


. où les y sont les matrices des vecteurs de la base ; 


et : 


d’où : 
di Bule, g) = 2e X C X NX vi X div, 
avec 1 =£ de soit: 
D Buy, D=2eXCXnXnXHeE. 
:— 26 X C X do, 
l'indice z, dans le premier terme du second membre, 


prenant maintenant ses quatre valeurs. 
Les équations de l’électromagnétisme s’écrivent 


. donc: 


Pi X COX Yi X Va + Pa X OX VI X Vs 
— (os X C°X Doi + ® X C X D pe) = Jifk. 


Les o étant des semi-spineurs de première espèce, 
elles s’écrivent aussi, avec K = y, X ÿ2 X Ys X Ya: 


D X COX NX Æ X Yo + Po X C X M X K X Vo 
— (1 X CO X Dos + D X C X D) — Jylé. 
Ve, yo) = ox C X n x Æ x ve, 

Sp, D) = ? X C XD: 


Or: 


formons alors une combinaison linéaire des quatre 
équations précédentes en prenant pour coefficients 
des dz, arbitraires d’univers, coordonnées d’un dM 
arbitraire de cet univers ; le système des quatre 
équations précédentes équivaudra à l'équation 
unique obtenue en annulant cette forme, qui 
s'écrit : 
AM.(V(o1 Ver) + Vips 92) 

.  (S(o1 des) + S(p2; dp2)) = AM.JJ4, (3) 


dans laquelle le . marque des produits scalaires de 
vecteurs, et où dy, et do, sont les différentielles 
semi-spinorielles correspondant à dW. 


IT, 2. — Utilisons, maintenant, l’indétermination 
de la décomposition (1), pour imposer aux spi- 
neurs ®., et w, une condition supplémentaire des- 
tinée à simplifier (3). Nous choisissons la condition : 


S(p1; dpy) + (Pr dope) = 0, (4) 
Cette condition (4) est satisfaite, par exemple, 


si l’on impose à w, et w, de se correspondre dans un 
retournement fixe arbitraire. Ce retournement 


se traduit par l’équation matricielle p, = b X :, 
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où b est une matrice bivectorielle telle que : 
b? 1 ; on sait que : 
bXC——C X b ; on trouve alors : 
S(p2; dpa) = — p1 X OC X D? X dpi = — S(p1 AP), 
qui est bien (4). 

Cette condition (4) est aussi satisfaite si o, et o 
sont de la forme : 

®1 = 3 fil M), Pr P2 f(M) 

où pet ©$ sont des semi-“spineurs fixes, et f.( M) et 
f( M) des fonctions scalaires arbitraires dans l’uni- 
vers ; les ondes semi-spinorielles utilisées en I,3 
dans la décomposition d’une onde plane réelle sont 
de cette forme. On trouve, dans ce cas : 


9 dfi ; d'où : S(pr, ds) = f1 df115(p9, p!) 
S(ot, pi) = 0; 


dp1 = 
OT : 
d’où (4). | $ 
Vérifions qu’il existe, lorsque B est un champ 
bivectoriel donné, une infinité à un paramètre de 
champs semi-spinoriels, solutions de (1) et de 
l'équation aux différentielles totales (4) : 
En conservant la représentation de E. Cartan, 
déjà utilisée, l’ensemble de (1) et de (4) s’écrit : 
aÿ u? + au? — P 
+= 
dur + Ga — À 
ai du, + a? dus = 0, 
dont les conditions de compatibilité pour les incon- 
nues af et aÿs’écrivent : 
P— Riu + uw) + Quiu = 0 
R(du, — du) + Q{u, du, — us dus) = 0. 


La première de ces équations détermine w, en 
fonction de w,, et la seconde, devient alors une 
équation aux différentielles totales pour %,, qui 
2% L 
s'écrit : 


ME Pr 2 , Ru / 
AL RP .dR’ + RP pe 
avec : PE=P)0; RERO; 


et dont la condition d’intégralité est : 
(dP' AdR }(u?—2R'uw + P') = 0. 
L’annulation du second facteur étant exclue par 
la condition w,  u,, elle se réduit à : 
LPENGRE 0 est PER"); 
soit encore : 
FX, E)Z)=0, (4°) 
F étant une fonction homogène arbitraire des 
coordonnées, X, Y, Z, du semi-bivecteur B-DB. 
Ce qui s’applique à des champs satisfaisant à (4). 


Il, 3. — Forme spinorielle des équations de 
l’électromagnétisme : 
L’équation (3) se réduit, par (4), à : 
Vu Von + Vos V?2) = J14. 
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Nous voulons montrer l’équivalence de cette 
dernière au système d'équations matricielles spino- 
rielles, symétriques en +, eb æ; : 

85(p1; Pa) X V X Pi = J X (5) 
85 (pa; Pa) X Ÿ X Pe = J X Pi. 


Nous poserons : 
V(ou VP1) = Vs V(p2 VP2) one 
on a donc : V, + Vo: = j. ne 
V, étant le vecteur extrait du produit tensoriel 
des semi-spineurs +, et y ®., estisotrope et commun 
aux biplans totalement isotropes de ces semi- 
spineurs ; de même pour V, ; on a donc : 
LR Rens 


TJ& = 7j; 


V, X vi = 0, d’où: 7 X o1 = Va X 1. 


Mais V, X ©, est la matrice d’un semi-spineur 
de seconde espèce dont le biplan contient V, ; ce 
biplan est donc celui de yo, ; et, par suite : 


IX p=PX VX Pp et: 1 X Pr XV X 15 
p et g étant des scalaires ; ces équations sont de la 
forme (5) ; il ne reste qu’à déterminer p et q. 

Pour le faire substituons les expressions ci-dessus 
de vo. et de vo, dans l’équation : V, + V, = 7; 
on trouve : 

= 1 Etre L 
p1 X C X Gr ET EUNEN x) Kobe 
P 
ce qui équivaut à : 
Q 4 à 11 ë 
DC Cr x x) x a 


pour tout vecteur X, arbitraire. 
Il suffit qu’elle soit vérifiée pour tout vecteur 
orthogonal à 7 et pour 7 lui-même : 
pour X orthogonal à j elle s’écrit : 
LES ) 
- += XX CN AO ERN pe 0; 
Ë F : @1 ] Pa 
d'où : pr g= "0; 
pour ] elle s’écrit : 

SI X Gi X OX @ù = à 
d’où : - 
q = 2p1 X C X Pa = 21, Pa), 
ce qui achève de démontrer les (5), compte tenu de 
l’antisymétrie de S(p1, po). 

Réciproquement : (5) implique que 7 est commun 
aux 2 triplans déterminés, d’une part par les 
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biplans de o, et de yo», d'autre part par les biplans 
de w, et de wo, (voir figure ci-contre représentant 


FTes 1. 


les images, dans (m), des éléments considérés). Îl en 
résulte que 7 est dans le biplan déterminé par V, 
et V;, donc : j = x V, + y V,; revenant-aux(S) 
ON Prouve que : &4 = y. — 14 


Conclusion. — Za solution générale des équations 
de Maxwell-Lorentz et de (4°) se décrit donc sous la 
forme (1) ou (1'), à partir de deux semi-spineurs, 
®, et +, solutions de (4) et de (5). 

(5) équivaut à un système de quatre équations 
aux dérivées partielles, du premier ordre, à quatre 
fonctions inconnues qui sont les composantes de w; 
et ®: ; ces équations ne sont pas linéaires par suite 
de la présence du facteur S(w;, ®). 


IT, 4. — Dans le vide, où J = 0, le système (5) se 
réduit à : 


X p1 = 0 TR 
Fe * . SEE F1 
qui équivaut, par (2) à : 
VX V;= 0 dog) 
VX, =0 El 


Chacune des équations (5”) est une équation de 
Dirac pour un corpuscule à masse nulle. 
(4) équivaut à : 


S( UPR d d) + S( Vo, d UP) ir 0, 


Conclusion. — Par (1'), {oute onde électromagné- 
tique dans le vide satisfaisant à (4') peut donc se décrire 
(d’une simple infinité de façons) comme duale de la 
somme des densités de moment électromagnétique de 
deux corpuscules de Dirac, à masse nulle (ou d'un 
ensemble de tels couples, puisque l'onde n’est pas 
normée), confondus avec leurs lacunes (ou ce qui 
revient au même : à densité de spin nulle) liés par (4). 


Manuscrit reçu le 8 décembre 1956, 
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1. Le problème fondamental de la mécanique 
statistique des phénomènes irréversibles. — Le 
problème fondamental de la mécanique statistique 
des phénomènes irréversibles consiste à éclaircir la 
contradiction apparente entre la réversibilité par 
rapport au temps de l’équation de Schrôdinger qui 
régit l’évolution des systèmes matériels et l’irré- 
versibilité du comportement macroscopique de ces 
systèmes, constatée expérimentalement. Le pro- 
blème fondamental est de prouver, en partant de 
l’équation de Schrôdinger, que ces systèmes macro- 
scopiques tendent vers un état d'équilibre. 

La première étude importante sur l’évolution 
irréversible d’un système matériel macroscopique 
isolé S vers son état d'équilibre en mécanique quan- 
tique est düe à W. Pauli [1]. Pauli écrit l’hamil- 
tonien J€ du système S sous la forme 


R=K+V (1) 


K est l’hamiltonien non perturbé, V est une per- 
turbation induisant des transitions entre les 
niveaux M, n, ... , p, ... de K. L'énergie du 
système S reste constante, À et V sont indé- 
pendants du temps. Si le système S considéré est un 
gaz par exemple, XÀ représentera l’énergie ciné- 
tique et interne de toutes les molécules gazeuses, 
et V tiendra compte des interactions (chocs) entre 
ces molécules. Le problème fondamental est alors 
de déduire de l’équation de Schrüdinger 
d®(t) 


Da ner q() (2) 


les équations «irréversibles par rapport au temps » 


D = 5 Wap Pal) + 2 Wom Pol (3 
qui donnent une description dynamique de la rela- 
xation. P,(t) désigne la probabilité pour que le 
système S soit à l’instant { dans l’état m, W,,, est 
la probabilité de transition de l’état m vers l’état p 
par unité de temps. Suivant les principes de la 
mécanique statistique, Pauli introduit un ensemble 
de Gibbs de systèmes S, chacun des systèmes S de 
l’ensemble étant représenté par une fonction ®(1) 
particulière. Mais Pauli ne parvient à déduire 
l’équation (3) à partir de l’équation (2) que moyen- 
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nant l’hypothèse que les phases des amplitudes 
complexes de probabilité des diverses fonctions ®(6) 
dans la représentation m, n, ..., p, sont 
réparties au hasard à l'instant initial { = 0 et à 
toute une série d’instants £1, to, ... t;, ... très 
nombreux, intermédiaires entre l’instant initial 
t — 0 et les temps de relaxation 7, Cette hypo- 
thèse revient à admettre que 


(mlo(t)|p) = Smplmle(t;)|m) (4) 


c’est-à-dire que la matrice de densité relative à 
l’ensemble de Gibbs de systèmes S est diagonale 
aux divers instant t, (7 = 0, 1,2, ...). Mais. 
l'équation (4) était vérifiée à l’instant t;, elle ne le 
serait plus à l'instant f,,, comme on le voit aisément 
en intégrant l’équation d’évolution de l’opérateur 
densité 


EU) À (sol) — p(16€). (5) 

L'hypothèse des phases réparties au hasard, 
invoquée par Pauli pour déduire (3) à partir de (2) 
ne peut.donc être considérée que comme une hypo- 
thèse ad hoc, au surplus en contradiction avec les 
lois mêmes de la mécanique quantique. 

Un progrès important fut réalisé par L. Van 
Hove en 1955 [2]. Van Hove remarque qu’il existe 
des perturbations V qui certainement n’engendrent 
pas des phénomènes de relaxation : un exemple 
d’une telle perturbation est constitué par l’inter- 
action entre le champ de l’électron libre et le champ 
de radiation électromagnétique en  électro- 
dynamique quantique. Les perturbations V qui 
déterminent une évolution irréversible de S vers un 
état d'équilibre doivent par conséquent vérifier 
certaines conditions, qui naturellement ne peuvent 
pas être imposées ad hoc, mais qui au contraire 
doivent avoir une origine physique. Van Hove 
parvient à montrer, sans faire un usage répété de 
lhypothèse des phases réparties au hasard, que les 
perturbations V qui satisfont à la relation 


(m|VAV|m) = S(Em — Em) fa(m) + (m'|g(4)lm) (6) 


(où À est un opérateur diagonal dans la repré- 
sentation ", m', ... p,...eboù Ey, Em, .:. En, .. 
désignent les valeurs propres correspondantes 
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de Æ) engendrent effectivement des phénomènes 
de relaxation à partir de certains états initiaux de 
non-équilibre. Van Hove affirme que la condi- 
tion (6) est réalisée par exemple pour l'interaction 
entre phonons dans les solides (conduction de la 
chaleur). Il est vraisemblable qu'il existe des per- 
turbations V qui n’obéissent pas à la condition (6) 
et qui pourtant engendrent également des effets 
dissipatifs irréversibles. 

Enfin, en 1956, l’auteur a développé la théorie 
de la relaxation par mouvements moléculaires 
internes (par exemple, semi-rotations du chlore 
dans des cristaux moléculaires renfermant des grou- 
pements CCI.) en résonance quadrupolaire, théorie 
dans laquelle on prouve le retour des spins 
nucléaires à l'équilibre thermique avec le réseau 
sans faire usage de l’artifice des phases réparties au 
hasard. Un résumé suceinct de cette théorie a été 
publié [3] (1). 


2. Les origines de la relaxation en résonance 
nucléaire. — Le présent travail est avant tout 
consacré au. problème fondamental de l'irréver- 
sibihté en résonance nucléaire, mais la méthode qui 
sera développée pour étudier ces problèmes pourra 
permettre de traiter d’autres questions. 

En résonance nucléaire magnétique, les spins 
nucléaires sont placés dans un champ magné- 


— 
tique H, constant et acquièrent ainsi une énergie 
(?) 


où y est le rapport gyromagnétique de ces spins 
— 


et I, l'opérateur vectoriel spin correspondant 
au j°me spin (7 = 1, 2, ..., N). N est le nombre 
total de spins. Un champ magnétique de radio- 
fréquence excite des transitions entre les niveaux 
de À et permet ainsi d’écarter les spins de leur 
état d'équilibre. Le retour des spins à leur état 
d'équilibre s'effectue grâce au fait que ces spins 
sont couplés à leur environnement moléculaire. 
S'il s’agit de la résonance nucléaire magnétique de 
spins dans un liquide par exemple, les spins 
nucléaires sont couplés au mouvement brownien 
des molécules par l’interaction dipolaire 


F — 


yat re Ga 0) Gi) | (8) 
T% 


È 
sx r& 


(2) Contrairement à la théorie de la relaxation exposée 
ici et contrairement à la théorie de Van Hove dans les- 
quelles la cause de la relaxation est un couplage « faible» Y, 
la théorie de la relaxation par mouvements moléculaires 
internes traite (pour la première fois à la connaissance de 
l’auteur) un cas où le couplage n’est pas faible et où par 
conséquent une méthode basée sur le développement de 
l'opérateur d'évolution en série de puissances de V n’est 
pas valable. 
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ou par l'interaction quadrupolaire 


—2 
px > of 
i B=2 


ou encore par d’autres types d'interaction. r}y; repré-\ 


: : LS 
sente le vecteur joignant les spins 7 et X, rx = ral. 
Les DS (B — — 2, 


santes appelées_F?* dans [4] (égalités (3)). Les I 
= — 2, 
irréductible d'ordre deux construit à partir des 


composantes Z.,;, L,,, L; du spin 7, et qui sont éga- 


lement définies dans [4] (égalités (4) et (5)). En 
résonance nucléaire, on pourra considérer le réseau, 
c’est-à-dire le système moléculaire auquel sont 
couplés les spins, comme un thermostat dont la 


capacité calorifique est infinie par rapport à celle 
des spins (sauf peut-être à très basse température) . 


et qui n’est donc pas affecté par la résonance de 


ces spirs et par l’apport d'énergie qu’il reçoit de ces … 


spins. 


y 


3. Relaxation de systèmes isolés et relaxation | 
de systèmes en contact avec un thermostat, — Pour » 
étudier la relaxation spin-réseau, on peut se placer 


à deux points de vue. 
A) On peut considérer le système { spi:s+réseau } 
comme un système isolé. Le réseau possède une 


énergie 4. Dans le cas de la résonance nucléaire. 

dans un liquide, Æ, représenterait l'énergie des 
molécules (du liquide) en mouvement. L’hamil- 
tonien du système isolé { spins + réseau } est alors M 


(10) 


=K+HR+V 


où À est donné par (7) et V par (8) cu (9) par 
exemple. #€ ne dépend pas du temps. Si l’on admet 
que le couplage V satisfait à une condition du 


type (6), on pourrait adapter l'argumentation de 


Van Hove au cas de l’hamiltonien (10) et démon- 
trer, en suivant une méthode analogue à celle de 
Wangsness et Bloch{5], l’irréversibilité des échanges 
d'énergie entre les spins et le réseau. Brout et 
Prigogine [18], dans un article consacré à la trans- 
cription au cas de la mécanique classique de la 
théorie quantique de Van Hove (en s'appuyant 
pour transerire (6) sur le principe de correspon- 
dance), ont donné une interprétation physique 
de (6) qui est probablement correcte dans son 
essence. Cependant, en résonance nucléaire dans 
les liquides, principal domaine d’application de la 
théorie développée ici, la forme exacte de Æ, est 
inconnue et il est par suite diflicile de vérifier diree- 
tement si la condition (6) est réellement satisfaite. 
C’est l’une des raisons pour lesquelles nous n’adop- 
terons pas ici le point de vue A, ni la condition (6). 
Ce point de vue sera néanmoins exposé dans une 
autre publication de l’auteur. 


(9). 


, 2) sont les composantes. 
du tenseur (irréductible d'ordre deux) gradient de, 
champ électrique à Pendroit du noyau 7, conpes 


. 2) sont les composantes du tenseur 


DT", JS ETES PEN LITE ct 


B) Dans la relaxation spin-réseau, ce sont les 
Spins qui se relaxent. Le réseau, lui, est un système 
stationnaire en équilibre thermique (thermostat). 
Au lieu de porter notre attention sur le système 
isolé { spins + réseau }, nous pourrons nous attacher 
… au système { spins } seulement. Ce système { spins } 
n’est plus alors un système isolé, son énergie ne 
reste pas constante, c’est un système en contact 
avec un thermostat. L’hamiltonien de ce système 
_ { spins } en contact avec un thermostat sera 
3€(t) = K + V(i) (11) 
où le couplage spin-réseau V est cette fois une 
fonction explicite du temps. Alors que Van Hove 
traitait la relaxation des systèmes isolés, nous 
adopterons ici le point de vue B et traiterons la 
relaxation d’un système en contact avec un ther- 
mostat. Au paragraphe suivant, nous énonçerons 
les conditions auxquelles devra satisfaire la fonc- 
tion V(t), et qui joueront ici le rôle que la condi- 
tion (6) joue dans la démonstration de Van Hove. 

Un examen critique du point de vue B sera 
présenté en paragraphe 12. Néanmoins, remar- 
quons tout de suite que dans le point de vue A (10), 
le système { spins + réseau } est considéré comme 
un système quantique, alors que dans le point de 
vue B, le système { spins } est bien quantique, mais 
le réseau intervient « classiquement ». Ce point de 
vue semi-quantique B est valable en| résonance 
nucléaire parce que les énergies fom:, = Em — En 
mises en jeu sont très faibles devant AT (T étant la 
température du réseau). 


4. Perturbations fonctions aléatoires station- 
paires du temps. — Dans le présent travail, nous 
développerons la théorie de la relaxation nucléaire 
dans les liquides. En ce cas, le point de vue B est 
beaucoup plus avantageux pour le calcul effectif 
des temps de relaxation car dans le cas d’un liquide, 
. la forme de Æ7, dans (10) n’est pas connue en 
général, alors que l’on dispose de l’image « clas- 
sique » du mouvement brownien. Dans le point de 
vue B, V(t) sera par conséquent une fonction 
aléatoire stationnaire du temps. 

Considérons donc AN spins nucléaires de rapport 
gyromagnétique y dans un liquide. L’énergie 
du jeme spin (j = 1, 2, ... AN) dans le champ 

— 


magnétique H, sera 


> 
HO = — yr H, L. (12) 

Nous désignerons par V@(f) l'énergie d’inter- 
action du j°%e spin avec le réseau, de sorte que 
lhamiltonien total du système des NW spins consi- 
dérés sera : 


N Je à N 
(tn = D (1) — © HG + E VA(t. (13) 
j=1 fi] î=1 
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Nous supposerons que chacun des N spins 
nucléaires interagit avec le réseau indépendam- 
ment des autres et « se relaxe » par conséquent 
indépendamment des autres. Cette hypothèse, éga- 
lement faite par Wangsness et Bloch, est vérifiée 
lorsque VW(4) est le couplage quadrupolaire (9) ou 
encore lorsque VÜ(f) représente l’interaction directe 
du spin 7 avec un ion paramagnétique dans le 
liquide. Dans tous ces cas, VW({) ne contient que les 
variables de spin du J°me spin, c’est-à-dire les 

— 


composantes de I; et ne contient pas les compo- 


santes des autres variables de spin E Notre hypo- 
thèse n’est pas vérifiée strictement lorsque le méca- 
nisme de relaxation prédominant est le couplage 
spin-spin (8). En ce cas, le comportement d’un spin 
donné dépend de celui des spins immédiatement 
voisins, et VÜW(?) porte sur les variables de spin de 
tous les V spins du liquide. Le problème général 
de la relaxation dans une assemblée de N spins 
nucléaires tous couplés entre eux par l’interaction 
spin-spin (8) est un problème difficile qui ne sera 
pas abordé ici. Par contre, on pourra généraliser 
notre théorie au cas d’une assemblée de spins dans 
laquelle chaque spin 7 + f n’interagit qu'avec les 
Sp JL je (0<T< 8; I= TL 
g + 1, 2g + 1, ...). Telle est la relaxation intra- 
moléculaire des protons dans l’eau, le spin de 
chaque proton n’interagissant qu'avec le spin de 
l’autre proton de la molécule d’eau. Dans le cas 
général d’un système de N spins couplés par l’inter- 
action spin-spin (8), il faudra pour appliquer notre 
théorie calculer le champ magnétique fluctuant 
produit par les V—1 moments magnétiques 
122,09 47 +1, :.., N à l'endroit où, se 
trouve le spin 7 ; ceci pourra se faire par exemple 
en remplaçant dans l'expression de VÜ(i) les 

—> 


> —> > > 


variables quantiques I;, I, ..., 15, 21, ... 1y 


par des variables classiques! I, (0), fonctions 
explicites du temps. Cette approximation du 
« champ local » sera valable pour des spins nuclé- 
aires suffisamment dilués dans un liquide ou pour 
des spins nucléaires dans un gaz. 

Soient mt, nù, ..., p®, ... les états propres 
de A, Désignons par +, le temps de corrélation 
de la fonction aléatoire V({). Pour des temps {> =, 
on aura 


(m| (en) (plW(0)la) — (mlV(s)1n).(p1V(0)lg) (14) 
où les moyennes portent sur les V spins nucléaires 


du liquide, c’est-à-dire que 


CRETE x 
PT = À (mt2| PU (s)1nt) 

US" 
par exemple. Tout au long de ce mémoire, la valeur 
moyenne d’une quantité X, indiquée en surlignant 
cette quantité, devra être prise sur les V spins 
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nucléaires du liquide : 
_ 1 : : 
RS MP D A ET (D) NE xX(N 
N ( ) 


(5) 


(14) pourra s’écrire symboliquement 
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V{t) V{0) = V{i).V(0) pour H > x (16) 


N°3 À 
y 


où 4 peut naturellement être positif ou négatif. Plus 


généralement, on aura 


la) Pl 02 Pl Fle) Ve) ele) 


= Via) Plon) 


à condition que tous les intervalles 
El Re 


soient grands devant % 

La condition que chaque spin interagit avec 
le thermostat indépendamment des autres et 
la condition (17) jouent dans le point de vue B 
(relaxation d’un système couplé avec un thermostat 
traitement semi-quantique) le même rôle que la 
condition de Van Hove (6) dans le point de vue A 
(relaxation d’un système isolé, traitement quan- 
tique). La condition (17), jointe à la condition de 
stationnarité (18) du thermostat, permet de démon- 
trer rigoureusement l’évolution des N spins du 
liquide vers l’équilibre thermique avec le réseau. 

La signification physique de la condition (17) 
est claire. Considérons tous les spins nucléaires 
j, k, ..., lpour lesquels, à instant { = 0, on a 


(mD| VE (0)[n D) = (MB | (0)|n®) = 
= (mO|V@(0)[n®) 


quels que soient m et n, soit symboliquement 
V&(0) — 7 (0) 


— p&(0). 


Les environnements de ces spins 7, #, ..., [seront 
en général différents. Entre les instants 0 et t, la 
molécule 7 (par exemple) subira par consé- 
quent des collisions qui ne seront pas iden- 
tiques à celles subies par les molécules 4, 
de sorte que pour des temps { > +, les pertur- 
bations V@(t), V(r), V@(#) n'auront plus 
aucune corrélation entre elles et ne « se rappel- 
leront » plus qu’elles avaient des valeurs identiques 
à l'instant initial & — 0. La condition (17) n’est 
donc pas ad hoc, elle doit être considérée comme 
une propriété physique du thermostat, c’est-à-dire 
ici comme une propriété physique des liquides 
formés d’un très grand nombre de molécules. La 
démonstration de (17) relève de la théorie molé- 
culaire des liquides et n’a pas à être abordée ici. 
On remarquera enfin que la condition (17) est très 
générale : elle n’implique nullement que les pro- 
cessus au sein du thermostat soient de Markoff par 
exemple. 

A la condition (17), nous devons adjoindre la 
condition 


Vin — 1) Vita —t)...-V(t, — 1) 
V{t) 


= Vi) Vin)... (18) 


lALPSIR ALP A DEPE 
Vlan) : Vitps) Vip): Vie): -.. 


V{trn) 


DATA AC 7) 2) 


stationnaire, invariant dans les translations paral- 
lèles à l’axe des temps. Dans (18), { est arbitraire, 
positif cu négatif. De (18), il résulte que la fonction 


Ê 


J 
: 


PS Ms VAI, 


È 4 
qui exprime le fait que lethermostat est un système \ 


qui figure à droite dans (18) ne dépend que des 


différences 4, — 1,2, ln "ba TOR 


Nous pourrons supposer que 


V{t) = V(0) = 0. (19) 


En effet, si V(0) 0, nous pourrons l’inclure 


dans À et écrire 


d@(t) HD + V(0) +[VA() — 7]. (20) 


La perturbation sera V4%(4) — V(t) dont la valeur 
moyenne est évidemment nulle. 


L'égalité V(t) — V(0) 
de (18). 

Nous admettrons que la fonction (18) est ure 
fonction continue de #,, t,_, ..., t,. Il résulte alors 
de (17) et (19) que 

V0 (0) — 0 pour 


RH > <e (21) 


Plus généralement, la fonction 
© Vtu) Vin) 


tend vers zéro dès que tous les intervalles (4, —14 |, 


LAS TE A e ne trs fé — EE 
lt — 1, sont grands devant r.. 


V(t:) 


sl 


5. La réversibilité de l’équation de Schrôdinger. - 


— I] nous faut maintenant intégrer l'équation de 
Schrüdinger relative au système des AW spins 


ja D = act) @() (22) 
dt 

où € (+) est l’hamiltonien (13). L'intégration exacte 
de l’équation (22) nécessite l'emploi d’un appareil 
mathématique assez compliqué. En raison de 
lPimportance de la question traitée ici, et pour 
permettre au lecteur de suivre le raisonnement sans 
trop de difficultés, nous exposerons l'argumentation 
et les calculs en détail. L’équation (22) possède des 
solutions de la forme 


D(t) = FU)(+) PO (s) . 2: FO)... Fr) 


où W(t) représente la fonction d’onde relative 
au jee spin et vérifie l'équation 


iñ RS — I (1) FA (1) (23) 


est une forme particulière | 


N°3 


en vertu de l'hypothèse faite que 4€(t) ne porte 
que sur les variables du spin 7. Pour une assemblée 
formée d’un très grand nombre AN de spins (dans 
toutes les expériences de résonance nucléaire, 
N > 102), la réversibilité par rapport au temps 
. de l'équation de Schrüdinger (22) résulte du fait 
que le thermostat (que nous désignerons par ©) est 
un système stationnaire et qui obéit lui-même à une 
équation d'évolution réversible par rapport au 
temps. 

En effet, inversons, à l’instant initial & — 0, le 
sens d'évolution du thermostat O — c’està-dire 
remplaçons, à l'instant { — 0, les vitesses (0) de 
toutes les molécules du liquide par leurs opposées 
— v®(0) sans changer leur position, etc... — Si 
l’hamiltonien 7, qui représente le thermostat © 
est réel, cela revient à remplacer la fonction 
d’onde U(t) du thermostat © par U*(— 1). Nous 
obtenons ainsi un thermostat inversé O,, évidem- 
ment identique au thermostat initialement consi- 
déré O, et pour lequel par conséquent 


(24) 

LV) V(tn_1) dy V(t) Jer = VU) V (tn) pires V{t)le 

où le membre de gauche est calculé pour le ther- 

mostat O, et le membre de droite pour ©. En vertu 
du fait que l’équation 

. dU(t) 

(or ren 


qui détermine l’évolution du thermostat © est 
réversible par rapport au temps, on aura 


[VO (ter = [VO(— ile (25) 


parce que V ne porte que sur des variables réelles 
d’espace du thermostat (en l’occurrence, les 
variables qui fixent la position des V molécules du 
liquide) et non sur les dérivées par rapport au 
temps de ces variables. Si alors on effectue dans 
l’équation (23) la transformation 


= HR U(t) 


FD 


| Enar d 
il vient l'équation 


. dFO(— à) 


ET = (200 (ber PA (— 1) 


(26) 
où [2% (t)lor, calculé pour le thermostat O, «iden- 
tique » à ©, est une fonction aléatoire stationnaire 
«identique » à [4% (+)], qui figurait dans (23). C’est 
en cela que consiste la réversibilité par rapport au 
temps de l’équation de Schrüdinger (22) (?). On 
voit sur (23) et (26) qu'il n'existe pas a priori une 


(2) On peut obtenir l’invariance de l’hamiltonien 4€(t) 
dans le changement de t en — + de la manière suivante : 
nous adjoignons à nos ! spins 7 initialement considérés 
(j = 1, 2,..., N) un autre système de AN spins numérotés 
de N +1à2N et tels que : 


JUN (4) = HU (— t) (= 1 N). 
Le système des 2N spins ainsi obtenu possède les mêmes 
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direction privilégiée d'évolution du système des N 
spins. Notons enfin qu’il résulte de (24) et (25) que 
Vita) Vins) Ve) = Ve tn) Vi tmx) ce V(—) 


où les deux membres portent sur le même ther- 
mostat ©. 


6. Les conditions initiales suffisantes pour qu’il 
y ait relaxation. — L’équation (22), réversible par 
rapport au temps, et les équations (3), irréversibles 
par rapport au temps, ne peuvent évidemment pas 
être équivalentes entre elles. Pour pouvoir 
déduire (3) de (22), il faut que la fonction 
d'onde ©®(t) satisfasse à certaines conditions à 
l'instant initial { — 0. Pour formuler de telles 
conditions, nous développerons la fonction (1) 
suivant les fonctions propres |m%) de l’hamil- 
tonien /70, soit 


PO (1) = E aŸ(1)|m®). (27) 
01 


Considérons, parmi les N spins du liquide, tous 
ceux pour lesquels 


Em < Am(0) < Dm + dm 
bp < a»p(0) <a bp + db 
ta res at he rene AD Ua (28) 
by < agl0) < by + db 
RO D up NN SD SR She 2 Sont 


des nombres donnés. Nous dirons que tous ces 
spinsappartiennentàlaclasseC(b,, b,, ..., b,, ...). 
Des conditions initiales suffisantes pour qu'il y ait 
relaxation sont alors les suivantes : 


a) Tout spin appartient à une classe € et à une 
seule, 


nombre 
. de spins très 


b) Chaque classe © renferme un 
MC) Jb,(C) 90,(C) ... db4(C) 
grand, par exemple 


AC) SbmlC) 369(C) .… SbglE) ... > 108 (29) 


c) Les différents intervalles 
Sbu(C), DEC), - -…; BBC), 


peuvent être choisis suffisamment petits pour que 
l’on puisse prendre dans les calculs 


; 5b 
a (0) & bm + nn 


56 2 5 
aP(0) & bo +7, ag (0) © ba + ce 


(30) 
pour tous les spins 7 qui appartiennent à la même 
ASE A OR RS LPO À 


moyennes (18) que le système initialement considéré des N 
gpins, il lui est identique pour N — 0. Si lon fait 
abstraction de l'étiquetage des différents spins individuels 


(c’est-à-dire de leur numéro d'ordre 1, 2,... 7,... AN, 
N +1,... 2N) qui n'intervient pas dans la relaxation on 
aura d0(t) = 8(— 1) pour le système des 2N spins. Celui-ci 


se relaxe pourtant ! 


12 
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Pour que ces trois conditions puissent être satis- 
faites, il faut que le nombre total 


= E a(e ) SbmlC) 36D(C) ... SbglE) : 
de spins soit très grand. Dans (31) la sommation 
porte sur toutes les classes €. La signification 
physique des conditions initiales a), b) et c) est 
la suivante : à chaque spin j, arbitrairement choisi 
et d'état initial W®(0) correspondent un très 
grand nombre de spins #, {, ... dont les états 
iritiaux F®(0), WW(0), ... sont très «voisins» de 
(0) ; la répartition des états initiaux est donc 
supposée très «dense». Il semble «intuitif» qu’une 
telle condition est toujours réalisée dès que le 
nombre total NV de spins est suffisamment grand. 
Pourterminer, notons qu’il est possible qu’il existe 
des conditions initiales autres que celles que nous 
venons de mentionner, et qui conduisent également 
à l’équilibre thermique, mais nous n’avons effectué 
aucune recherche de telles conditions. 


7. La matrice de relaxation 6(#). — La valeur 
moyenne, à l'instant {, d’une grandeur physique G 
relative au 7eme spin est égale à 


(FO (GO FO (D) = 2 a *(r) aŸ (1) (n@ | GA mt) 
nm 

de sorte que la valeur moyenne totale de la gran- 

deur G, calculée pour tout le système de A spins, 

est égale à 

Ÿ te t) a) (+) 


< Gt) > w| (n|Gm) 


> Imist (t)|n) (n| Glm) = À trace de [o(t) G]. 
(32) 


pee 
=NS 


L'étude de la relaxation revient par conséquent à 
l’étude de l’évolution de la matrice de rela- 
xation o({) définie par 


ntotèin) = % EE a*(1 dl = ax) ant) (83) 


j=1 
et qu’il ne faut pas confondre avec la matrice de 
densité. La moyenne dans (33) porte en effet sur 
les N spins nucléaires du liquide, tandis que pour la 
matrice de densité, la moyenne porterait sur un 
ensemble de Gibbs de systèmes de Æ spins. Nous 
reviendrons sur cette question au paragraphe 14. 
Nous appellerons relaxation énergétique celle qui 
est régie par les éléments diagonaux de o(t) : C’est la 
relaxation de grandeurs G qui sont diagonales dans 
la représentation m, n, .. . Telle est par 
exemple la relaxation ‘de l'énergie ‘totale E(t) du 
système des W spins : 
ED =N> 
L/02 


(mlo(t)|m) E (34) 


où 


| 


En = (MmÔ|H®|m) = (m|H|m). (35) 
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Nous appellerons relaxation transversale [6] celle . 
qui est régie par les éléments non diagonaux de o(t): 
C’est la relaxation de grandeurs G qui ne sont 24 
diagonales dans la représentation m, n, ...,p, ..…. 
comme par exemple la composante transversale My x 
de l’aimantation (la direction Ox étant perpen- 
diculaire à la direction Oz du champ magnétique 


— 
constant H,). La relaxation transversale ne met pas 
en Jeu d'échanges d’énergie spins-réseau. 

Nous devons former maintenant l'équation 
d'évolution de la matrice o(t). Introduisons à cet 
effet la matrice o(t) définie par : 


ErQ ot] QD) = an * lt) am (1) (36) 
et qui obéit à ” 
in _. — [8e (1) oi (1) — of (+) JEU) (4)]. (37) 
L'intégrale o°(t) de (37) peut s’écrire 
oû(t) = TD(t) (0) TH *(1) (38) 


où T(t) est l'opérateur d'évolution relatif au 7eme 
spin et est donné par : 
TO (:) = , deu, (39) 

D’après la condition initiale c) du paragraphe 6, 

tous les spins 7 qui appartiennent à une même 

classe € (nous les noterons par j « C) ont initia- 

lement la même matrice o(0) = «t©(0) de sorte 

que 

> (md | (t)n@) 

ieC. , 

= E 2 (MIT (Hp) bare 
JeC »,a 


In@) (p@]0t? (0)[g@) 


= (plot@(0)[g) E (OT (t)|p@) (a IT #(e )In@) (40) 
p,a jeC 


en vertu de (38). De plus, on aura : 
1 


AC) SbmnlE) dbplE) -.: 06 (CE): 
Z (MD ITA (1)[pD) (gAITA* (IRD) = (m[T (1) |) (gl 7 #(0)1r) 
eC 
à (4) 
avec 
(m|T(t)lp) (gl T#(0)1») 
LS ( (mOITO(t)|pA) (g#IT#(HIn@). (42) | 
AY ÿ=1 


En effet, le membre de gauche de (41) représente 
une moyenne prise sur la classe €. Mais comme en 
vertu de l'hypothèse b) du paragraphe 6, toute 


classe € renferme un très grand nombre 
(EC) Sb,(C 2) 5B,(C) . DC) :.. de spinsnt 
cléaires, qu’on de supposer répartis un peu 


partout dans le liquide, la moyenne sur € et la 
moyenne sur les V spins de tout 1 liquide sont 


No3 


N°3 


sensiblement égales (8). C’est Ià le sens de l’hypo- 
thèse b). L'égalité (41) serait strictement vérifiée 
au cas où toutes les classes renfermeraient un 


nombre infini de spins, ce qui nécessite N — co. 


à z (m@ | (1)|n) 


Sommons (40) sur toutes les classes ©, il vient : 
N{mlo(t)|n) 


— > {5 (Ce) 5 (C) Sb(C) .….. SblC) .…. (pl ©)(0)[g) | 

?.4 | C | 

| (m|T(t)]p) (g|T*(0)n) 
RUESN 2 Cf (ip) WIE* 0) (ple(0)le) (43) 


\ 


Telle est l'équation d’évolution de la matrice o(t), 
compte tenu des conditions initiales du para- 
| graphe 6. 


8. Calcul de perturbation au deuxième ordre en V 
de la matrice o(t). — Pour calculer o(t) au moyen 
_ de (43), nous développerons l’opérateur d’évo- 

lution 7(t) en puissances de la perturbation V(+) : 


î (1 LA 2 
T{t) = a, requt € To (+) 


en (44) 
où 7(t), d'ordre n en V, est donné par 

? _ jn ft t te 
T(t) = (=) L dt, 1. pans £ dt, (45) 


eiE—hn V(ty) e—An-tn—1V8 Vins) 


et 


Vu) ei 


TO(+) = e—iAU, (46) 


Dans son domaine de convergence, la série infinie 
qui figure à droite dans (44) représentera exac- 
tement l’opérateur d'évolution T(t). La signi- 
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fication physique du développement (44) est très 
claire. La terme (45) d'ordre » en V de ce dévelop- 
pement représente des processus au cours desquels 
la perturbation V provoque n transitions succes- 
sives aux temps f1, lo, ..., l, alors que pendant les 
intervalles (0, #1), (f, ta), «.., (6, t) évolution du 
système est déterminée par l’hamiltonien À non 
perturbé. 

Lorsqu'on porte (44) dans (43), on obtiendra un 
développement de la matrice ot) de relaxation en 
série de puissances de V. Nous calculerons dans ce 
paragraphe tous les termes jusqu’au deuxième ordre 
inclus de ce développement ({). Au second ordre 
en V,on aura 


(mlo(t}|n) = Z ZX 


(RIT (s)1p) (17% (11) (plo(0)|o). 
D,Q k+I<2 


(47) 


Les termes d’ordre zéro en V dans (47) s’écrivent 


moltn)® = © (raje-@|p) (gl iÆain) (plo(0)|a 
— eionmt(mlo(0)|n) à 
avec 
D ner En, (49) 


ñ 


Les termes du premier ordre en V dans (47) sont 
nuls en vertu de (19). Les termes du deuxième 
ordre en V dans (47) donnent 


(mlo(t)1r)® — es ÈË (m[T@ (+)[p) (gl D X(1)1n) (plo(0)|g). 


(50) 
On a, en vertu de (44) et (45) : 


E UT Ep) IT * (1) (lea 2 EF “dés [des e=iBntt Url) Ur PP) 
?,a DM 0 0 


ei Etat) à e—iEptlh e+iEntlh (plo(0)|n). (51) 


Pour calculer (51), prenons comme variables d'intégration t, — {, = u et t,. Il vient, en vertu de (18) 


> (TE lp) (TO % ein) (plo(0)1e) 


ñ 


rh ARE à à 
26 cionnt fl de, eiompts le du (mn Plu)l) (ral V(0)|p) eiopmt (plo(0)|n). 


(5) En fait, il y a ici une hypothèse supplémentaire, de 


… même nature que celle du chaos moléculaire, à savoir que 
- les fonctions d’onde des différents spins à l’instant initial 


ne dépendent pas des opérateurs TG), C’est le « chaos des 


- spins » dans le liquide à l’instant initial. Certains considèrent 


\ 


une telle hypothèse comme étant de nature statistique. En 
réalité, cette hypothèse n’introduit aucun postulat statis- 
tique dans la théorie, elle n’est pas statistique au sens que 
nous donnons à ce terme au paragraphe 14. 


dis 1 ct 


(52) 


(4) Y. Ayant, dans sa théorie de la forme des raies 
de résonance nucléaire [8] a évalué les expressions 
(m|T{(t)|p) (g/T*(t)in) jusqu’au deuxième ordre en VW. La 
raison pour laquelle ces expressions apparaissent dans la 
détermination du spectre d’absorption est expliquée au 
paragraphe 15. Nous ne suivons pas ici l'exposé d’Ayant 
parce que sa méthode de calcul est très particulière et ne 
saurait servir d'introduction à un calcul général de (43) à 
tous les ordres en V. Notons aussi que 'es moyennes (indi- 
quées en surlignant les expression en question) n’ont pas le 
même sens chez Ayant qu'ici. 
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Prenons t >. Nous pourrons alors rem- 
placer {, par co comme limite supérieure de l’inté- 
grale en du, puisqu’en vertu de (21), la quantité 
sous le signe somme (en du) est négligeable dès 


(RIT ®()]p) (IT O* (In) (pls(0)|9) 


—2 


Om 


en tenant compte du fait que les niveaux de À ne 
sont pas dégénérés, et en écrivant séparément les 
termes (53 b) pour lesquels p = m, les termes (53 a) 
étant ceux pour lesquels p Æ m. Pour t > o;}, 


D: (mL (t)|p) (gl2 2% {t)[n) (plo(0)|g) = E (nILE(t)|g)* (pIT O*#(e)|m)* (glo(0)|p)* 


D,4 


Mai 


puisque (54) se déduit de (53) en interchangeant 
les rôles de m et x et en prenant l’imaginaire 
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que u > + et que l’on a f, > +, pour la grande | 
majorité des valeurs de f, comprises entre 0 et t. 
Il vient ainsi : 


(7) SE tetomot 1) au Ta rai eicom (ple(0)n) 


+ > le ) eiSnmt re du Tr PUu)mes) (ral (0) fm) eiommau (alo(0)|n) 


ce ÉL Ai [© du (n]P{u)[m)* [re] V(0)in)* e—iommau (n|o(0)[me)* 


(53 a) | 


(53 b) 


les termes (53 a) serort négligeables devant les w 
termes (53 b) proportionnels au temps £. 
De même, on a : 


(54) 


2 Em]TO(s)1p) (IT %(&)1») (plo(0)|g) 


HS () () f° dti dés etEmthlä (m]V(t1)|p) (r]V(63)19)* 


D,Q 


Prenons comme variables d'intégration t, — t, = uet f,. Il vient, en vertu de (18) : 


E (m]T@(1)[p) (gTO*(1)1n) (plo(0)|g) 


2.4 


D,q 


Prenons { > +. Nous pourrons-alors remplacer 
les limites — a(f:) et + «(f,) de l’intégrale en du 


É PRIT (2)]p) (IT O#(0)1n) (ple(0)lg) 


qe elomp—@nqlt — 1 
En CON 
ñ mp — na] 


TAN ?7 


Om ON 


£ 53 eionnt t Je LT du Pr) + n — miV(0)In) eiompu (plo(0)|p + n — m) 


où il a été fait usage de ce que V(t) est un opérateur 
hermitien. (57 b) renferme tous les termes, pro- 
venant de (56) et pour lesquels o»» = on. Nous 
posons 


En = — ÿà Hom (58) 


= + eionnt 3. | dé, eitomp + Gr 1h CONIOUONE CT (ple(0)la). 


[2 du Ont Ip) Cal PO) emo (ple(0)) 


conjuguée, (54) est, tout comme (53 b), valable À 


pourt> tett >o,. Enfin, on a: 


eEphlh eiEn(t—ts)/ñ ei Egtalh | plo(0)|g). (55) 


(56). 


— (ts) 


respectivement par — co et + co, et il vient 
ainsi : 


de sorte que om — &ng entraine 
M—p=n—g 


(59) 
d’où (57 b). Pour t > [om — onl; 1e88 


# D rines (57 a) seront négligeables devant les 
… termes (57 6) proportionnels au temps f. 
3 Faisons alors la somme des termes (48), (53 b), 


Emlottiln) = Emlois)|n}® + (rmfaitiin)® 
— eiSnmt (m\c(0)|n) 


5 RE co - 
—_ comm du (m| 
= | 


LL 


1 
+ jaeisunt eZ [UT du 


| 9. Les équations de Boltzmann auxquelles obéit 

la matrice de relaxation. — Supposons que la 
« perturbation V soit suffisamment petite pour 
que (60) représente une valeur approchée de l'élé- 
_ ment de matrice (m|o(t)\n) dans l'intervalle 


LEE (61) 
avec 
PR 1,1 es. - (62) 
Introduisons alors la matrice 8(t) définie par 
(m|g(t)in) = —_— (mlo(t}[n). (63) 


- La matrice £(f) représente la matrice de rela- 
_ xation dans le trièdre de coordonnées tournant à la 
… fréquence angulaire de Larmor 6 = ;A,| autour 
de Oz (direction du champ magnétique cons- 
. tant H,). Pour des temps ft satisfaisant à (61), on 
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1 RP 
RE eienmi t sf du (n|V{u)|m;) (m,|V(0){n) ist (]0(0)[n) 
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(54) et (57 b). Nous obtenons ainsi, en vertu de (47) 
et (50) : 


Mm|V(0)|m) eiemmt (m|lo(0)|n) 


(60) 


(| Vu)|p} (p + n — m|V{(0)\n) eiempt (plo(0)|p + n — m). 


pendante du temps. Supposons alors que les condi- 
tions initiales a), b) et c) du paragraphe 6, valables 
à linstant t — 0, soient encore valables aux 
instants é, 2f, ...,rt, .... La matrice { — Af qui 


* transforme (0) en &(f) transformera alors 
. B(Er — 116) en B(ri) : | 
; Bret) = (1 — Ai) 8([r —1j2 (65) 
d’où l’on tire 
S(rt} = (1 — Arjr 8(0). (66) 
Posons ré — u, il vient 
Blu) — (1 LE A} 8(0) & e—-Au 8(0) (67) 


pour # > +, u > ©, lorsque r — co. Déri- 
vons (67) par rapport à x, et remplaçons x par t, 
il vient 


D À at) = — Age) (68 
> 8 =H—AQ 80) (64) DU Pi (68) 
“où À est la matrice explicitée en (60) et indé- ou encore, en explicitant la matrice A (60) : 
3 {mB(c){n) = — er = du (m|V{u)lm:) (mel V(0){m) em (m|8(t)Îr) 
dé CE 0 1 z 4 ET - Ca 4 
— 2 È E du Pam) Oral TON) et (m6) \ (69) 
TE > 1 du (m|V{u)|p} (p + n — m|V(0){n) eismpu (p|6(t)|p + n —m) 
TE 


…— (68) et (69) constituent les équations de Boltzmann 

ur la matrice de relaxation 8(t). Les équa- 
Fons (67), (68) et (69) résolvent complètement le 
_ problème de la relaxation énergétique et trans- 
sale. On pourra les comparer avec les équations 


a ti) = 55 SR 


obtenues par Wangsness et Bloch [5] pour leur 
- «matrice de distribution ». 
Pour étudier la relaxation énergétique, faisons 
n — m dans (69), il vient : 


du lu) Oral PTON) eme (la () me) 


(70) 


D . 32 sf du Unie) (PP (ON) eicmse (pla(s)(p). 
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Ces équations (70) sont bien du type (3) avec 
P,(t) = (mlo(t)|m) et 
1 


Wap == 


= (71) 


+ du (m] V{u)|p) (pl (0)|m) eicmpu. 
—CO 

Rappelons que (69) et (70) sont valables à condi- 
tion que les inégalités (61) et (62) soient ous 
tibles entre elles, il faut pour cela que +,, w et V 
soient suffisamment petits. On voit sur (70) que la 
relaxation énergétique est en ce cas indépendante 
de la relaxation transversale régie par les éléments 
non diagonaux de {(t). 

Notons enfin que des cas particuliers de (71) ont 
été obtenus par Bloembergen [7] et par Ayant [8], 
qui ont utilisé pour les établir la méthode de varia- 
tion des constantes de Dirac, c’est-à-dire une 


méthode ad hoc. 


10. Calcul de perturbation au i‘* ordre en V 
de la matrice o(t). — Dans les paragraphes pré- 
cédents, nous avons déduit les équations (67), (68), 
(69) et (70) à partir de l’équation de Schrô- 
dinger (22) en faisant l’hypothèse que les condi- 
tions a), b) et c) du paragraphe 6 étaient vérifiées 
non seulement à l'instant initial £ — 0, mais aussi 
à toute une série d’instants successifs très nom- 
breux 


LISA RRET De 


Il en résulte que notre démonstration de (69) est 
encore incomplète. En effet, si à l’instant initial 
t — 0, la fonction d’onde ®(0) du système des NV 
spins est donnée, et vérifie les conditions a), b) et c) 
du paragraphe 6, toute l’évolution ultérieure de ®(t) 
est déterminée par l’équation (22). Pour rendre 
corrects les développements du paragraphe pré- 
cédent, il nous faudrait donc prouver que les condi- 
tions a), b)et c) du paragraphe 6 sont vérifiées à tous 
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les instants (72) si elles le sont à l'instant & = 0. 
Mais il paraît difficile de prouver cette hypothèse, " 
bien qu’il semble qu’elle soit toujours vraie. C’est 
pourquoi nous ne poursuivrons pas dans cette voie, « 
mais étendrons le calcul de perturbation de o{t)« 
à tous les ordres en V. 

L'expression (60) représente le développement de 
la matrice o(t) en série de puissances de V jusqu’au 
second ordre en V inclus. De ce fait, cette expres- 
sion n’est valable que pour des temps f inférieurs à 
un certain temps {, qui apparaît dans (61). Ce 
temps £, est généraiement faible devant les temps « 
de relaxation T,. De même, les équations de « 
Boltzmann (69) et (70), qui ont été déduites au M 
paragraphe précédent de l’expression (60), n’ont M 
ainsi été établies que pour des temps é compris 


entre t, et {.. Afin de démontrer que ces équations « 
restent valables pour des temps { > #,, il nous " 


faudra calculer une expression exacte de o(t), c’est- » 
à-dire évaluer tous les termes du développement 
de o(f) en puissances de V et non plus seulement 
les termes jusqu’au deuxième ordre en V. Nous 
allons voir qu’un tel calcul, encore qu’assez com- 
plexe, est considérablement simplifié par le fait 
que la perturbation V(t) est une fonction aléatoire 
stationnaire du temps au lieu d’être une fonction 
arbitraire. Cette particularité permettra de pousser 
l'analyse exacte de la relaxation beaucoup plus loin 
qu’on ne le faisait jusqu’à présent et d'obtenir une 
image physique simple pour les processus de rela- 
xation d’ordre supérieur à deux en V. 

D’après (43), le terme d’ordre L en V dans le 
développement de (m]o(t)|\n) est donné par 


Pate = SE DT Ep ITR (ple(0)lg). 
(734 


Les T®(f) étant donnés par (45), on aura 


lb t U + \ 
dt [au f dt xs 2 À 
[ 1 0. L 5 1-1 5 k+1 


) (m|V(tx)|mx-1) 
Mk 1 Mk — 9... 
M1 M — Done 


(mal Vtx=21) |M%—0) : - (mal V(t 


eiEm(t-o/i—iEmy_ Ut à hi... 


Pour calculer l'intégrale (74), il nous faut étudier 
la moyenne sur le produit d'éléments de matrice 
de V qui y apparaît et que nous écrirons symbo- 
liquement 


V{tx) Vis) +. V{ts) Vu) Vus) <<: Vitres) 
selon une abréviation déjà utilisée au paragraphe 4. 
Nous effectuerons le calcul pour des temps { véri- 
fiant 


(75) 


1 +. (76) 


lp) Vu) ru) 


—iBptli eiEn(t-tli+ Em _ tt ait …. 


*(my1|Vts)frue)* .. (me sal V(tx+3)|g)* 


. + iEgtg+a/h, (74) l 


Désignons par 
V(t) V(ts) Vity) -» 


une fonction égale à 


(77) 


Ver) V{ts) V{tu) (78) 


lorsque toutes les variables £,, £,, 4, ... qui y 
figurent sont comprises dans un intervalle d’éten- 
due = et égale à zéro lorsque les variables f,, 4,1... 

sont éparpillées dans un intervalle d’étendue supé- 
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rieure à 7. + désigne un intervalle de l’ordre de 
grandeur du temps de corrélation +. Nous appel- 
lerons une telle fonction (77) un complexe. On aura 
alors 


V{tx) V(tx_1) .. PE] Vu) Viu-1) +. V(ix+3) 
Cie Mie) 


la somme Z’ portant sur toutes les manières pos- 
sibles d'associer les V4), Vi), ..., V(t+1) qui 
apparaissent dans (75) en complexes. Ainsi, dans le 
membre de droite de (79), nous avons explicité 
trois complexes : le complexe (4) V(4,) . .. %(1,), le 
complexe Ÿ(4) V(t,) ... et le complexe (4) ...…. 
Toutes les variables #, #_;, PP EN 
qui figurent à gauche dans (79) figurent également 
AidrObren des ln... Le. inais leur 
ordre peut avoir changé. Le terme à gauche 
dans (79) est de degré / en V, tous les termes de la 
somme Z’ sont également de degré ! en Ÿ. Enfin, 
un intervalle toujours supérieur à + doit séparer les 
variables de temps de deux complexes quelconques 
distincts appartenant à un même terme de la 
somme 2’ ; dans le cas contraire, le terme corres- 
pondant doit être pris égal à zéro ; cette restriction 
dans la sommation (79) est indiquée par le prime 
qui figure à côté du Z. Il résulte de cette restriction 
que pour une valeur donnée quelconque des va- 
PO E ,  U L...., 41, seul Un 
terme de la somme ZX’ est différent de zéro, tous les 
autres termes de cette somme sont nuls. L’éga- 
lité (79) est sensiblement vérifiée en vertu du fait 
que la fonction (75) obéit à (17) et compte tenu 
de (19). On a naturellement 


Lo h UT 


(80) 
Un complexe ; 
Vi) Vo) Viti) +. Vs) V(tr) 


(81) 
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de degré 3 en Ÿ ne dépend en réalité que des d — 1 
variables 
(82) 


t — Lys To — is... te — ty 


ceci à cause de la stationnarité du thermostat, 
exprimée par l'égalité (18). Portons alors la décom- 
position (79) de l’expression (75) dans (74). Seuls 
interviendront les termes de la somme X”(79) pour 
lesquels 


DR A le QUE CE) 
D CNE 0 | 


et nous indiquerons cette restriction par le signe X”. 
L'intégrale (74) devient alors une somme d’inté- 
grales portant sur des produits de complexes et des 
exponentielles. Pour calculer l’intégrale portant sur 
le produit de complexes explicité à droite de (79), 
nous prenons comme variables d'intégration 


(83) 


NU ge Di Un el} Ur, 


ee eo pe 
et 


lis Tps es lp, A... 


Les limites d'intégration relatives aux intégrales 
portant sur les différences (84) seront soit zéro et 
l'in fini (pour les différences t; — t, pour lesquelles t, 
et t, figurent tous deux sur une même ligne de (83)),. 
Soit — co et + co (pour les différences {, — #, pour 
lesquelles t, et t, figurent respectivement sur la 
ligne supérieure et sur la ligne inférieure de (83)). 
En effet, par définition même des complexes (77), 
ceux-ci sont nuls dès que l’une des différences (84) 
dont ils dépendent devient supérieure à +. Comme 
nous avons pris { > r, la plupart des valeurs prises 
par les différences (84), à l’intérieur du domaine 
d'intégration sont en valeur absoluesupérieures à 7, 
etil n’y a par conséquent aucun inconvénient à 
étendre les limites d'intégration soit à 0 et + co 
soit à — co et + co suivant les cas. On obtient 
ainsi 


; 0 t t : \ 
eionnt [+144 f ASE (tique Fntn+ + Fat) 
0 0 0 


orne ses screen ess ressens ses ess 


nt eltials etes ee ip alele ele ann lohela ee as sets 5 » ee ve 4 0 jee à 0 


DURE... Pi. Fu... Fe... 
sont des constantes réelles, indépendantes du temps 
(ce sont des combinaisons linéaires des énergies Æ,, 
Eng) +: + Es En, : +. Ej) que nous calculerons à la 


fin de ce paragraphe, Parmi les intégrales portant 


+ ir dusr (my| (45) |my—1) (mal V(to)|Mg-x) » - « 


im IV (t)[tn, 1) ei Frotigt + Ferusr) 


sur t, tn, ...,t,, -.. figurent soit l’une, soit l’autre 
des deux intégrales 
u 
dt 
0 


! 
[ dty, 
0 
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soit même les deux. (Ces deux intégrales figuraient 
en (74)). En effet, la variable #, apparaît une fois 
dans chaque terme des sommes Z” (79) et Z” (86). 
Nous rangerons f; parmi les variables du type (85), 
ce qu’il est toujours possible de faire puisque parmi 
les variables t,, t,, ...,t, dont dépend un complexe, 
on peut choisir arbitrairement celle qui sera placée 
parmi les variables (85). De même, on pourra 
toujours ranger { parmi les variables (85), sauf 


si #, et t font partie du même complexe V(#) (4)... 
auquel cas seule l’une des deux variables f; et b 
pourra figurer en (35). 

On vcit alors que, puisque les intégrales figurant 
dans (86) et portant sur les variables u,,, 
pis ++ Us = + SONt des constantes indépendantes 
du temps #, l'expression (86) est un polynôme en t 
dont les coefficients sont des exponentielles en # à 
argument imaginaire, c’est-à-dire un polynôme en t 
dont les coefficients sont des fonctions de t qui 
restent bornées quel que soit £. 

Nous allons déterminer maintenant le degré en t 
de ce polynôme (86). Considérons le terme de la 
somme X” explicité à droite de (86). Le monôme de 
plus haut degré en { qui figure dans ce terme 
s’obtient en y faisant 


Fe Fe sers 


. =-0; 


(87) 


Le degré en t de ce monôme est donc égal au 
nombre de complexes qui apparaissent dans le 
terme explicité (86) de la somme XŸ”, puisque chaque 
complexe fournit une intégrale respectivement 
sur t}, lp, ..., ls -.. Le nombre maximum de 
complexes qui peuvent apparaître dans un même 
terme de la somme 2” (86) est égal au degré / en? 
de (86), mais comme ®(t) — V(t) — 0, le nombre 
maximum de complexes différents de zéro est 
seulement égal à //2 lorsque / est pair, et à ({— 1) /2 
lorsque l'est impair. Pour ! pair, le monôme de plus 
haut degré en t{ dans (86) contiendra //2 complexes 
d'ordre deux en ?, du type 


Fe du (ml Vus) |my1) (mal V(0)[my-1) eiPrauto (88) 


Pour l'impair, le monôme de plus haut degré en t 
dans (86) contiendra ({ — 3)/2 complexes d’ordre 
deux en Ÿ du type (88), et un complexe d’ordre trois 
en Ÿ. ‘ 

Le développement de la matrice de rela- 
xation f({) en puissances de Viest donc de la forme 


B(t) — B(0) + $ (= vien (1) | LS CE 


FRERE ) | 
1=1 Ly=0 1=1 Ly=0 


(89) 


où les Va) t) sont Le fonctions bornées de { d'ordre x 
en V (ce sont des combinaisons linéaires d’expo- 
nentielles en { à exposant imaginaire). 

Nous calculerons (89) pour des temps { sufi- 
samment grands afin de n’avoir à retenir dans 
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chacun des polynômes figurant entre crochets que 
son monôme de degré le plus élevé en t. Le terme 
avec « — 2] pair fournira donc la contribution 


v{2n 41 (90) 


tandis que le terme avec à = À +1 impair four- 
nira la contribution 


Vert à (91) 
négligeable devant (90) lorsque V est suffisamment 
petit. Supposons que V®)f soit une quantité finie. 
Il en résulte que seuls les termes avec « = 2} pair 
fourniront une contribution sensible à (89). Pour # 

suffisamment grand, on aura 
se 1 
Bite B(0} Een 7m, (92) 

1=1 

Pour interpréter physiquement les résultats 
obtenus, rappelons-nous que les instants t, et f, qui 
apparaissent dans un complexe d'ordre deux en V 
du type (88) sont séparés par un intervalle de temps 
inférieur à T, tandis que les instants qui appa- 


- raissent dans deux complexes (88) distincts sont 


séparés par des intervalles supérieurs à +. On peut 
donc dire que seuls contribuent à la relaxation des 
processus dans lesquels la perturbation V agit à des 
instants #, , ..:, . Lists VU restent groupés 
deux par deux: (#,;, ,), ...,(4,t,), .... Les instants 
dans un même groupe, tels que {; et 4 par exemple, 
sont très voisins, l'intervalle qui les sépare étant 
inférieur au temps de corrélation +. Les instants 
dans deux groupes différents, tels que t, et t, par 
exemple, sont séparés par un intervalle de temps 
appréciable, supérieur au temps de corrélation .. 
Partant de (6), Van Hove obtient des résultats 
différents dans la forme, mais identiques quant au 
fond, puisqu'il trouve également le groupement 
deux par deux des instants auxquels agit la per- 
turbation. Cette concordance des résultats sur ce 
point indique qu’il doit y avoir une relation entre 
la condition (6) de Van Hove et les conditions de : 
l’auteur, relation qui sera élucidée dans un travail 
ultérieur. 

Nous devons maintenant expliciter les condi- 
tions (87) et pour cela, calculer F;, F,, ... Les 
exponentielles qui se rapportent au complexe 
(nm V(t)im,_.) (m,|V(t)1m,_,) et qui renferment les 
variables t, et {, sont, d’après (86) et (74) 


eiFfott-to) eiFitf eitEm-Emy_)tf/ñ eitEm y Em, 


D'talf 


d’où l’on tire 


FF; = On y et Er — 3 Pny- 1 Eng Æ En = 

(3j 
La condition (87) fournit alors 

4 My = My + Magny — My. (94) 


No3 


_No3 


De même, les exponentielles qui se rapportent au 
complexe 


(ms|V{ts)[ms_3) (rexl (tr) [mer—1)* 


= (ms|V(t)|(ms-1) (ml V(t) rer) 


et qui renferment les variables f, et t{, sont d’après 
_ (86) et (74) 


etFsr(ts—tr) eiFsts er eilEms—Ems_) ts/h e Em, Em,_,) tr 
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d’où l’on tire ; 
Pi = Onpm.… et Pi ne : Ent En, 
: (95) 
La condition (86) fournit alors 
My = Mai EE My — Ms. (96) 


Finalement, compte tenu de tous les résultats 
obtenus dans ce paragraphe, nous pourrons 
écrire (86) sous la forme 


\ À /1\ 21% 
(nTT® (1) [p) RTS = (—) (;) 2" S eionnt [ LE Fr jé Ad 
0 


ñ ñ 


co ; 
à ! duyg (ml Vue) ms) (my + My — my] V(0)[mg) ro 


(97) 


+ co s : 
# L  duss (ms] V(usr)|ms-1) (Ms + my — ms|V(0)1m;) CORSA 
CO 


CRE CCC ET 


La somme ZX” porte sur toutes les manières 
possibles d’associer les V(#), V(#_:), ..., 


_Vitx) ... V(&x+,) qui apparaissent dans (74) en 


complexes d’ordre deux, compte tenu des condi- 


- tions (83). On a naturellement remplacé / (figurant 


dans (74)) par 21. Seuls apparaissent des complexes 


_ dutype 


a] Vlu)|6) (b + e — a] V(0)jc) eioapu 


figurant déjà dans l'expression (60) de o{t), expres- 
sion exacte jusqu’au second ordre en V inclus. 
Cette circonstance laisse présager le fait que les 
équations de Boltzmann (69) obtenues au para- 


(98) 


graphe 9 à l’aide d’un calcul de perturbation au 


- deuxième ordre en , sont en réalité valables à tous 


\ 


L Leiomt > PTE Gp) MIT PIE (plo(0) [gi] 
D,4 


RS 


les ordres en V. C’est ce que nous allons montrer 
maintenant. 


11. Démonstration générale des équations de 
Boltzmann. — Nous sommes maintenant en mesure 
de passer à la démonstration des équations de 
Boltzmann (69) en n’utilisant les conditions a), 
b) et c) du paragraphe 6 qu’à l'instant initial { = 0. 
La méthode que nous allons suivre consistera à 
établir d’abord une relation de récurrence 
entre B(t)2 et B(t)2-%, ou plus exactement entre 
dg(#)2? Jdtet B2-2(4). Rappelons que B(t)®? contient 
dans le développement de (t) en puissances de V, 
tous les termes d’ordre 21 en V. On a : 


nous calculerons le terme (99 b). Pour évaluer 
l'intégrale (99 a), nous devons d’après le para- 
graphe précédent, décomposer le produit des V{(t), 
V(&-), - Vi), Vu), Vars), + +, V(te+3) en 
complexes d’ordre deux en V. Portons notre atten- 
tion sur les complexes dans lesquels figure (+). 
Parmi ceux-ci, seuls peuvent apparaître dans (99 a) 
les complexes NE OCR 


PO Va. 7, (oo) 


— jE Jak t la i t 4 
ES (3) (ple(o)la) f du f pe (US d'a L au / deu [dun es fran | (99 a) 
(mleiAt/ñ 40) eriH(t-tx—1)/à 02) ER eTiHt/ñ|p) (nleiÆtau/à 40) e—iH(toy-tar1)/h ..., e—tAtk+r/i]g)* 
: NE Ji\ 21 t tk la t bla 4, [ +2 
+ E ( JE) wletoilo fan fran. … fan fau fra … fran La 
(mleiAtx/h V (tx) e—iHig=tgr e—tHt'ñ]p) (nJeiAt/à V{r) etant /n Vitry) ettx+1/|q)*. 
Nous allons d’abord calculer le terme (99 a), puis et F0 7). (101) 


En effet, supposons ab absurdo qu’il puisse appa- 
raitre un complexe V(t) V(#;_,) par exemple. Ce 
dernier complexe n’est différent de zéro qué pour 

t— ty LT. (102) 

Mais lorsque l'inégalité (102) est satisfaite, il 

résulte de (83) que l’on aura a fortiori 


tn <= tp < Tr (103) 
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Mais (103) est en contradiction avec le grou- 
pement des temps L'REN- llen cr Re 
par deux tel qu'il a été expliqué au paragraphe 
précédent, puisque le temps { du groupe (f, #3) 
ne serait pas séparé du temps be du groupe 
(te t) par un intervalle supérieur à +. On peut 
encore dire que lorsque l’inégalité (102) est vérifiée, 
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il apparaît le complexe - | 
ae V1 à V2) V(t3) 


que la contribution à la relaxation est res bA 
Notre affirmation est donc établie. Le terme (99 a 
est donc égal à 


t : S = = z 
» > , (plel0)lel dy (m|V(t) ) my — :1) Dar eee eiPmt/ñ( (Mx— LA tx) |m) e—Emtx—1/# 


HE dt, 1 dis Je dre 


DA My 


Bee 21—X ftx—1 # 
= = () À de ie EU 


(04 a) M 


k+2 
Si dép +1 


(mletEtk=21% Vi) ego sf, 


. e7#h/ñ|p) (n| eAtauñ Vi) …. 


ee iHix+1/ñ|g)* 


+ Z ZX (plo(0) 


D, M 


= ARJNA-R ft Dpt 
X (=) () Î dns f di …. 


t x : 2 
la) [ digg ePmtlh tal V(t}{mx) (nl Via) mr + n — m)* e Em tr e+ÜEmy 5 tn—mEntauh 
0 


t2 121 b21—1 
L dt, Faut f dires re 
0 0 


X (my let —vr Vis) AE 


Pour établir (104), nous avons tenu compte de ce 
que seuls pouvaient apparaître des éléments de 
matrice du type (98). De plus, nous n’avons expli- 
cité, dans la décomposition en complexes d’ordre 
deux en V de 


V{t) Vins) +. Vu) V{iy) AU ... V(x+) 


(qui apparaît en (99 a)) que les seuls com- 
plexes (100) — qui apparaît en (104 a) — et (101) — 
qui apparaît en (104 b) —. L'intégrale 


É 
1. DR on 
0 


qui figure dans (104 a) peut être remplacée par 


t 
. SEM E 


puisque l’expression (104 a) est sensiblement nulle 


(105) 


(106) 


e tHtu/A|p) (Mg + n — m| etAtey—1/à Vi) et 2 tour 


RUE 
(104 8)/ 


le+ 2 À 
1: dip+s ‘ 
0 j 


dès que { — ty … > +. Il en résulte que la substi: À 
tution de (106) à (105) ne modifiera pratiquement * 
pas (104 a) puisque nous effectuons le calcul pour 4 


des temps { tels que 
1D re (107) 


L'expression (104 a) est donc égale à 


eiômnt È 2e dul (m]V{u)|a) (a[ V(0)[m) eiomau h 
2 à | 20 

5 UT (D 1p) (ITEE-# (119) * (plo(0)|g) (108) 
comme on le vérifie aisément. On a posé 


t— 1, = u et l’on a remplacé m;_, par a” De 


même, on pourraremplacer l’intégrale k feu diy=08 
qui figure dans (104 b) par be DURE: 3 


que l’expression (104 b) devient égale à 


. de sorte 


eiomnt © 22 f © du (m|V{u)|a) (a + n — m|V(0)]n) eiomau 
D, a À 


à) (+ à) TE D pl le + MIT 


où l’on a posé — 1, — uet où l’on a remplacé my, 
par a. 

Nous pouvons calculer de manière identique le 
terme (99 b). Dans la décomposition de 


Vi) Vs) + Vies) VO Vis) 2e Vlm+3) 


(qui figure en (99 b)) en complexes d’ordre deux 
en V, fixons notre attention sur les complexes qui 
contiennent V(+ (t). Ceux-ci ne sont qu’au nombre de 


(110) 


*{plo{0)lg) (109) 


deux, à savoir 


V{(t) V{tx) (111) 


et 


V{t) V{ar-1) (112) N 


les autres n’apparaissent pas. La démonstration est | 
identique à celle exposée pour le calcul de (99 a) et . 
ne sera pas réexposée ici. Le terme (99 b) est donc 
égal à ‘4 

. 
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t I Dr PES Le ; 
D © (ple(0)|g f° dé em El PG)fma) Qui P rs nm ent of nn 0 
D,Q me 1 


> f—i AY 16 + ‘ 
X + () ATTES Le dt 5 lens 1e du f° dt A modes a fr die (113 a) 


X (myiletÆ1/à Vty=) ea te2 he iBh/ñ]|p) (my + n — m] eiltar1/ù Vis) ... ete +1/i]g)* 


nt 


t Se ee SA Fos 
“+2 ZE (lo(ole f due Pati CPE ma) Ur Pise ent rétine all 
D,4 m 


nn” 


2-1 


A: ET ANCr Er t : c 
x | —) (à d'y 1 dia ef dt, ea SE ALTO UE (113 b) 
0 


X (re iltylñ Val e—tHulh|p) (n| etHtoy-olh Vi) e A (toto tot sh... etHtk+1/2)g)*, 


__ Parunraisonnement identique à celui qui permet de passer de (104 a) à (108) et de (104 b) à (109), 
on démontre que (113 a) est égal à 


com D 22 [7 du PalPRat je + n— mi P On ion 
?,q a À" J—0co 
É = 2) (0) (QT 9 (t)|p) (a + 0 — mITÉEE D (1)]g)* (plo(0)g). (114) 
On démontre de même que (113 b) est égal à Finalement, la quantité 
| ciémt © DE : du (n|V{u)|a) CAO eionau [oo ZE (m|T® (:))|p) (R| TER (1) 19) * (plo\0, a)| (116) 
p,q a À J—co dt D,q 


est égale à la somme des expressions (108), (109), 
(114) et (115). Si nous sommons cette égalité sur 
indice k de À — 0 à k — 21, il vient 


(2 (al a lp) GITE 219 (plo(0)|g. (115) 


Goom TOO MT (plo(o)i] 
=  Segope 
= — == f © du (miel) (al {0)fm) eiéman (m|B(0)|n) C1 (117 a) } (117) 
+ 2 [1 du QnPtpfa) da + nm TO re etc mat (alB(r)la + n—mj8-9 (117 à) 
— 2 [du PT talPTONr) ei nav (Be) 1) et. (117 0) | 


Le terme (117 a) est obtenu en sommantl’expres- tion (117) constitue la relation de récurrence 
sion (108) Sur l'indice k. Le terme (117 b)est obtenu annoncée entre B2(4) et 82-2(7). Sommons alors la 
en ajoutant (109) à (114), puis en sommant letotal relation (117) sur l'indice ! de ! — 0 à ! — oo, il 
sur l'indice 4. Enfin, lé terme (117 c)est obtenuen vient : 
sommant l'expression (115) sur l’indice k. La rela- 


d< 


LE, UnlBtonen = "À (mit 


= — 22 [du EnPTaa) CaiP TON eèma% (ral) (n) 


22 [7 du WT GPO e*one% (m]8(0n) (59) 


+52 07 du GPqufa) Ca + n — mIP CON) eème* (a|B(s)la + n — m) 
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12. Discussion des équations de Boltzmann. — 
En démontrant les équations (69) à partir de l’équa- 
tion de Schrôdinger, nous avons établi l’irréver- 
sibilité du comportement d’un système de A spins 
non seulement en ce qui concerne la relaxation 
énergétique (déjà traitée par Van Hove dans 
d’autres cas et d’une manière différente), mais 
également en ce qui concerne la relaxation trans- 
versale. Les équations (69) sont valables pour 


1>r et 1>/yH| 


et à la condition que la perturbation V soit sufii- 
samment petite, ou plus exactement à la condition 
(moins restrictive que la précédente) que 


(118) 


Jde OnlPtuila) CIVTON re cime | € [vHol 


<|YH| 
(119) 


= f du (m|V{u)la) (a +n—m|V(0})in) ei®mau 


quels que soient m, a, et n. Physiquement, (119) 
signifie que la largeur d’un niveau d’énergie £,, est 
faible devant l'intervalle (ny ,| qui sépare deux 
niveaux voisins. 

Les équations (69) sont du type (68). Les temps 
de relaxation sont donc les inverses des valeurs 
propres de la matrice À, plus exactement les parties 
réelles de ces inverses. Les parties imaginaires 
correspondent à des déplacements des niveaux 
d'énergie £,, E,, ... dus à la persistance des inter- 
actions spin-réseau. Ces déplacements sont en 
général trop faibles pour être observables actuel- 
lement. La matrice À est du type 
A) 


At-1 
At+1) 


At-o) 
At+ 9) 


I 


A 


AtC-2D 
At+2D 


(120) 


A est la matrice W,,, (71) correspondant à la rela- 
xation énergétique. Cette matrice a été étudiée en 
détail dans la thèse de Lurçat [9]. Les autres 
matrices. AD, .: SAT AE). SAUTER) corres- 
pondent à la relaxation transversale. En examinant 
les coefficients dans (69), on remarque que 


À (MOI) s’exprime à l’aide des seuls élé- 


ments (a|8(t)|b) pour lesquels b — a = n — m. 
A désigne donc la sous-matrice de A reliant entre 
eux tous les éléments de matrice (m|8(t){n) (et leurs 
dérivées) pour lesquels 7 — m = g. Aucune étude 
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approfondie n’a été faite de ces matrices A. La 
matrice A est une matrice carrée de dimension 


21 + 1— gl, il lui correspond par conséquent ; 
214 | g| temps de relaxation transversale 
T9, TP es TP ar 


Mais tous ces temps de relaxation transversale 
n’ont pas le même intérêt physique. Pratiquement, 


toutes les mesures de temps de relaxation trans- M 
versale effectuées jusqu’à présent portaient sur 
transversales 


la relaxation des composantes 
CM > = Yh<L(O> et <M,(E) > = y LD 2 
de l’aimantation totale. À cause des règles de sélec- M 
tion des opérateurs Z, et 1, (Am — + 1) et en 
vertu de l'égalité (32), la relaxation de < M.,(t) > et 
<M,(t) > ne fait apparaître que les temps 


TN TDR, TE (121) « 


avec g — 1( ceux avec g — — 1 sont les mêmes). 


Ce sont les temps (121) qui présentent le plus grand 
intérêt. 

Comme nous l’avons expliqué au paragraphe 4, 
la théorie de la relaxation développée ici est valable 
pour un système de spins nucléaires dans un liquide” 
à condition que les interactions spin-spin entre 
spins nucléaires identiques ne soient pas trop 
importantes, La théorie de la relaxation éner- 


 gétique présentée ici reste valable pour la rela- 


xation quadrupolaire énergétique dans un cristal 
moléculaire ou ionique en résonance nucléaire 
magnétique ; avec quelques modifications, la 
théorie reste également valable pour la relaxation 
quadrupolaire énergétique dans un cristal molé- 
culaire en résonance quadrupolaire. De même, elle 
est encore valable pour la relaxation énergétique 
par interaction avec les électrons de conduc- 
tion en résonance nucléaire magnétique dans les 
métaux, etc... Mais la théorie de la relaxation 
transversale présentée ici ne sera valable pour les 
cristaux moléculaires et ioniques et pour les métaux 
que dans certains cas. Prenons l’exemple d’un 
cristal moléculaire. Alors que dans un liquide, 
l’interaction dipolaire (8) est nulle en moyenne 


(V(t) = 0), elle ne le sera plus en général dans un 


cristal (V(t) -£ 0), où il subsiste par conséquent une: 


interaction dipolaire statique V(t). Celle-ci, à cause 
de son caractère statique. ne pourra pas en général 
(à condition de n’être pas trop importante) induire 
des échanges d’énergie entre les spins et le réseau, 
mais par contre, elle contribuera à la relaxation 
transversale qui n’exige aucun échange d’énergie 
entre les spins et le réseau. Cette interaction dipo- 
laire statique se manifeste encore par l’élargis- 
sement des raies de résonance nucléaire qu’elle pro- 
voque, élargissement qui vient s’ajouter à celui dû 
à la relaxation énergétique (durée finie des niveaux 


Na 


… d'énergie £,, quatrième relation d’incertitude !). 
On peut admettre que notre théorie de la relaxation 
transversale reste valable tant que l’élargissement 
dipolaire statique de la raie demeure faible devant 
celui dû à la relaxation énergétique. Ces restrictions 
ont semble-t-il été perdues de vue par Wangsness 

- et Bloch qui ont postulé pour leur théorie de la 

. relaxation un domaine de validité beaucoup plus 
étendu qu’il n’est en réalité. 

Nous avons remarqué au paragraphe 3 que dans 
Ja théorie de la relaxation présentée ici, le ther- 
mostat apparaissait comme un système classique 
et non pas quantique. Il en résulte que 


Win DE Wym 


comme on le vérifie sur (71) et qu’à l’équilibre 


. thermique, tous les niveaux Æ, sont également 


peuplés de spins : (M]o(co)im) = (n|o(co)|n) quels 
que soint m et n. On peut montrer que l’équation 


de Schrôdir ger (22) d’où nous sommes partis est la 


limite vers laquelle tend l'équation quantique rigou- 
reuse 


lorsque l’on fait tendre %o,,/kT vers zéro, ou 
encore lorsque l’on fait tendre la température T du 
thermostat vers l'infini. Les rapports entre théorie 
quantique et théorie semi-quantique, entre point 
de vue À et point de vue B, ont été étudiés dans un 
article de l’auteur. (La résonance quadrupolaire 
dans les liquides, paragraphe 3, J. Physique Rad., 
1956, 17, 876). La théorie présentée ici est la théorie 
générale de la relaxation d’un système de « parti- 

_cules » non interagissantes (ou faiblement inter- 
agissantes) en contact avec un thermostat à une 
température T'telle que les énergies fo, mises en 
jeu soient faibles devant AT. Il convient de noter 
que les matrices 


of) — (oo) et. B(1) — B(or) 


obéissent aux mêmes équations respectivement (43) 
et (69) que les matrices o(t) et 8(t), ainsi qu’on le 
vérifie aisément. En particulier AN{[(ml£B(t)|m) 
— (mn|B(co)|m)] représente la différence entre les 
populations du niveau Æ,, à l'instant t et à l’équi- 
libre thermique. Les équations de Boltzmann (69) 
dans lesquelles on aura remplacé £(t) par 
[B(t) — B(co)] représentent correctement la rela- 
xation (abstraction faite des termes en %oy/kT 
que l’on néglige puisque %oy, € AT) et elles per- 
. mettent d'exprimer correctement les valeurs ini- 
tiales N[(m|8(0)|m) — (m|B(oo)|m)] des populations 
relativement à leur valeur d’équilibre. Nous pré- 
senterons ultérieurement une extension de la 
théorie exposée ici, extension dans laquelle le ther- 
mostat aura une température finie et qui fournira 
des probabilités de transition W,,, et des popu- 
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lations d'équilibre N(m|B(co)im) convenablement 
pondérées par les facteurs de Boltzmann. 
Wangsness et Bloch traitent le réseau comme un 
système quantique possédant une température 
finie, mais leur travail n’est pas une « théorie de 
l’irréversibilité » puisqu'ils font appel à l’artifice des 
phases (du réseau) réparties au hasard. 


13. Irréversibilité, « Umkehreinwand » et « Wie- 
derkehreinwand ». — Notre démonstration des 
équations (69) a clairement prouvé qu'il n’y a 
aucune contradiction entre les équations (69) de 
Boltzmann qui décrivert la relaxation énergétique 
et transversale et l’équation (22) de Schrôüdinger 
réversible par rapport au temps. En effet, (69) n’est 
équivalent à (22) que compte tenu des conditions 
initiales a) b) c) du paragraphe 6. Ces conditions ont 
été imposées pour que l’on puisse en déduire les 
égalités (41) et (42) essentielles pour établir. (69). 
Pour que (41) et (42) puissent être vérifiées en toute 
rigueur et non seulement approximativement, il 
faudrait que le nombre N de spins nucléaires dans 
le liquide soit infini. Par conséquent, on ne peut 
rigoureusement déduire la relaxation (69) des spins 
nucléaires à partir de l’équation de Schrôdinger (22) 
que pour. un système d’une infinité de spins. Pour N 
très grand, mais fini, les équations (69) constituent 
une bonne approximation, mais elles ne sont pas 
rigoureusement vérifiées puisqu'elles excluent les 
fluctuations des grandeurs macroscopiques (32) 
pendant la relaxation. Pratiquement, on opère 
toujours sur V > 1012? spins nucléaires, 1l est alors 
impossible de mettre en évidence les fluctuations 
parce qu’elles sont trop faibles, et l’on constate que 
les équations (70), conséquences de (69), sont 
toujours bien vérifiées par l’expérience. 

_Il est nécessaire de souligner que la condition (17) 
à laquelle doit satisfaire la perturbation V ne peut 
être rigoureusement vérifiée que pour un ther- 
mostat comportant une infinité de degrés de liberté 
(donc dans le cas où le thermostat est un liquide, 
pour un liquide composé d’une infinité de molé- 
cules). Une assemblée formée d’un petit nombre N 
de spins en contact avec un thermostat infini 
pourra-t-elle se relaxer ? (On pourrait « réaliser » 
expérimentalement une telle situation en mettant 
dans un solvant infini les V molécules porteuses 
des N spins nucléaires en question). Notre théorie 
semble montrer qu’il n’y a pas relaxation des spins 
en ce cas, malgréle fait que le système ! spin + ther- 
mostat } comporte toujours une infinité de degrés 
de liberté. (A première vue, il semble qu'il y ait là 
une contradiction entre les résultats de Van Hove 
et les nôtres). En effet, les égalités (41) et (42) ne 
pourraient être satisfaites pour un nombre N 
de spins qui ne serait pas très grand. Il est alors 
vraisemblable que le système des AN spins décri- 
rait indéfiniment des « cycles de Poincaré » [10, 
#1] tout à fait irréguliers, mais de plus en plus 
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rétrécis lorsque N diminue. On pourra comparer 
ces résultats avec ceux obtenus par J. Eisens- 
tein [12]. Celui-ci a montré par intégration directe 
de l’équation de Schrôdinger qu’un système de 
deux (!) spins couplés par leur interaction dipolaire 
ne se relaxe pas. 

L'idée que seuls pourraient se relaxer des sys- 
tèmes macroscopiques, à très grand nombre de 
degrés de liberté, n’est pas une idée nouvelle. En 
fait, cette idée avait été mise en avant pour per- 
mettre d'échapper au « Wiederkehreinwand » de 
Zermelo [10, 11] basé sur un théorème de Poincaré 
selon lequel tout système fini (à nombre fini de 
degrés de liberté) est quasi-périodique. (Ce sont ces 
parcours quasi-périodiques qui constituent les 
cycles de Poincaré). Mais malgré la reconnaissance 
de ce que la relaxation ne pouvait être que le fait 
d’un système macroscopique, connaissance acquise 
depuis plus de cinquante ans, personne n’avait 
réussi avant Van Hove à prouver, en partant des 
équations de la dynamique (classique ou quan- 
tique), qu’il y a effectivement relaxation en faisant 
tendre N (nombre de degrés de liberté du système) 
vers l'infini. On voit clairement ici comment le fait 
que N doive être très grand s’insère dans le raison- 
nement en tant que condition indispensable. 


14. Le rôle de la mécanique statistique dans la 
théorie des phénomènes de relaxation. — Nous 
désirons attirer l’attention sur le fait que le pro- 
blème de l'irréversibilité de l’évolution du système 
des À spins n’a pas été traité ici par les méthodes 
de la mécanique statistique. On justifie habituel- 
lement [13] l'introduction de méthodes statistiques 
pour traiter les problèmes de relaxation par le fait 
que l’on ne possède de l’état initial de non-équilibre 
du système à étudier qu’une connaissance incom- 
plète, «macroscopique». Ceci signifie que seules sont 
connues à l’instant initial un certain nombre de 
grandeurs macroscopiques < G(0)> (32) relatives 
au système étudié. Dans le cas traité ici d’un sys- 
tème de  spins nucléaires, nous pouvons supposer 
par exemple que seules sont connues à l'instant 
initial les composantes <,(0)>, <M,(0 Lee et 
<M,(0)> de l’aimantation totale. Quant à la 
* fonction d’onde ®(0) du système étudié, elle n’est 
pas connue et ne peut être déduite de la seule 
connaissance des < G(0)>. C’est pourquoi on intro- 
duit un ensemble de Gibbs de systèmes, tous de 
nature identique à celle du système étudié. Chaque 
système S; de l’ensemble possédera une fonction 
d'onde ®,(0) qui, à linstant initial £t — 0, est 
choisie de telle façon à ce que 


(DA(0)|GIDA(0)) = < G(0) > 


c’est-à-dire de façon à être compatible avec la 
connaissance macroscopique que nous avons du 
système étudié à l’instant initial. On admet alors 
que le comportement du système étudié sera vrai- 


(122) 
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semblablement proche du comportement «moyen » « \ 
des systèmes S, de l’ensemble. Mais les conditions » 
du type (122) ne suffisent pas à elles seules à déter- 


miner la répartition dessystèmesS, dans l’ensemble, 
c'est-à-dire à choisir un opérateur densité 
fixer p, il faut introduire des postulats de nature 


statistique, par exemple égalité des probabilités 


a priori de tous les états m compatibles avec (122) . 
et répartition au hasard des phases de tous les 
états m compatibles-avec (122). La discussion de 
ces postulats statistiques est très clairement pré- 
sentée dans le livre de Tolman : The principles of 
Statistical Mechanics [13] auquel nous renverrons 
le lecteur pour une étude approfondie de cette 
question. Seul nous importe ici le fait que l’intro- 
duction de méthodes statistiques pour traiter les 
problèmes de relaxation, c’est-à-dire l'introduction 
d’un ensemble de Gibbs de systèmes S; et d’un 
opérateur densité p? nécessite le recours à des 
postulats statistiques qui ne sont pas réductibles à 
ceux de la mécanique (classique ou quantique). 
L'introduction d’arguments statistiques dans 
l’étude de la relaxation est due à Boltzmann.-Son 
point de vue statistique (évolution du système S 


vers les états les plus probables, etc...) est trop. 


connu pour qu’il soit nécessaire de s’y arrêter ici. 
Ce point de vue garde encore aujourd’hui toute sa 
valeur intuitive, mais il nous semble juste de 
reconnaître qu’il représente une étape dans l’étude 
de la relaxation, étape où il manquait la claire 
conscience de la possibilité de déduire la relaxation 
uniquement à l’aide des lois de la mécanique et 
sans faire appel à des arguments statistiques. 
Certains auteurs, devant la difficulté de déduire 


la relaxation d’un système macroscopique à partir | 


de l’équation de Schrüdinger, et sous l'influence des 
travaux de Boltzmann, Gibbs, Tolman, etc... ont 
admis précisément que ce problème fondamental 
ne pouvait être traité uniquement dans le cadre 
des lois de la mécanique quantique, mais néces- 
sitait en plus pour être résolu l'introduction d’hypo- 
thèses statistiques [14]. La déduction de la rela-' 
xation à l’aide des seules lois de la mécanique 
quantique que nous avons présentée ici, de même 
que la déduction de Van Hove, s'inscrivent en faux 
contre une telle opinion. Non seulement l’usage des 


ensembles statistiques était parfaitement inutile. 


ici, mais il aurait certainement cbscurci la recherche 
des conditions initiales suffisantes pour qu’il y ait 
relaxation. Le théorème 7 déduit ici n’est pas un 
théorème Æ «statistique», c’est-à-dire relatif au 
comportement moyen des systèmes S; dans un 
ensemble de Gibbs, c’est un théorème À «exact» (5). 


(5) Il convient de noter ici qu’un thermostat n’a pu être 
défini jusqu’à présent que moyennant le recours aux pos- 
tulats statistiques mentionnés. De même, les conditions 
(16) et (17) n’ont pu jusqu’à présent être clairement déduites 
à partir de la seule mécanique par la théorie cinétique des 
liquides. Il est cependant possible que le travail de VAN 
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1h Cependant, notre démonstration des équa- 
Ltions (69) n’épuise pas la question de l’irréversibilité 
kde l’évolution du système des N spins nucléaires. 
LE subsiste la question de savoir si d’autres condi- 


tions initiales que les conditions a) b) c) du para- 


faible pour que l’état de régime qui s’établit soit très 
proche de l’équilibre thermique, le spectre d’absorp- 
tion est donné à un facteur multiplicatif près par 


A(o) ju Trace de [Z,{(t) 1,]e—i®t dt (123) 


graphe 6 pourraient également conduire à la rela- 


| xation. Il subsiste également l’importante question 
de savoir pourquoi de telles conditions finissent 


toujours par être réalisées pratiquement (‘). On 


hpeut tirer des conceptions de Boltzmann une argu- 


mentation statistique pour expliquer ce fait, mais 
il reste à en fournir une justification plus précise. 


15. Relaxation transversale et forme des raies de 


résonance nucléaire. — Le champ magnétique de 
.radiofréquence dont on désire étudier l'absorption 
_par le système des NW spins nucléaires étant supposé 


de polarisation rectiligne suivant Ox (Ox étant 
perpendiculaire à la direction Oz du champ magné- 
tique constant H,;) et d'intensité suffisamment 


d io 77! 
LE (IT le) (alT #10] 


1 co ; : 
38 2 ff du EnfP Ca Lai P TON cit mat [ont TT TP IT ] 


Ft sf. du (n|Vlu)la) a] P{0)pn) eérat Cebmnt Gn]T (op) (alT* (ir) ] 
| + CO 
T2 6 J-c0 


Ces équations permettent de calculer la 


| trace (124). À cause des règles de sélection de 


l'opérateur Z,, on devra avoir 
| g = p + 1, Nom EI: 


On voit donc que le problème de la détermination 
de la forme des raies de résonance nucléaire est 


. Hove ouvre la voie à une définition purement mécanique 


. d’un thermostat. Quoi qu’il en soit, l'équation réversible 
. (22) une fois posée est résolue à l’aide des seules lois de la 


. mécanique quantique, aucun ensemble statistique n’étant 


introduit pour représenter les spins. 
(5) Comme nous l’avons vu, l’état initial du système dont 
on désire étudier la relaxation n’est pas complètement 
- connu (parce que le système possède un très grand nombreN 


- de dégrés de liberté). On est ainsi tenté d’attribuer la rela- 


xation à cette connaissance partielle de l’état initial. C’est 
- ce que semble faire R. Brout [16] d’après lequel notre 
ignorance de l’état initial est « incorporée dans une hypo- 
thèse initiale » et « permet ainsi de prendre la moyenne sur 
* un grand nombre de variables ». Un tel point de vue est 
tout à fait inacceptable, La relaxation est un processus 
‘physique objectif, réellement observé, et indépendant de 
notre connaissance du système à l’instant initial. Nos condi- 


_ tions a, b) et c) du paragraphe 6 suffisantes pour qu’il y ait 


_ relaxation ne sont nullement incompatibles avec une con- 
naissance microscopique complète de l’état initial (aussi 
complète que le permet la mécanique quantique). Si nous 
possédions d’un système macroscopique une connaissance 

* microscopique compiète à l’état initial, pourvu seulement 
que celui-ci vérifiât les conditions a), b) et c) du para- 
graphe 6, il y aurait relaxation. Nulle part ici, nous n’avons 
moyenné sur notre ignorance | 


où J, désigne la composante du vecteur spin sur 
l’axe Ox, et où Z,(t) désigne cette même compo- 
sante en représentation de Heisenberg 
Talt) = T*(:) To T(t) 
T(t) étant l’opérateur (44). Mais 
Trace de [7,(11,] 


= nn TOP QT) (alZulre) (plLle (124) 


On voit ainsi réapparaître les expressions 
(mT(t)|p) (g\T*(t)n) qui intervenaient dans la 
théorie de la relaxation. En vertu des égalités (43), 
(63) et (69), on aura (avec q9 = p + n-— m) (?) 


du (m] P{u)[a) (a + n — m7 (0)pn) mal etémnt (a T()|p) (glL* (bla + n — m)] 


équivalent à un cas particulier du problème général 
de la relaxation transversale, à savoir à la rela- 
xation transversale de quantités <G(t)> G ayant 
la règle de sélection Am = + 1. 

La raie, associée aux transitions Am = + 1, 
sera donc en général une superposition de courbes 
de Lorentz, les paramètres T, des diverses 


courbes étant les temps de relaxation T7, 
SN ACTA EE 

A Ÿ - 126 

ETES — mi { 2 ) 


Les T peuvent donc être interprétés indif- 
féremment comme des temps de relaxation trans- 
versale associés aux grandeurs G dont la règle de 
sélection est Am — - 1 ou comme des para- 
mètres 7, relatifs à la largeur des raies (126). Nous 
renverrons pour l’étude approfondie des formes de 
raie à l’article d’Ayant [8] qui est le premier à avoir 
donné une théorie correcte de la forme des raies de 
résonance nucléaire dans les liquides. 

La théorie de la relaxation développée ici est un 
cadre général qu’il faut remplir en fixant d’une 


(7) Les expressions (m|T(t)|p) (g/T*{t)In) ne sont sensi- 
blement différentes de zéro que pour g = p + n — m. 
Cela a été démontré par Ayant [8] au second ordre en W. 
Cela est également visib'e sur (60) au second ordre en W. 
Nous laissons la démonstration générale de ce fait (elle 
n’est pas difficile !) aux soins du lecteur. 
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part la nature du couplage V et en explicitant 
d’autre part la nature des mouvements molé- 
culaires dans le thermostat. Dans un prochain 
article de cette série, nous développerons la théorie 
de la relaxation quadrupolaire dans les liquides 
dont un résumé suceinct a paru : [17]) et nous 
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confronterons la théorie avec les résultats expé- 
rimentaux (5). 

(8) La théorie exposée en détail dans cet article a été 
brièvement résumée dans trois notes de l’auteur: [6] et [15]. 


Manuscrit reçu le 28 décembre 1956. 
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RE 


REVUE DES LIVRES 


L’énergie atomique dans ses répercussions sur la vie et la 
santé. (1 vol., 13,5 X 21 cm, 254 pages, Édition de 
l’Expansion, Paris, 1956, 1 000 F.) 

C’est le compte rendu d’une conférence scientifique qui 
s'était réunie au Muséum National d'Histoire Naturelle 
les 127 et 2 juillet 1955. Des physiciens atomistes, des 
géophysiciens, des météorologistes, des biologistes, des 
généticiens et des médecins y ont mis en commun les 
* résultats de leurs recherches et de leurs travaux de mises 
au point, sur des questions vitales pour l’avenir de l’espèce 
humaine. 

Les principaux rapports sont les suivants : Dangers liés 
à l’équipement atomique, par F. Netter. — Effets radio- 
actifs immédiats et différés des explosions atomiques, par 
R. Chastel. — Analyses des cendres radio-actives projetées 
lors de l’explosion de Bikini du 1€* mars 1954, par J. Laber- 
rigue-Frolow et T. Yuasa. — Propagation à grande dis- 
tance et caractéristiques des particules radio-actives issues 
des explosions atomiques, par H. Garrigue. — Influences 
éventuelles des explosions atomiques sur l’évolution du 
temps, par J. Roulleau. — Mesure de la radio-activité de 
l’air et de la pluie, par O. Tanaevsky et E. Vassy. — Les 
nuages radio-actifs, par J. Orcel. — Exposé général sur les 
effets biologiques des radiations ionisantes, par A. Lacas- 
sagne. — Problèmes médicaux posés parles effets immédiats 
des explosions atomiques, par P. Vigier. — Lésions ocu- 


laires consécutives aux explosions ou aux recherches 
d'énergie atomique, par M. A. Dollfus. — Radiations 
atomiques et vie aquatique, par M. Fontaine. — Danger M 
biologique des poussières radio-actives émettrices de rayons 
bêta, par M. Lefort. — Le problème de l’action biologique. 


des faibles doses de radiations ionisantes, par H. Marcovich. M 


— L'énergie ionisation et l’avenir du patrimoine-héréditaire 


humain, par R. Turpin et J. Lejeune. — Protection contre LM 


les rayonnements, par P. Savel. — Les aérosols et gaz! 
radio-actifs dans l’industrie atomique, par J. Labeyrie. — 
Protection thérapeutique des radio-lésions, par H. Marco- 
vich. 

En outre sont reproduits les exposés de A. Haddow 
(Grande-Bretagne) et E. Rabinowitch (U. $S. A.) sur les 
responsabilités humaines du savant, ainsi qu’une partie ” 
substantielle des discussions qui furent soulevées à propos \ 
des effets météorologiques des explosions, de l’efficacité des 
moyens de protection et des mutations génétiques induites 
par les radiations atomiques. 

A part l’étude des périls spectaculairés causés par Îcs 
explosions nucléaires, il faut signaler les développements M 
importants qui ont été consacrés aux dangers liés à la 
manipulation des radioéléments ou des autres sources de 
radiations ionisantes et qui ne manqueront pas d’attirer 
l’attention de nombreux physiciens. 
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STRUCTURES HYPERFINES DE RAIES DU SPECTRE D'’ARC DU MERCURE 
ET MOMENT QUADRUPOLAIRE DE ?01Hg 


Par JEAN BLAISE et 


Henrt CHANTREL, 


Laboratoire Aïmé-Cotton, G. N. R. $., Bellevue (Seine-et-Oise). 


Introduction. — Nous avons entrepris en 1952 de 
vérifier les performances des interféromètres Fabry- 
_ Perot à couches diélectriques ; depuis cette époque, 


_ chaque amélioration apportée aux spectromètres 


photoélectriques à un seul étalon [1] et à double 
étalon, a été contrôlée en étudiant la structure 
hyperfine d’une raie du mercure. Etant donné le 
grand nombre d’auteurs qui se sont intéressés à la 
structure hyperfine des raies du mercure, (il y a eu 
70 publications sur ce sujet depuis 1930), nous 
avons été supris de constater le besoin d’une révi- 
sion systématique de certaines structures déter- 
minées antérieurement. Grâce aux nombreux enre- 
gistrements de raies ainsi effectués, nous avons 
maintenant à notre disposition toutes les données 
expérimentales nécessaires à la détermination des 
moments nucléaires des isotopes impairs du mer- 
cure et, en particulier, du moment quadrupolaire 
de l’isotope 201. 

Le moment quadrupolaire de ?1H5 a été déter- 
miné pour la première fois en 1935 par Schüler et 
Schmidt [2] à partir des constantes de couplage 
quadrupolaire B ‘des termes 6p %PS et 6p $P? du 


… spectre d'arc du mercure ; ils obtinrent un moment 


quadrupolaire Q = 0,5 X 10-24 cm?. En 1950, 
Murakawa ‘et Suwa [3], utilisant le terme 6p 5P9, 
ont obtenu Q = 0,2 X_ 10-24 cm2. De leur mesure 
de la résonance quadrupolaire nucléaire d’un cristal 
de HgCl,, Dehmelt, Robinson, et Gordy [4], en 


… 1954, ont déduit un moment quadrupolaire de 


- 0,6 X 10-24 cm?. Enfin, récemment, Murakawa [5], 
ayant amélioré sa mesure de la structure hyperfine 
du terme Pet tenu compte dans le calcul de Q des 
interactions de configurations électroniques, a 
conclu à une valeur de 0,42 X 10-24 cm?. 

De telles divergences n’ont rien d’étonnant 
quand on songe d’une part à la difficulté de la 
détermination des constantes de couplage B qui 
sont très sensibles aux erreurs de mesure et d’autre 
part aux hypothèses intr duites dans le calcul de Q 
à partir de B, que ce soit en spectroscopie optique 
ou, surtout, en radiofréquence. Un des buts de ce 
travail a donc été la détermination aussi précise 
que possible des constantes de couplage B des 
termes P° de la configuration 6s6p de Hg J. 
L'accord entre les valeurs de Q, déduites de ces 
constantes B, est susceptible d’être amélioré par 

_ un raffinement des calculs. 


Dans tout ce qui suit, les notations utilisées sont 
celles de Kopfermann [6]. 

L'énergie W# des niveaux hyperfins d’un terme 
spectral d'énergie W; est donnée par la formule de 
Casimir [7] : 

? CIC HA NE 1) AJ HA) 
21(21—1) J(27 +1) 
où CPR LAN TIEAE 4) JUNE) 


F est le nombre quantique hyperfin, J le nombre 
quantique interne et J le spin nucléaire. Le facteur 
d'intervalle magnétique À et la constante de cou- 
plage quadrupolaire B seront exprimés en milli- 
Kkaysers (1 mK = 107$ em-1. 
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Fic. 4. — Schéma de niveaux du spectre d’arc du mercure. 


La figure { donne le schéma des niveaux qui 
fournissent les transitions les plus importantes 
de Hg Z. On voit que pour la détermination du 
moment quadrupolaire, les termes intéressants sont 
ceux de la configuration 6s6p. Parmi eux, il n’y à 
que le terme #P9 qui soit indépendant du couplage 
et qui permette donc un calcul simple de Q à 
partir de B. C’est pour cette raison qu'il a été 
choisi par Schüler et Schmidt. Nous verrons qu’on 
peut déterminer directement les trois intervalles 
hyperfins de ce terme à partir de la raie À 5 461 
(6s6p 3P9 — 6s7s 3%$,) et qu’on obtient ainsi une 
valeur précise de B. L'interprétation de la raie 
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verte nécessite au préalable la détermination de la 
structure hyperfine du terme 3S, que l’on déduit 
de la raie À 4 047 (6s6p PS — Gs7s 3S,). 

La structure hyperfine du terme %P? est 
mesurée directement sur la raie À 4 078 
(6s6p #P? — 6s7s 1S,) dont toutes les composantes 
sont parfaitement résolues. 

Le raie À 4 358 (6s6p $P, — 6s7s 389) permet un 
contrôle très satisfaisant des déterminations 
séparées des structures hyperfines des termes 3P9 
et %5:. 

Quant au terme 6s6p 1P9, la détermination 
de sa structure hyperfine devrait pouvoir se 
déduire immédiatement de l’étude des raies 
À 4 916 (6s6p 1P9 — 6s8s IS) ou À 4 108 (6s6p 1P9 
— 6s9s LS;) qui sont identiques, au déplacement 
isotopique près. Malheureusement, la composante 
médiane de l’isotope 201 est masquée par la plus 
intense de l’isotope 199. Nous indiquerons un pro- 
cédé basé sur des considérations de déplacements 
isotopiques qui permet d'estimer la position de 
cette composante et par suite d’évaluer B. 

Murakawa et Suwa ont essayé de déduire la 
structure du terme 1P9 de la raie À 5 790 (6s6p 1P9 
— 6s6d 1D,). Cette raie est perturbée par À 5 789 
(6s6p 1P? — Gs6d #D;) et de plus, les composantes 
intéressantes pour le calcul de B sont très faibles 
(1,10 et 0,22 %) et leur distance est inférieure 
à 50 mK, soit moins de 1 /25 de la structure totale, 
ce qui limite sérieusement la précision de la mesure. 

Nous avons mesuré la raie À 4 347 (6s6p 1P9 
— 6s7d 1D,) qui présente une structure hyperfine 
du même type que A5 790. 

Il s’est avéré que le terme 7d 1D, était perturbé 
par le terme 7d 1D;, cette perturbation étant ana- 
logue à celle des termes 6d 1D, et 6d 1D;, bien que 
plus faible. Dans ces conditions le contrôle de la 
structure du terme 6p 1P9 à l’aide de cette raïe est 
assez aléatoire. 


Résultats expérimentaux. — Les enregistrements 
des raies ont été effectués avec des épaisseurs 
d’étalon comprises entre 1,132 et-65,250 mm (1). 
Les grandes épaisseurs ont été mesurées à l’aide 
d’un comparateur, les faibles à l’aide des raies 
À 4 358 et À 4 347 de Hg. 

Nous avons utilisé sept types de sources de 
lumière dont la puissance allait de 1,5 W à 2 KW. 
Les positions des composantes de structure hyper- 
fine sont identiques, quel que soit le type de lampe, 
mais leurs intensités relatives ne sont les intensités 
théoriques que pour les lampes fonctionnant sous 
faible pression. Nous ne parlerons dans la suite que 
des enregistrements effectués avec une cathode 
creuse de Schüler refroidie à l’azote ou à l’hydro- 


(2) Sur l'enregistrement de la raie À 4 358 effectué avec 
une épaisseur de 1,132 mm, les composantes g et e dis- 
tantes de 69 mK sont encore séparées, ce qui indique une 
finesse effective de 65. 
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gène liquides, le néon ou lhélium servant à produire 
la décharge. Nous avons également obtenu de bons - 
enregistrements de raies faibles avec un tube scellé 
contenant de la vapeur de mercure pure et excité M 
au moyen d’un générateur H. F. à magnétron « 
(« Raytheon microtherm unit ») de 2 450 Me/s, en 
n’utilisant que 10 à 15 % de la puissance disponible 
et en refroidissant le tube avec un ventilateur. 

D’après Nier [8], les abondances relatives des 
isotopes du mercure sent les suivantes : 


TABLEAU I 
Isotopes .. 196 198 199 200 261 202: 204 
Abondances ÿ 
relatives.. 0,146 10,02 16,84 23,13 13,22 29,80 6,85 


Chaque fois que cela a été possible, nous avons 
pris comme origine des structures hyperfines la 
composante de l’isotope 200. Presque toutes les 
raies étudiées ayant déjà été analysées par Schüler 
et ses collaborateurs [1], [9], [10], nous avons 
adopté leurs notations. ; 

Nous indiquons pour chaque raie la structure 
chservée, la position des centres de gravité des 
isotopes 199 et 201 et les schémas de niveaux 
hyperfins correspondants. 

À moins d’indication contraire, l’erreur probable 
sur la position des composantes est inférieure M 
à 1 mK. 


Rate À 4047 À (6s6p 5P9 — 6s7s 35). — Un 
enregistrement fait avec une épaisseur d’étalon: 
de 3,836 mm permet la détermination directe 
du terme 7s %S,, toutes les composantes desisotopes 
impairs étant parfaitement résolues. Au cours 
d’un enregistrement fait avec une épaisseur de M 
5,146 mm, nous avons substitué à la cathode 


-739,6 389,4 “à S8&8 274,0 335,8 672,0 6 
2 29,4 
C.g. 199 =_22,7 
C. g: 201 = + 8,6 
Hg 199 Hg 20! 
%k V 
7s $s, 398,0 
3 
72 1061,4 
Ve 
B a b c 
5,61 11,23 6,61 4,41 2,20 
V2 6p P, 


Fic. 2. — À 4 047 À (6s6p 5P, — 6578 #S;). 


N° 3 


- creuse une lampe à Hg”, ce qui nous a permis de 
- repérer la position de cet isotope par rapport aux 
composantes a et B. Enfin avec une épaisseur 
- d’étalon de 22,105 mm, nous avons pu séparer les 
isotopes 200 et 202 et mesurer leur distance à la 
composante B. Aucune expérience ne nous ayant 
permis de séparer l’isotope 204, nous avons admis 
que les isotopes 200, 202 et 204 étaient équi- 
distants. 
|: Les figures 2 rassemblent les résultats ainsi 
obtenus. On voit que le terme 7s7s %S, obéit parfai- 
tement à la loi d'intervalle de Landé et qu’il a donc 
la symétrie sphérique, comme l'exige la distri- 
bution de sa densité électronique. 


Rate À 4 078 À (6s6p 3P9 — 6s7s 18,). — Toutes 
les composantes impaires étant résolues sur un 
enregistrement fait avec le tube excité en haute 
fréquence et une épaisseur d’étalon de 4,770 mm 
(fig. 3a), on a immédiatement la structure hyper- 

| fine du terme 6p #P£. 

En raison de l'intensité plus faible de la raie 
à 4 078, le déplacement isotopique 200-202 n’a pas 


Do me 


Fic. 3a : Tube excité en H. F. Épaisseur d’étalon : 
e = 4,770 mm. 


-434,S -270)-182,7 Je\s 283,7 467,7 OC 


-28,3 31,4 


C.g. 199 = _ 24,2 
C.g.201 =+ 8,5 


755, % 


Ant Dose C © 


À B 
11,23 5,61 2,20 4,41 6,61 
* ; 2518 
É 
737,8 6p °R° le me Âres 
V2 . 


| Fic. 36 et 3c : Structure fhyperfine et schémas de niveaux. 
Fc. 3. — À 4 078 À (6s56p SP, — 6s7s 150). 
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été mesuré avec une précision aussi grande que 
pour À 4 047. 

Comme pour À 4 047, la substitution en cours 
d'enregistrement de la lampe à “SHg à la cathode 
creuse a permis de repérer la position de cet isotope 
par rapport aux composantes impaires 


(4, 198 = 242 + 1 mK). 


Raïe À 4 358 À (6s6p SP9 — 6s7s 38). — Cette 
raie a été enregistrée avec des épaisseurs d’étalon 
très diverses et seule la composante d (0,59 %) n’a 
pas pu être isolée. 


198 C 
= b 204 
A a 
I î 196 d % QUE D 
— 986,6 -580,4 -248,8 co =] a a À 229,2) 480 5494 266 G 
C- g. 199 =-23, 140 270-300 88.8 
C.g. 201= 9,0 
fe F 
PA SM ln 2 
398,0 
3, 4 D 
10754 7SNS,) | 1061,4 
663,4 | ë 
Ve #2 
ARE AE ER atMbmretidWe ta 
1,98 4,63 184 0.59 1,98 087 84 
1,87 3,74 9.36 1,87 | 
3 % 
’e M, 251,8 [ 
737,8 6p “A dé 4664 | 718,2 
Ve ee % 


Fic. 4. — à 4 358 À (6s6p 5P, — 6s7s 8S;). 


L'accord avec la structure calculée à partir des 
termes 3P© et %S, est excellent et nous consi- 
dérerons donc comme très sûre la valeur de B tirée 
du terme %P9. 

L'isotope 196 (0,15 %) est visible sur un enre- 
gistrement effectué avec la cathode creuse refroidie 
à l'azote liquide (intensité du courant : Z — 8 mA) 
et une épaisseur d’étalon de 22,105 mm. 


Rare À 5 461 (6s6p PS — 6s7s %S,). — Cette 
raie a été enregistrée avec un étalon double, la 
cathode creuse étant refroidie à l'hydrogène liquide 
(1 = 7 mA). Nous avons constitué ce double étalon 
d’un étalon à 7 couches diélectriques d’épasseur 
4,850 mm et d’un étalon à 5 couches d’épaisseur 
34,800 mm ou 65,250 mm, soient les rapports, 8 ou 
15: 

D’autres enregistrements ont été faits avec un 
étalon simple d’épaisseur 22,105 mm et une cathode 
creuse refroidie à l’azote liquide (7 = 8 mA). 

Seules n’ont pas été résolues les compo- 
santes B (10,11 %) et e (1,41 %). L’intervalle (e, f) 
étant identique à l'intervalle (a, b) la position de e 
$e calcule immédiatement. 


196 


La distance de B à la composante 202 est au 
maximum de 0,5 mK, B se trouvant du côté de la 
composante 200. 


e+rc B 


Xg 2 AUTRE A 


DAC A 


Fic. 5a: Cathode creuse refroidie à l’hydrogène liquide 

INT nn ASTRES 

Double étalon : (es — 65,250 mm et e, — es/15.) 

Les composantes paires ont été enregistrées avec une 
amplification 4 fois plus faible que celle utilisée pour le 
reste de l’enregistrement. 

Les 2 composantes désignées par « sont des ghosts des 
composantes 200 et 202 + B. 


a b d el 
=724,7 ATX] -2876 -231,3 a] tea 
C.g. 199 =-_ 25,0 =574" gs 10 
C.g. 201=+ 9,1 


1075,4 % 


EE O AELO AREA TC 
Of8 119 £,29 0,22 1,41 2,78 110 1/0 


| 


2 
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Fic. 5b et 5c : Structure hyperfine et schémas de niveaux. 
Fic. 5. — À 5 461 À (6s6p P; — 6s7s S;). 


POSITIONS RELATIVES DES ISOTOPES. — Les raies 
précédentes ont été mesurées avec une très grande 
précision et nous permettent de tirer des con- 


TABLEAU II & 


VALEUR ABSOLUE 


À IDU DÉPLACEMENT 196 198 199 200 201 202 
| (200-202) 
| 047 29,4 — 0,90 — 0,77 0 0,29 1 
4 078 30 — 0,90 — 0,79 0 1027 1 
4 358 30,0 4,64 — 0,90 — 0,79 0 0,30 1 
5 461 31,0 — 0,904 — 0,807 0 0,29 1 
[Valeur 
moyenne — 1,64 — 0,90 — 0,79 0 0,29 1 
TABLEAU II b 
VALEUR ABSOLUE 
À DU DÉPLACEMENT 198 199 200 201 202 204 
(200-202) 
15 295 143,1 — 0,89,  — 0,804 0 0,30 1 2,0 
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clusions sur les déplacements isotopiques relatifs 
En prenant comme unité le déplacement isotopique 
(200-202) dans chaque raie, on obtient les résultats 
suivants (Tableau IT a). J 
Pour la raie À 15 295 À (6s7s 35, — 5d6s26p SP9}, 
les déplacements relatifs mesurés étaient les sui 
vants [11] (Tableau II b). . 
C'est cette constance des déplacements 1is0 
topiques relatifs dans les raies non perturbées quM 
nous à permis par extrapolation d'interpréter les 
raies-à 4 916 et à 4 108 À dont il est impossible de 
résoudre la structure si l’on ne dispose pas d’iso= 
topes enrichis. FES 
14916 À (6s6p 1P9 — 6s8s 1S,) et à 4 108-ÀN 
(6s6p 1P9 — 6Gs9s 15,). — Ces deux raies ont leuw 
terme inférieur commun et leur terme supérieur“ 
n’a pas de structure hyperfine. Elles ne peuvent, 
donc différer que par leur déplacement isotopique 
Ces raies sont peu intenses en cathode creuses 
mais s’excitent très bien à l’aide du géné” 
rateur H. F., comme le montre la figure 6& où I& 
résolvance effective est de 930 000, l'épaisseur de« 
l’étalon étant de 18,205 mm. | 
De l’intensité et de la largeur de la composante 
complexe B + b, on peut conclure que B et hr 
coincident à mieux que mkK près, dans la 
raie À 4 916. On en déduit le schéma de niveaux de 
la figure 6c. 
Les structures de la figure 6b ont été mesurées à 
partir du milieu à mi-hauteur de la composantes 
complexe formée par les isotopes pairs. C’est pratis 
quement le centre de gravité des isotopes pairs qui 
d’après le tableau Ila, doit se trouver à une dis 
tance y de l’isotope 200 telle que y — 0,49x 0 
x = (200-202). | 
Pour À 4 916, la position du centre de gravité» 
de Hgl® par rapport au centre de gravité des 
isotopes pairs est — 7,60 mK. On a done# 
(0,79 + 0,49) x = 1,28 x = 7,60. A 
D'où x = 5,93mKety = 2,99mK. F 
Le centre de gravité de Hg*1, en supposant 
comme nous l’avons fait, que les composantes b 


w 
z 

L 
D 


N°3 


et B coïncident, se trouve à — 1,17 mK du centre 
de gravité des isotopes pairs On a donc: 
DD 2018=0043 — (0,79 + 0,29) x — 1,08x. 
bous=5,95 mK. 


LR RLENERN 


0 B+#b c 


Fic. 6a : Tube excité en H. F. Épaisseur d’étalon : 
e — 18,205 mm. Résolvance effective : R — 930 000. 
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FIG. 6. — À 4 916 À 


He) 4108 À 
(656p 1P, — 6s8s 1). 


(6s56p 1P, — 659s 15,). 


L'accord entre les deux valeurs de æ semble 
justifier le procédé, qui est confirmé par ailleurs par 
le fait suivant : l'intervalle (198-202) calculé est 
de 1,90x, soit 11,3 mK, alors que le même inter- 
valle mesuré directement par McClung et 
Holmes [12] à l’aide d’isotopes séparés, est de 
13 + 1 mK. 

La structure calculée de À 4 916 est donnée par 
la figure 6d. La structure de la raie À 4 108 s’en 
déduit facilement en observant que l’inter- 
valle (A, a) est de 55,2 mK au lieu de 58,2 mK 
dans À 4916. Il en résulte que l'intervalle 
(199-201) = 6,43 — 3,0 — 3,43 mK et par suite 
que (200-202) — 3,2 mK. On en déduit la structure 
indiquée par la figure 7d, structure qui est bien 
confirmée par un enregistrement effectué avec un 
étalon de 39,31 mm d’épaisseur. 
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Trois autres raies bleues ont pour terme infé- 
rieur 6p !P1, ce sont : À 4 347, À 4 343 et À 4 339. 

À 4 347 (6s6p LP9 — 6s7d 1D,). — La figure 8a 
montre un enregistrement effectué avec un étalon 
d'épaisseur 5,146 mm (?). Outre les composantes À, 
B,C,aet h déjà mesurées par Schüler et Jones [13], 
nous avons pu déterminer les positions de d et g. 
La composante d est très bien résolue sur un enregis- 
trement effectué avec une épaisseur de 18,145 mm, 
Ces sept composantes permettent, connaissant la 
structure hyperfine du terme 6» 1P9, de calculer 
celle du terme 7d 1D, et par suite la structure de la 
raie. 


nee 


a B bO dc Cegfh 


+= 
A 


Fic. 8a : Cathode creuse : 1 — 25 mA. Épaisseur d’étalon: 
e — ,146 mm. Résolvance effective : R — 815 000, 
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F1G. 8b et 8c: Structure et schémas de niveaux. 


Fic. 8. — à 4 347 À (6s6p 1P, — 6s7d 1D;). 


À 4 343 (6s6p 1P9 — 6s7d D). — Cette raie est 
malheureusement masquée par un ghost de À 4 358, 
dû au réseau utilisé comme monochromateur, et 
d'intensité sensiblement double. La composante 
complexe des isotopes pairs a la même largeur que 
celle de À 4 347 et une composante impaire intense 
se trouve à environ + 40 mK. 


(2) Nous tenons à remercier M. le P' D. A. Jackson de 
nous avoir prêté cet étalon dont les lames sont remar- 
quablement planes. 


198 
À 4 339 (6s6p 1P9 — 6s7d $D,). — Cette raie à 


été enregistrée avec des épaisseurs d’étalon 
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C.g. 201 =- 1,6 


Fic: 9. — à 4 339 À (6s6p 2P, — 6s7d $D;). 
En haut à droite lire æ 45 au lieu de = 46. 


de 5,105 mm, 18,145 mm et 39,31 mm. Sa structure 
n’a été qu'imparfaitement résolue et les résultats 
indiqués sur les figures 9a et 9b ne sont que pro- 


al 


cal. 
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visoires pour l’isotope 201 dont seules les compo- 


santes b et c ont été observées, donnant une sépa- 
ration de 75 mK entre les niveaux F = 38/2 et 
F = 5/2 du terme 7d%D,. Avant de terminer 


Pirterprétation de ce terme, il est nécessaire de. 


considérer les résultats obtenus pour le terme 7d1D, 
à partir de À 4 347. 


Perturbation des niveaux hyperfins des termes 
de la configuration 6s7d. — La théorie de l’inter- 
action quadrupolaire permet de prévoir pour les 
termes de la configuration 6s7d des facteurs de 


couplage B très petits, de l’ordre de 0,05 mK. On" 
devrait donc s’attendre à ce que les niveaux hyper-" 


N°3 


TI WT. 


ET PS TS a 


fins de 7d 1D, obéissent pratiquement à la règle 


d'intervalle. Or c’est loin d’être le cas. L’appli-" 
cation de la formule de Casimir aux niveaux hyper- 
70,2 mK et" 
B = 6,4 mK, valeur cent fois supérieure à celle“ 


fins de 7d1D, donne : A?01 = — 


prévue par la théorie. 
De plus, alors que pour les termes des confi- 


gurations 6s7s et 6s6p, le rapport des facteurs 


d'intervalle des deux isotopes impairs est sensi- 
blement constant et égal à — 0,370, il est nota- 
blement inférieur pour le terme 7d 1D,, comme 
l'indique le tableau III. 


TABLEAU III 


CONFIGURATION 6575 
Terme 355 DRE 
ANSE 716,93 302,92 
A201 — 265,35  — 111,85 
A201 F 
PTT — 0,37014 — 0,3692 
B 0 14,02 


On se trouve donc en présence d’une pertur- 
bation par le terme 7d D, distant de 20 K, anaa- 
logue à celle des termes 6d 1), et 64 3D, décou- 
verte par Schüler et Jones [13] et interprété par 
Goudsmit et Bacher [141 et par Casimir [6]. Suivant 
ces auteurs, lorsque la structure: hyperfine est de 
l’ordre de gran deur de la structure fine, des niveaux 
ayant la même valeur de F se repoussent l’un 
l'autre. Les niveaux F — 3/2 de l’isotope 199 et 
F=5/2, F = 3/2 et F = 1/2 de lisotope 201 
sont repoussés par les niveaux correspondants 
de °D, comme le montre la figure &8c. Si le 
terme 7d 1D, n’était pas perturbé, le rapport des 
structures totales des deuxisotopes Ao(201)/Ac(199) 
serait égal à — 1,11. La théorie de Goudsmit et 
Bacher permet de prévoir un rapport de — 1,08 et 
l'expérience donne — 1,07. La position anormale 
des centres de gravité des isotopes 199 et 201 par 
rapport au centre de gravité des isotopes pairs 
con firme cette interprétation. 

Le terme 7d%D, doit être perturbé par le 
terme 7d D, distant de 22 K et on peut également 


6s6p 6s7d 
Ge: LP: ÊDS $D; 
491,87 — 119,20 198,0 — 87,2 
— 181,89 43,91, - — 70,2 
— 0,3698  — 0,3684 — 0,355 
— 9,35 8,65 (6,4) 


prévoir une diminution du rapport Ac(201) /Ac(199). 
Pour le calcul de la structure de la raie À 4 339 
nous avons choisi Ac (201) /Ao (199) — — 1,07. 
L'expérience He a (199) 
D'où : Ao (201) — 233 + 5 mK. On mesure de 


plus G : = uk 


Les deux intervalles sont estimés. La structure 


— 218+2mK. 


ainsi calculée est justifiée par le fait que la position « 


anormale des centres de gravité des isotopes 


impairs est analogue à celle de la raie À 4 347, ce 


qui est conforme à la théorie de Goudsmit et Bacher. 


Rapport des moments magnétiques des isotopes 
impairs. — Ce rapport se déduit du rapport des 


facteurs d'intervalle du terme 6s7s %9, : (cf. 
-tableau IIT). 
u201 201 4201 
MEL = D 08 = 1540 SU En 11103. 


La valeur de Schüler et Schmidt était de : 
— 1,1089, 


_ Nos 
Er 


x 


? _ Calcul de ul*® à partir du terme 6s7s°S,, — 


On a : : 
1 
} A(S:) — 3 (dés + @7s) = 716,93 mK. 
; _ Sil’on admet que : 
à dna(75) 
Gzs _ Na°(65) dn . 
male ans 009725 
dn 


On a : ay = 1 306,8 mK. 
En portant cette valeur dans l’expression : 


4 De nel : 
RER A1 7) + ARCS 
| où la correction relativiste Æ,{1 /2, 80) — 2,2573, 
la correction de Breit-Crawford-Schawlow : 
1 — 5 = 0,88 et la correction de Bohr-Weisskopf : 
1 — € — 0,97. 
On obtient : ul® — 0,532 magnétons nucléaires. 
La valeur trouvée par Proctor et Yu [15] est 
- uw = 0,4994 + 0,001 sans correction diamagnétique 
» et p —= 0,5202 avec la correction de 0,96 % pro- 
posée par ces auteurs (5). 


Caleul du moment quadrupolaire de 2!Hg. — 

- Le moment quadrupolaire est lié à la constante de 

couplage quadrupolaire d’un terme spectral par la 
formule : 


B = eQ. ogr(0) 


- où oyy(0) est le gradient du champ électrique 
- produit par les électrons à l'endroit du noyau dans 
la direction du vecteur J. 
Nous avons déterminé B pour tous les termes de 
la configuration 6s6p (cf. tableau III). 
. Le terme °P5 est indépendant du couplage et 
on à : 


2 ENT 
AUDE 


RE, 7) 


| avec 
ces 5W 
PA à uB(2l + 1),2; H,(l, Z;) 


où 3w est la distance 3P2 — 3P9 avec la correction 
de Wolfe qui tient compte de l'interaction spin- 


_ orbite (3w = 5 622,6 K). 


o=8.+1)@+ LA tee INR 
| ,., 21 6) Sÿ.RAI, Zi) 
(2) — 0,37276.10—2 C. G. S. 
Zee 
H,(1.76) — 11632;  R’,11,76) = 1,2892 et 


B = 14,02 .10—$, 


() Avec la correction diamagnétique calculée par Lamb, 
! on obtient : u — 0,5043 magnétons nucléaires. 


STRUCTURES HYPERFINES DE RAIES 


D'où : 
Q'= 0,478,10—24 cm?. 
Pour les termes #P9 et 1P9, on doit appliquer les 
formules du couplage intermédiaire de Breit et 
Wills (16). On a : 


oJJ(0) = €. ee Cr Ra4, 76) — aus, (176) | 
D'où : 
0 2: Lane : 0,32276 OS 
SW CR, —92V3cc 8 
OÙUS 15084 


c, et c, sont les constantes de couplage qui per- 

mettent d'exprimer la fonction propre d’un terme 

J — I, en couplage intermédiaire, à partir des 

fonctions propres des électrons. Dans notre cas : 
43 pl 


de des) + a de(s) 


Dans le cas du terme #P9 : 


Din OS LOTS 
ce = 0,4419 : ec, — 0,8974 
Q0="0;,546.10=22çcm?, 


Dans le cas du terme 1P£ : 


B = 8,65.10—3 
c = 0,8974 ; © = — 0,4412 
ON MIS eme 


Discussion des résultats. — La valeur de Q 
déduite du terme 1P9 est nettement inférieure aux 
deux autres et l’erreur possible sur B ne peut pas 


expliquer cet écart. Comme le terme !P? a un 


déplacement isotopique inférieur à celui des autres 
termes de la configuration 6s6p, il est possible qu'il 
soit perturbé par le terme complexe 54°? 6s? Gp !P° 
qui a un B négatif. La diminution de la valeur de B 
du terme 6s6p 1P£ proviendrait donc en partie de 
cette perturbation mais ce n’est sûrement pas 
l'unique cause de la petitesse de la valeur de Q 
déduite de ce terme. 

Il semble que les divergences doivent être attri- 
buées au mode de caleul de 7,(0). En effet, d’après 
la théorie de Breit et Wills, les facteurs d'intervalle 
des térmes de la configuration 6s6p sont donnés 
per les formules : 


a 3 
A(P3) = % + © Gpair 


1 # = 
= 7 L(2c? — C2) ag “+ 50? Apapo + 203 Apype + EC Ca V2 a] 
+ 


où €, el c, ont respectivement les valeurs indiquées 
plus haut. On admet généralement que les facteurs 
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d'intervalle des électrons p sont liés par les rela- 
tions : | 


avec 
Fyuja = 20870 ; Fygo = 11338 ; Grl1,76) = 1,1480. 
(1 — 5), = 0,96 ; (4 — 5)po > 1. 


D'où : 
Gpsl2nes l 
Aile 8,882 
d'y Me 4] 
EN NT 


La validité de ces formules peut être vérifiée 
dans le cas de l’isotope 199 pour lequel les facteurs 
d'intervalle des termes de la configuration 6s6p 
ont été mesurés avec une:très grande précision. On 
obtient les équations suivantes : 


A(3P9) = , M . Epais— 302,9 


à 
A(8P9) = = (1,4159 a, + 0,9735 aps + 1,6106 ayayo 


| 
+ 2,2394 CS) = 491,9 
A(LP9) — :. (— 0,4159 a, + 4,0265 ap + 0,3894 Gp 


— 2,2394 ’») = —119,2. 


Avec les relations imposées entre les à, on arrive 
à un système d’équations incompatibles. On est 
obligé d’admettre que la valeur de a, croît nota- 
blement lorsqu'on passe du terme *P£ au terme 1P9 
ou, ce qui est plus vraisemblable, que le rap- 
port Gysye [pie St nettement inférieur à la valeur 
théorique 1 /8,5. Dans cette dernière hypothèse, on 
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est conduit à considérer comme trop grande, la 
valeur de Q déduite de 3P9 et beaucoup trop petite 
celle déduite de 1P9?. 

Une situation analogue se présente dans la confi- 
guration 6s6p du spectre d’arc de l’ytterbium. Les 
valeurs du moment quadrupolaire de l’isotope 173 
(I — 5/2) ont été déduites des termes 5P9 et 5P° 
par Schüler, Roig et Korsching [17] et des 
termes $P9 et 1P9 par Krebs et Nelkowski [18], à 
l’aide de la théorie de Breit et Wills. (Pour °P9, 
c, — 0,4742 et c, — 0,8804, valeurs comparables à 
celles du terme correspondant de Hg 7). Les valeurs 
de Q18 ainsi obtenues sont comparées dans le 
tableau IV à celles de Q?°1. 


TABLEAU IV 


3P, sp, 1, 
QI x 1024 SRI KSUISE TT 

K. N. 4,0 21 

Q201 x 1024 0,48 0,55 0,27 

Q1r3/Q01 FBtur 0S COU 


La variation de la valeur de Q en fonction du 
terme spectral considéré est donc tout à fait iden- 
tique pour les deux éléments. Comme c’est pour le 
terme #P£ que le calcul est le plus sûr. nous admet- 
trons en définitive comme valeur du moment qua- 
drupolaire de l’isotope 201 : 


Q = (0,50 + 0,05) x 10—?4 cm?. 


La correction de Sternheimer n’a pas été appli- 
quée. | 

Ce travail a été effectué sous la direction de 
M. le Pr P. Jacquinot, à qui nous voulons exprimer 
notre gratitude pour les utiles conseils qu'il n’a 
cessé de nous donner. Nous tenons également à 
remercier tout particulièrement MM. P. Giacomo et 
S. Roïzen qui ont préparé les couches diélectriques 
que nous avons utilisées. 


Manuscrit reçu le 6 décembre 1956. | 


Î 


_ 


[10] Scaüzer (H.) et Jones (E. G.), Z. Physik, 1932, 74, 
631. 

[11] Baise (J.), CHaBBAL (R.) et Jacquinor (P.), J. Phy- 
sique Rad., 1954, 15, 749. 

[12] McCzuwc (F.) et Hozmes (J. R.), Phys. Rev., 1955, 
100, 1792. 

[13] He (H.) et Jones (E. G.), Z. Physik, 1932, TT, 

4 


[14] Gounswrr (S.) et Bacner (R. F.), Phys. Reo., 1933, 
, 894. 
[15]' Procror (W.) et Yu (F.), Phys. Reo., 1949, 76, 1728 ; 
Phys. Rev., 1951, 81, 20. 
[16] IBreir (G.) et Wizzs (IL. A.), Phys. Rev., 1933, 44, 
470 


17] ScuÿLen (H.), Rouc (J.) et Konscwinc (H.), Z. Physik, 

1938, 111, 165 el 386. 

[18] Kress (K.) et Nezkowskr (H.), Z. Physik, 1956, 145, 
543. 


1 
1] 
| 


LE JOURNAL DE PHYSIQUE ET LE RADIUM 


TOME 18, MARS 1957, PAGE 201. 


LETTRES À LA RÉDACTION 


DIFFUSION ÉLASTIQUE DES PHOTONS 
PAR LE CUIVRE 


par J. B. Bercicarp, J. Mizcer et C. TzARA, 
Service de Physique Nucléaire, C. E. N., Saclay. 


La diffusion des photons par les noyaux dans la 
région de la résonance géante a été observée par une 
méthode d’activation [1] et par des méthodes de 


. spectrométrie des photons diffusés [2, 3, 4]. 


diffuseur _et délecteur 
de cuivre 


Chambre d'ionisefian à parois cpaisses 


Principe de la méthode. — Un spectre de frei- 


nage N(E, En) frappe une cible de cuivre dont la 


section eflicace différentielle de diffusion élastique 


est ee (E, 6). Un détecteur annulaire de cuivre 
s'active avec une section efficace 6,,(E) (fig. 1). L’acti- 
vation pour une énergie maxima Æn est donc : 


E j ni 
M(En) = À [NE En) SyrlE) cynlE) dE. 


Dans À rentrent les facteurs de géométrie, d’absorp- 
tion, de rendement de détection du rayonnement du 
noyau activé et d'intensité incidente, ainsi qu'un 
facteur dépendant de la distribution angulaire du 
rayonnement diffusé que nous choisissons en Î + cos? 0. 

Nous supposerons que cette distribution ne varie pas 


_ avec l’énergie. Nous négligeons la contribution des 


photons inélastiques que nous ne pouvons distinguer. 


Corrections. - Les diffusions Rayleigh et Thomson 
ont une section eflicace o < 5.10-%%em°?. 

Les photons diffusés par effet Compton simple n’ont 
pas une énergie suffisante pour activer Cu. La diffusion 


Utilisant pour la détection de l’activité induite par 
les photons diffusés un cristal de INa de 10,16 cm 
X 10,16 em nous avons pu obtenir des résultats quan- 
titatifs sur la diffusion des photons par le cuivre 
(mélange isotopique naturel). Ce genre de mesure, 
quoique indirect, peut être utilisé devant des sources 
intenses mais pulsées de rayonnement de freinage dans 
des cas où la spectrométrie par scintillation échoue, à 
cause des empilements de photons diffusés de basse 
énergie. 


Paraffine 


ÉREer 
20 cm 


Cible de freinage 


Emplscement de lecren de Plomb 


F1G. 1. — Dispositif expérimental. 


Acfivafion 
en unife 
arbitraire 


G 


N 


= 


See 
15 20 25 


F1. 2. —- Courbe d’activation. 
Cercles vides : activations sans correction de neutrons 
rapides. 
Cercles pleins : activations corrigées pour les neutrons 
rapides. 
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Compton multiple a une section efficace négligeable. 
La section efficace d'émission de photons de freinage 
suffisamment énergiques par les électrons de la cascade 
photons-électrons dans le diffuseur est de l’ordre 
de 2.10—?° cm?. Nous négligeons ces corrections devant 
la section efficace nucléaire 6,4, = 10? cm?. 

Les neutrons rapides émis par le collimateur de 
plomb peuvent activer le détecteur de cuivre. L’inter- 
position d’un écran de plomb dans le faisceau réduit 
l’intensité de photons dans un rapport connu. Le résidu 
d'activité est dû aux neutrons rapides formés dans le 
collimateur et dans l’écran. Nous en déduisons la limite 
supérieure de la correction de neutrons rapides. 


Résultats. — Nous obtenons deux courbes d’acti- 
vation pour Ex — 18, 20, 21, 22 MeV + 100 keV 
(fig. 2). L’une correspond à l’activité sans déduction 
de l’activation par réaction (n, 2n), l’autre à celle une 
fois déduite cette activation. Pour déduire de l’acti- 
vation 6,,(£), nous prenons une courbe moyenne lissée. 
De la connaissance;de 6,,(Æ) (cf. 5) et de N(E, En) nous 


7 


| Org 10 “cm 


15 + 


++ 


F1c. 3. — Courbe 1 : section efficace. 
En grisé : zone d’imprécision. 
Courbe 2 : section efficace déduite de la relation de 
dispersion. 
Points 3 : résultats de Fuller et Hayward. 


€ 


déduisons la section efficace o,4(E) (fig. 3). La zone 
grisée indique l’imprécision provenant de la statistique. 
Une incertitude de + 10 % provenant de l’évaluation 
des facteurs correctifs s’y ajoute. 


Diseussion. — Notre mesure est en accord général 
avec celle de E. G. Fuller et E. Hayward [4] et le 
résultat déduit de la relation de dispersion de Kramers- 
Heisenberg, sauf du côté des énergies élevées, où elle 
est plus faible. 

Nous avons supposé comme Fuller et Hayward que 
la distribution angulaire est en À + cos? 0, caracté- 
ristique de la diffusion cohérente par interaction dipo- 
laire. Dans nos mesures, où l’angle solide de détection 
est de 0,75 X 4x, o,, est peu sensible à une variation 
de distribution angulaire. Le problème de la distri- 
bution angulaire reste à approfondir. Finalement nous 
n'avons aucune donnée sur la part de diffusion iné- 
lastique qui, par comparaison avec les résultats pré- 
cités, semble très petite. 


Lettre reçue le 12 janvier 1957. 
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ÉMISSION ÉLECTRONIQUE SECONDAIRE 
PROVOQUÉE PAR IMPACT D'IONS SODIUM 
SUR DES ÉCHANTILLONS MÉTALLIQUES 

CONTRAINTS 


par MM. R. GourtTe et C. GUILLAUD, 


Laboratoire d’Optique électronique ; 
et de Physique du Métal, Université de Lyon. 


Les publications précédentes de notre laboratoire sur 
l'émission électronique secondaire [1] ou photo- 
électrique [2] de surfaces métalliques soumises à des 


Alumine frittée 
24 


Fi. 1. —— Dispositif expérimental. 
F : filament.— C: collecteur. — W : wehnelt. 
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déformations mécaniques ont montré que le coefficient 


d'émission croît en général avec la contrainte, la valeur 
de cette exaltation pouvant atteindre 10 à 20 %. 
Nous avons voulu rechercher si un effet analogue était 
observable dans le cas où l'émission électronique secon- 
daire est provoquée par un bombardement primaire 


ionique. 


Sources d’ions. — La source utilisée est un émetteur 
d’ions solide. Il s’agit d’un mélange de carbonate de 


| sodium, d’alumine et de silice, fondu vers 1 650° pour 
| réaliser le seltriple Na,0, AL,0;, 2810, [3]. C’est princi- 


palement l’ion sodium qui est émis [4]. Le produit 
refroidi est broyé finement. Nous en déposons une 
petite quantité dispersée dans du collodion, à l’extré- 
mité d’un filament de tungstène en V qui sert d’anode. 


Un préchauffage sous vide élevé permet de stabiliser 


l'émission. 


Dispositif expérimental. — L’anode chaude est 
placée dans un wehnelt, convenablement polarisé, pour 
concentrer les ions émis sur l’échantillon. Le collecteur 
des électrons secondaires émis est protégé du bombar- 
dement ionique primaire par un cylindre relié à la 
dispositif (fig. 1) a été utilisé par 
G. Couchet [3]. L’échantillon constitué par un ruban 
métallique est mis sous traction jusqu’à rupture grâce 
au dispositif décrit précédemment [11]. L’ensemble est 
placé sous la cloche d’un bloc à évaporation donnant 
un vide de 2 X 105 mm de Hg [5]. Le dispositif de 
mesure électrique est représenté sur la figure 2. 

Le filament et la batterie S, assurant son chauffage 
sont portés à la haute tension ainsi que le wehnelt et 
ses piles de polarisation S.. Le réglage du chauffage 


F1G. 2.— Schéma du montage électrique. 


F : filament. — W : wehnelt. — C: collecteur. — G : écran. 
E : échantillon. — P, : potentiomètre de chauffage. — I : 
intensité chauflage.— $, : batterie 6 V.—P,: polarisa- 
tion du wehnelt. — $, : pile 180 V. — R, : résistance 
de sécurité. — R, : résitance de débit. — KV : mesure 
de la H. T.—$S, : pile 150 V. — P, : polarisation du 
collecteur. — G, : mesure du courant secondaire. — G, : 
mesure de la somme des courants primaire et secondaire. 
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du filament s'effectue à l’aide de la résistance 
variable P;. La polarisation du wehnelt est ajustable 
avec le potentiomètre P.. R, est une résistance de 
protection en série avec la haute tension. Le collecteur 
est relié à sa polarisation S;, également ajustable avec 
le potentiomètre P;, par l'intermédiaire du micro- 
ampéremètre G:. G, mesure le courant secondaire £. 
L’échantillon est relié à la masse par le micro- 
ampèremètre G,. G, mesure la somme des courants : 
ionique primaire et électronique secondaire, soit 1 + 7. 


Résultats expérimentaux. —- Les valeurs des coef- 
ficients d'émission secondaire y = 1/1 sont données 
pour une tension d’accélération des ions primaires de 
1 500 volts. 

1) Pour un échantillon de nickel préalablement 
recuit à À 3000 au four à hydrogène, nous trouvons : 


DOC ASS A PIE T= 166 K 107$ 2AÀ 


ce qui conduit à un coefficient d’émission secon- 
daire y — 0,344. Cette valeur est en accord avec les 
résultats publiés par Couchet [31et concernant le même 
métal, les ions primaires étant alors des ions lithium. 
Un effort de traction provoquant la rupture provoque 
une augmentation sensible du coefficient y qui devient 
yr = 0,362. L'augmentation relative est de 6 %. 
2) Avec une cible d’argent, nous avons : 


1 — 35,5 x 108 À I +i—3150 x 10-8 A. 


Le coefficient d'émission secondaire est donc 
y = 0,31. Après traction yr = 0,327. L'augmentation 
relative est de 5,5 %. 

3) Dans le cas d’une cible d’or 


USD D AM OS EEA REbET =MR oc MTORSUAE 


Le coefficient y de l’or est y = 0,35. Après traction 
le coefficient d'émission secondaire devient Yr = 0,36. 
L'augmentation relative qui est de 3 %, est plus faible 
que dans les deux cas précédents. 


Conelusion. — Ces résultats sont évidemment très 
incomplets, bien que régulièrement reproductibles. Ils 
téndraient néanmoins à confirmer que l’accroissement 
d'émission provoqué par la contrainte se produit éga- 
lement si le bombardement primaire est de nature 
ionique. | 

L'importance de cet accroissement (3 % à 6 %) reste 
dans l’ordre de grandeur de ceux que nous avons 
obtenus avec les photoélectrons [2]. 

Malheureusement, les ions alcalins utilisés polluent 
très vite la surface et l’influence de cette pollution 
reste indéterminée. Pour éviter ce grave inconvénient, 
nous comptons reprendre ces essais avec des ions 
gazeux légers qui offrent au contraire l’avantage de 
décaper les surfaces qu’ils bombardent. 

Lettre reçue le 20 janvier 1957. 
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SUR LE PRINCIPE DE MACH 
ET QUELQUES QUESTIONS CONNEXES 


Par O. Cosra DE BEAUREGARD, 
Institut Henri-Poincaré, Paris. 


Une intéressante suggestion de D. Park [1] est à 
l’origine de remarques que nous avons exposées ici 
même [2] ainsi que d’un modèle de cosmos expansif 
récemment proposé par nous [3] avec une applica- 
tion [4] aux idées de Dirac sur les constantes univer- 
selles [5]. Ë 

Toutefois, nous n’avions pas su bien mettre en évi- 
dence le point essentiel, et c’est ce que nous nous pro- 
posons de faire ici. 

L'idée de D. Park est d'interpréter l'hypothèse bien 
connue de Mach relative à l’inertie, en disant que la 
constante universelle / de la loi de Galilée-Newton 
> — 

F — Im représente en réalité la valeur du fond non 


nul du potentiel de gravitation ; introduisant donc la 
masse 7 du cosmos (supposée finie), une longueur L 
caractéristique des dimensions globales de l’espace 
et la constante de Newton G ou, mieux, la constante 
d'Einstein y — 8x Ge, D. Park propose d'écrire 


T= xMIZ ; (1) 


le second membre est bien (comme on le sait) un 
nombre pur, et il doit se trouver rendu égal à 1 du fait 
du choix même des unités classiques. 

Une confirmation frappante de l’idée de D. Park 
est fournie par la théorie cosmologique d’Einstein, 
où l’on a, À désignant le rayon de la sphère spatiale, 


xXMIA = &r?; (2) 


des relations ne différant que par la valeur du second 
membre valent dans plusieurs modèles de cosmos 
expansifs, et notamment celui récemment proposé par 
nous (3). Ceci incite à un examen plus serré de la 
question. 

La constante universelle de Ja formule F — 1mY 
est certainement un scalaire d’espace-tmps, car cette 
formule admet la transcription 4-dimensionnelle bien 
connue FAdr — ImoduÀ. 

Par ailleurs, la formule (2) est conséquence très 
directe de l’équation d’Einstein contractée 

—R=yxT, ou —y%;T[]R—1, (3) 
comme on le voit au moyen des relations (où 
2r?A8 est le volume de l’hypersphère spatiale) : 


M = 27? A6, AR = —%, 
pole = rapport simple. 


T = Po 
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C’est donc l’équation d’Einstein contractée (3) qui 
assure automatiquement la satisfaction du principe de 
Park : le nombre — yTR-1, homogène à un potentiel 
gravifique, se trouve, dans le cas d’un cosmos sphérique 
homogène, égal à l’expression 

XM,187? À =1, avec M, = 27? A8 T(A). (4) 

Avant de changer de sujet, faisons encore deux 
remarques. Au second membre de la formule d’Eins- 
tein ei 

Rav. Ce ou — y Tu 5 
9 5 X (5) 


figure la somme de toutes les densités d’impulsion- 
énergie inertiques, potentielles, etc.…, sauf une: la 
densité d'impulsion-énergie potentielle de gravitation. 
Or, l’énergie potentielle de gravitation figurant dans 
la théorie classique doit bien apparaître quelque part 
enrelativité. Jl est tout naturel d'admettre que l’expres- 
sion 


R 
m2: AU == plu 
(2 2 8 ) 


représente la densité d’impulsion-énergie potentielle de 
gravitation ; l'équation (5) d’'Einstein exprime alors 
en particulier que la somme de toutes les densités 
d’impulsion-énergie est nulle. 

L = vD/èv désignant maintenant la constante de 
récession de Hubble, Jordan [6] a remarqué que, 
d’après les valeurs expérimentales des trois facteurs 
intervenant, le nombre ypZ? semble avoir une valeur 
de l’ordre de Punité. Or, théoriquement, il en doit bien 
être ainsi en conséquence de la formule d’Einstein 
contractée (3) ; en effet, le rapport p0/0 = T/Ta est 
toujours un nombre petit, et 1l en est de même du 
rapport L/A dans la plupart des modèles de cosmos 
expansifs. En bref, la remarque de Jordan peut être 
considérée comme l'expression d’une vérification (fort peu 
précise il est vrai) de la théorie de la relativité générale. 


Manuscrit reçu le 5 février 1957. 
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Structure hyperfine dans le spectre de l’indium II, par 
D. A. Jackson, Laboratoire Aimé Cotton, CG. N.R. $. 
Bellevue (S et O) (J. Physique Radium., 1957, 18, 145). — 
Les mesures de structure hyperfine des raies de In 77 
émises dans une lampe à cathode creuse chaude, faites par 
Paschen et Campbell avec un réseau concave donnaient 
des rapports d’intervalles montrant de grandes divergences 
avec la règle des intervalles. Dans le présent travail, la struc- 
ture hyperfine du multiplet 5% Di,9,1 — 4% Fiso (4 638 à 
4 685 À) du même spectre a été redéterminée en prenant 
comme source une cathode creuse refroidie à l’eau et comme 
instrument résolvant un étalon de Fabry-Perot, traité avec 
des couches multiples diélectriques. Grâce à la plus faible 
largeur des raies émises par la source et à la plus grande résol- 
vance effective de l’étalon, une résolution plus complète de 
la structure a été obtenue, rendant possible une évalua- 
tion nettement plus précise des rapports d’intervalles dans 
les termes mis en jeu. Bien que ces nouvelles mesures ren- 
dent nécessaires quelques changements dans la corrélation 
des transitions et dans la structure hyperfine des compo- 
santes, la nouvelle analyse est, dans l’ensemble, en accord 
avec l’analyse originale de Paschen et Campbell, et les 
grands écarts à la règle des intervalles sont confirmés ; 
la cause de ces divergences est discutée. 
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Résolution de l’équation de Schrôdinger des atomes à 
États S symé- 
triques, par G. Munscay et P. PLuvinace, Institut de 
Physique de la Faculté des Sciences de Strasbourg (J. Phy- 
sique Rad., 1957, 18, 157). — Il est possible de définir une 
triple infinité de fonctions de base telles que l’équation de 
Schrôdinger des atomes à deux électrons soit résolue rigou- 
reusement par des développements suivant ces fonctions. 
Les coefficients des développements sont solutions d’un 
système linéaire et homogène de rang infini mais dont 
chaque équation se compose d’un nombre fini de termes. 
En égalant à zéro le déterminant, on trouve les racines 
carrées des valeurs propres de l’énergie changées de signe. 


Nouvelles mesures des rapports gyromagnétiques du fer 
et du nickel, par A.J. P. Meyer etS. Brown, Labora- 
toire Pierre-Weiss, Institut de Physique, Strasbourg, (J. 
Physique Rad., 4957, 18, 161). — Les mesures décrites ont 
été entreprises en vue de découvrirl’origine du manque de 
concordance entre les résultats expérimentaux antérieurs. 
On a pu montrer qu’une erreur systématique faussait les 


Pre "+ mec % 
rapports gyromagnétiques voisins de 1,00 (7) mesurés 


par.effet Einstein-de Haas. Après correction, les résultats 
expérimentaux des différents auteurs entrent en bon accord 
entre eux. On décrit ensuite le nouvel appareil de mesure 
de l’effet Einstein de Haas à résonance forcée et à bobine 
d’aimantation mobile. Les expériences donnent pour le fer 
o (5) — 1,035 —+- 0,004 en accord avecles résultats d’autres 


auteurs. Pour le nickel on trouve 9 (5) = 1,078 + 0,005 ; 
valeur proche de celles de Scott. On montre queles valeurs 
voisines de 1,04 données antérieurement pour le nickel sont 
erronées pour des raisons diverses. 


Une forme spinorielle des équations de l’électromagné- 
tisme, par G. Bopiou, Faculté des Sciences de Marseille. 
(J. Physique Rad., 1957, 18, 169). — II s’agit d’une descrip- 
tion d’un champ électromagnétique général, etd’une trans- 
cription des équations de Maxwell-Lorentz, en termes de 
spineurs du premier rang, dans l’espacetemps : Tout bivec- 
teur réel est le dual de la somme des densités de moment 
électromagnétique de deux spineurs, confondus avec leurs 
« lacunes » respectives. La solution générale des équations 
de Maxwell-Lorentz est de cette forme, ses deux spineurs 
générateurs étant solutions d’équations convenables. Dans 
le vide, et pour une classe étendue de champs, ces deux 
spineurs sont des solutions, convenablement liées, de l’équa- 
tion de Dirac à masse nulle. (Ces résultats ont été énoncés, 
sans preuves, aux C. R. Acad. Sc., 29 octobre 1956, 248, 
p. 1287.) 
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New Measurements of the Gyromagnetic Ratios of Iron 
and Nikel, by Anpre J. P. Meyer and SHELDON BROWN, 
Pierre Weiss Laboratory, [Institute of Physics, Strasbourg, 
France (J. Physique Rad., 1957, 18, 161). — Gyromagnetic 
measurements undertaken to reveal the source of discor- 
dance between previous experimental results are des- 
cribed. Itis shown that a systematic error was responsible 
for the incorrect gyromagnetic ratio of approximately 
1.00(mc/e) found in some previous measurements by the 
Einstein-de Haas method. After correction of these values 
the results of different authors are found to be in good 
agreement. An apparatus with a moving coil for the 
resonant measurement of the Einstein-de Haas effect is 


described.The new experiments give for iron p (2) = 
1.035 + 0.004 in agreement with the results of other 
= 1.078 + 0.005, 


a value close to that found by Scott. It is shown that 
the values in the neighborhood of 1.04 found previously 
for nickel are erroneous for various reasons. 


authors. For nickel one finds p (£) 


A spinor form of the electromagnetice equations, by 
G. Bopiou, Faculty of Science, Marseille (J. Physique Rad., 
1957, 18, 169). This paper contains a description of the 
general electromagnetic field and a transcription of the 
Maxwell-Lorentz equations in terms of first-rank spinors 
in space-time : Any real bivector is the dual (or ortho- 
complementary) of the sum of the electromagnetic densities 
determined by two spinors compounded with their res- 
pective ‘‘ lacunal ”. The general solution of the Maxwell- 
Lorentz equations is of this form, its two generating spinors 
being solutions of suitable equations. In vacuum, and 
for a large classe of fields, these two spinors are solutions, 
suitably interrelated, of the Dirac equation for zero-mass 
corpuscules. (These results were stated, without proofs, 
in C. R. Acad. Sc., 29 octobre, 1956, 248, 1287.) 


Hyperfine Structure in the Spectrum of Indium Il, 
by D. A. Jackson, Laboratoire Aimé-Cotton, G. N. R. S. 
Bellevue (S.-et-O.) (J. Physique Rad., 1957, 18, 145). — 
Measurements of the hyperfine structure of lines in the 
spectrum of In IT emitted by a hot hollow cathode discharge 
made by Paschen and Campbell with a concave grating, 
showed interval ratios with large divergences from the 
interval rule ; in the work described below the hyperfine 
structure of the multiplet 5 3D3,9,1 — 4 %F4,3,2 (4 638 to 
4 685 À) of the same spectrum has been re-determined 
with a water cooled hollow cathode as light source and a 
Fabry-Perot etalon, with multilayer dielectric reflecting 
surfaces as the instrument of high resolving power. Due 
to the smaller width of the lines emitted by the cooled 
hollow cathode and the greater effective resolving power 
of the etalon, more complete resolution of the structure 
was effected so that a substantially more accurate evalua- 
tion of the interval ratios in the terms 5 #D;,,,, and 
&3Fy,e,9 and 4!F; was possible. Although the new 
measurements make necessary some changes in the corre- 
lation of the” transitions and the observed hyperfine 
structure components, the new analysis is in the main 
in agreement with the original analysis of Paschen and 
Campbell, and the large departures from the interval rule 
are confirmed ; the cause of these divergences is discussed. 
The indium used in the hollow cathode was natural indium 
containing 96 % 45 In and 4 % 5 In ; no indication 
of resolved components due to the rarer isotope was found. 


Solution of the Schrôdinger equation of two-electron 


atoms. II. Rigorous method. Symmetric states, by G. 
Muxscay and P. PLuviNAGE, Institut de Physique de 
la Faculté des Sciences de Strasbourg (J. Physique Rad., 
1957, 18, 157). The Schrôdinger equation of two-electron 
atoms is rigorously solved by a development in series of 
orthogonal functions. The method is here explained in 
thesimplest case of the Ssymmetricstates. A determinant 
can be constructed, the roots of which are (— Æ) 42, In 
contrast to the variation method, there is no variation of 
argument and the number of terms in each row or column 
is finite. 
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Réversibilité et irréversibilité en résonance nucléaire. I. 
héorie de la relaxation nucléaire dans les liquides, par 
Joseph Seine, Laboratoire d’Électronique et de Radio- 
électricité, avenue du Général-Leclerc, Fontenay-aux-Roses 
(J. Physique Rad., 1957, 18, 173). — On démontre à partir 
de l’équation de Schrôdinger réversible par rapport au 
_ temps, et sans faire intervenir l’hypothèse des phases répar- 
ties au hasard, l’irréversibilité de l’évolution d’un système 
de spins nucléaires couplés à un thermostat. L’interaction 
spin-spin est supposée suffisamment faible pour que le 
comportement d’un spin donné 7 puisse être considéré 
comme indépendant de celui de ses voisins. Le thermostat 
.est considéré comme un système classique (et non quan- 
, tique), ce qui est légitimé par le fait que les énergies mises 
en jeu sont faibles devant AT. L’énergie de couplage spin- 
thermostat 2 VW(t) est considérée comme une fonction 


ÿ 

aléatoire stationnaire du temps, elle induit des transitions 
> > 
entre les états propres de À = — y% H, E I; énergie des 
ÿ î 

— 
spins dans le champ magnétique constant H,. Les spins 
sont écartés de l’équilibre thermique avec le thermostat 


par un champ magnétique de radiofréquence. Pour qu'il y 
ait relaxation à partir de l’état de non-équilibre ainsi créé, 


… les spins doivent vérifier certaines conditions. Des condi- 


tions suffisantes sont que les spins soient en nombre infini 
. (pratiquement en nombre très grand) et que les fonctions 

d’onde Y@{(0) des divers spins à l’instant initial de non- 
équilibre forment un ensemble dense dans l’espace des 
fonctions d’onde des spins individuels. La théorie de la 
relaxation est alors basée sur un développement de l’opé- 
rateur d’évolution en série de puissances de VW{r). On 
introduit une matrice de relaxation o{t), jouant ici un rôle 
analogue aux fonctions de distribution de Boltzmann. Le 
calcul de o{t) au deuxième ordre en VW sert d'introduction 
au calcul rigoureux de ot) à tous les ordres en VW) néces- 
saire pour n’avoir à invoquer la condition de densité des 4 
qu’à l'instant initial : — 0. Ce calcul rigoureux montre que 
les processus les plus généraux qui contribuent à la rela- 
xation consistent en une succession de processus du 
deuxième ordre en VU), au cours desquels la perturbation VO 
agit par deux fois dans des intervalles de temps inférieurs 


2: 
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au temps rt, de corrélation des V{(1). On parvient ainsi à 
former une équation de Boltzmann vérifiée par o(1) à tous 
les ordres en V( et qui fournit une description complète 
de la relaxation énergétique et transversale. La déduction 
ne fait appel qu’aux seules lois de la mécanique quantique, 
on n'introduit pas d'arguments statistiques, dont le rôle 
dans la théorie des phénomènes de relaxation est stric- 
tement circonscrit. La théorie est applicable à la résonance 
nucléaire magnétique dans les liquides et les gaz, et, avec 
certaines restrictions, à la résonance dans les cristaux molé- 
culaires, les métaux etc. On montre enfin que le problème 
de la détermination de la forme des raies de résonance 
nucléaire est équivalent à celui de la relaxation transversale 
d’une grandeur ayant pour règle de sélection Am = + 1. 


Structures hyperfines de raies du spectre d’are du mer- 
eure et moment quadrupolaire de 21Bg, par J. BzaIse et 
H. CHaAnTRrez, Laboratoire Aimé Cotton, CG. N.R. S., Bel- 
levue (S-O). (J. Physique et Rad., 1957, 18, 193). — Les 
structures hyperfines des transitions 6s6p — 6s7s et 6s6p 
— 6s7d de Hg Z ont été mesurées avec précision. Les 
termes 6s 7d1D,et 6s7d D, présentent des anomalies dues à 
ce que la structure hyperfine est du même ordre de 
grandeur que la structure fine de la configuration 6s7d 
Le rapport des moments magnétiques des isotopes impairs 
est : ue 20174 189 — —— 1,110 0,001. La valeur du mo- 
ment quadrupolaire de 22Hg varie avec le terme spectral 
utilisé pour son calcul. A partir des termes 6s 6p 5P5, Pet 
1 P9 on obtient respectivement 0,48, 0,55 et 0,27.10—?4cm?. 
Une explication de ces divergences est proposée. On est 
conduit à admettre Q — (0,50 + 0,05) X 1024 cm?. 
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intervals of the size of %, the correlation time of (5). 
A system of differential equations (Boltzmann équations) 
which must be obeyed by ot) is formed ; these equations 
give a complete specification of the relaxation process. 
We do not introduce statistical arguments, the demons-, 
tration being accomplished only with the help of the 
laws of quantum mechanics. The theory can be applied 
to nuclear magnetic resonance in liquids and gases and, 
with some restrictions, to nuclear resonance in molecular 
crystals and metals. It is shown that the determination 
of the absorption spectrum is equivalent to the problem 
of the relaxation of a quantity with selection rule 
Am= +1. 


Hyperfine Structure in the Are Spectrum of Mercury 
and Quadripole Moment of °° Hg, by J. BLaise and 
H. CHANTREL, Laboratoire Aimé-Cotton, G. N. R. S$., 
Bellevue (S.-et-O.) (J. Physique Rad., 1957, 18, 188). — The 
hyperfine structure of the lines corresponding to the 
transitions 6s6p — 6s7s and 6s6p — 6s7d of Hg I have 
been measured accurately. The 6s7d 1D, and 6s7d D, 
terms show anomalies due to the fact that the hyperfine 
structure is of the same order of magnitude as the multiplet 
structure of the 6s7d configuration. The ratio of the 
magnetic moments is: u?01/ul% — —_ 1,110 + 0.001. 
The values of the quadrupole moment of 24Hg — He 
calculated from the 6s6p PS, P9 and !P, terms are 0.48, 
0.55 and 0.27 X 10V!1 cm!, respectively. A possible 
explanation of these discrepancies is given. The best 
value is : Q — (0.50 + 0.05). 10—24 cm2. 
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Reversibility and irreversibility in nuclear resonance. 
I. Theory of nuclear relaxation in liquids, by J. SeiIDEN, 
Laboratoire d’Électronique et de Radioélectricité, Fontenay- 
aux-Roses. (J. Physique Rad., 1957, 18, 173). Starting. 
with the Schrôdinger equation reversible relative to time, 
we demonstrate the irreversibility of the evolution of am 
system of nuclear spins in thermal contact with à heat 
reservoir. The spin-spin interaction is supposed to be 
sufficiently weak for the behaviour of each spin 7 to bem 
considered as independent of the behaviour of its neigh= 
bours. The heat reservoir is considered as a classical 
system ; this is justified because the energies are € KT 
The interaction energy È V@(:) of the spins with the heat 


3 + 
reservoir is considered as a random function of time ; this 
energy causes transitions to take place between them 
stationary states of À — — y# H, Z L;, the energy of ‘# 


the spins in the magnetic field H,. The spins must obey 
certain conditions for relaxation to take place ; sufficientM 
conditions are that there are an infinity of spins, and that 
the wave functions W{(0) of the spins at the initial time 
form a compact ensemble in the space of the functions 
FO. The theory of relaxation is based on a development 
of the evolution operator in powers of the perturbations 
V@(1). We introduce a relaxation operator o(t), similars 
to the distribution functions of Boltzmann. It is neces-M 
sary to calculate oft) to all orders in V@ in order to use! À 
the hypothesis of compactness of the W((0) at the initial 
time only. This rigorous calculation shows that the only. 
processes contributing to the relaxation are those of second 
order in V@. The perturbation V® acts twice in time 
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